Hi,
I have adapted the navier_system.C example to solve the Euler equations. I
have tested it for a sample problem (flow through a channel without
disturbances, so that the solution should remain the same for all time steps)
without AMR, and things seem to work fine. So, if I initialize th
Hi,
In order to fix this, I ended up adding the following change
to UnsteadySolver::reinit().
Currently, it returns after resizing the vector. I copied the portion
after the if block in UnsteadySolver::advance_timestep after the resizing
part in reinit(). This forces it to copy the syste
One more observation:
UnsteadySolver inherits from TimeSolver, but both of them define a member
variable old_local_nonlinear_solution.
Since the UnsteadySolver operates and resizes this vector, it may be best
to remove it from TimeSolver to avoid confusion.
Manav
On Fri, Mar 1, 2013 at 1:46
Ok, so I looked at the code, and here is my assessment.
When using AMR, following the refine_and_coarsen_elements the
equation_systems.reinit is called.
This prolongs the vectors of the system, and reinitializes the values (not
just size) of the system vectors. However, it only resizes the
old_lo
