#NEWLINE# def linearize( self , params , unknowns , resids ) : #NEWLINE# #INDENT# ' Jacobian for our paraboloid.' #NEWLINE# x = params [ 'x' ] #NEWLINE# y = params [ 'y' ] #NEWLINE# J = { #NEWLINE# #INDENT# #NEWLINE# #UNINDENT# } #NEWLINE# J [ ( 'f_xy' , 'x' ) ] = ( ( ( 2.0 * x ) - 6.0 ) + J ) #NEWLINE# J [ ( 'f_xy' , 'y' ) ] = ( ( ( 2.0 * y ) + 8.0 ) + x ) #NEWLINE# return J