Changed how pardiso is called again. Now calling each time with phase =
13 (analyze factorize and solve) then phase=-1, resetting pt to 0 after
each call with phase=-1 So analysis is redone buy miraculously, it seems
to run super fast this way. Also added timing of full program run (in
main program) and full FRONT routine (from ARJUNA_2D). Timing done with
OMP_GET_WTIME() and written to ArjunAir.out
Fixed undefined variable in OpenMP timing of whole program run
Added extra output for inversion diagnostics. Now write,
delpar, and model parameters at first iteration to file.
Finished OpenMP matrix assembly using local variables for boundary
columns. Tested it and it doesn't work. Will continue debugging
sometime.
OpenMP parallelized matrix assembly working now. I wish I could get
around the reduction slowdown but that seems to be unavoidable.