You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Convert the original main program into a function that returns a result for comparison to the new code. The result would probably be declared something like double precision, allocatable :: positions(:,:,:), which would be allocated to shape [ncells, npositions, 3], where 3 is the number of spatial directions and where positions(:,:,1) holds x values, positions(:,:,2) holds y values, and positions(:,:,3) holds z values.
TBD: we can discuss in an upcoming call how best to compare the function result with that of the new code.
double precision, allocatable :: positions(:,:,:), which would be allocated to shape [ncells, npositions, 3], where 3 is the number of spatial directions and wherepositions(:,:,1)holdsxvalues,positions(:,:,2)holdsyvalues, andpositions(:,:,3)holdszvalues.