HiFiLES / HiFiLES-solver

High Fidelity Large Eddy Simulation Solver
Other
171 stars 131 forks source link

Solution not evolving #133

Closed burgreen closed 7 years ago

burgreen commented 7 years ago

I am running the euler cylinder test case (and several other test cases). For each case, everything appears to be executing fine, but the solutions for all test cases are not evolving. I am running single processor, no MPI, no GPU.

Here is my typical output. What could be the problem?

$ ../../../a input_cylinderinv ____ . | | | | | | | || | | | | | / | | || | | | | | | | | | | | | (----` | | | | | | | | |__| | | | | \ \ | | | | | | | | | | |`----.| |__.----) | || || || || |**| |**|||/

Aerospace Computing Laboratory (Stanford University)

---------------------- Non-dimensionalization ---------------------

----------------------- Mesh Preprocessing ------------------------ reading connectivity ... done reading connectivity reading vertices done reading vertices Setting up mesh connectivity Done setting up mesh connectivity reading boundary conditions done reading boundary conditions

---------------- Flux Reconstruction Preprocessing ---------------- initializing elements tris Initializing tris quads tets pris hexas done initializing elements setting elements shape ... done. pre-computing nodal shape-basis functions ... done. setting element transforms ... at solution points 0.00% 9.94% 19.89% 29.83% 39.78% 49.72% 59.66% 69.61% 79.55% 89.50% 99.44% at flux points 0.00% 9.94% 19.89% 29.83% 39.78% 49.72% 59.66% 69.61% 79.55% 89.50% 99.44% done. initializing grid velocity to 0 ... done. setting element transforms at interface cubpts ... done. setting element transforms at volume cubpts ... Setting initial conditions... Writing Tecplot file number 0 on rank 0 Done writing Tecplot file number 0 ....

Iter Res[Rho] Res[RhoVelx] Res[RhoVely] Res[RhoE] Fx_Total Fy_Total 1 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 100 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 200 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 300 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 400 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 500 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 600 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 700 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 800 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 900 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 1000 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 Writing Tecplot file number 1000 on rank 0 Done writing Tecplot file number 1000 .... Execution time= 4.750000 s

burgreen commented 7 years ago

I rolled my own compile script. Turns out that not including a compilation flag for one of the BLAS define options was the issue.