Removed original AA solver code. added call to mkl_free_buffers() at end
of each run through the frequency loop. We'll see if that gets rid of
the memory leakage.
Changed pardiso calling. Metis ordering computed in phase 11 is now
stored in array perm, and then later used as a user supplied ordering.
Pardiso called with phase -1 after each fact./solve call to clear
internal memory.
commented out reseting of pt after calling pardiso with phase -1
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.