#NEWLINE# def test_extrusion( self , backend , elements , dat_coords , dat_field , coords_map , field_map ) : #NEWLINE# #INDENT# g = op2 . Global ( 1 , data = 0.0 , name = 'g' ) #NEWLINE# mass = op2 . Kernel ( '\nvoid comp_vol(double A[1], double *x[], double *y[])\n{\n    double abs = x[0][0]*(x[2][1]-x[4][1])+x[2][0]*(x[4][1]-x[0][1])+x[4][0]*(x[0][1]-x[2][1]);\n    if (abs < 0)\n      abs = abs * (-1.0);\n    A[0]+=0.5*abs*0.1 * y[0][0];\n}' , 'comp_vol' ) #NEWLINE# op2 . par_loop ( mass , elements , g ( op2 . INC ) , dat_coords ( op2 . READ , coords_map ) , dat_field ( op2 . READ , field_map ) ) #NEWLINE# assert ( int ( g . data [ 0 ] ) == int ( ( ( ( layers - 1 ) * 0.1 ) * ( nelems / 2 ) ) ) )