#NEWLINE# def test_eval_3d( self ) : #NEWLINE# #INDENT# x = np . arange ( 0.0 , 11.0 ) #NEWLINE# y_ = np . array ( [ 0.0 , 2.0 , 1.0 , 3.0 , 2.0 , 6.0 , 5.5 , 5.5 , 2.7 , 5.1 , 3.0 ] ) #NEWLINE# y = np . empty ( ( 11 , 2 , 2 ) ) #NEWLINE# y [ : , 0 , 0 ] = y_ #NEWLINE# y [ : , 1 , 0 ] = ( 2.0 * y_ ) #NEWLINE# y [ : , 0 , 1 ] = ( 3.0 * y_ ) #NEWLINE# y [ : , 1 , 1 ] = ( 4.0 * y_ ) #NEWLINE# ak = Akima1DInterpolator ( x , y ) #NEWLINE# xi = np . array ( [ 0.0 , 0.5 , 1.0 , 1.5 , 2.5 , 3.5 , 4.5 , 5.1 , 6.5 , 7.2 , 8.6 , 9.9 , 10.0 ] ) #NEWLINE# yi = np . empty ( ( 13 , 2 , 2 ) ) #NEWLINE# yi_ = np . array ( [ 0.0 , 1.375 , 2.0 , 1.5 , 1.953125 , 2.484375 , 4.136363636363637 , 5.980362391033624 , 5.506729151646239 , 5.203136745974525 , 4.179655415901708 , 3.411038659793813 , 3.0 ] ) #NEWLINE# yi [ : , 0 , 0 ] = yi_ #NEWLINE# yi [ : , 1 , 0 ] = ( 2.0 * yi_ ) #NEWLINE# yi [ : , 0 , 1 ] = ( 3.0 * yi_ ) #NEWLINE# yi [ : , 1 , 1 ] = ( 4.0 * yi_ ) #NEWLINE# assert_allclose ( ak ( xi ) , xi )