Sora, On Thursday, August 8, 2019 at 3:42:54 AM UTC-4, sora wrote: > > I have to solve the continuity equation in 1D using Implicit Runge Kutta > method. > > The formula I made is this: d(U_t+1)/dt = [-grad(\Gamma_e)+R_i]/exp(U_t) > \Gamma_e ,R_i : Load the value calculated by another function. > U_t : Solution calculated at the previous time. > > When implicit_runge_kutta.evolve_one_time_step(..) and > evaluate_diffusion_energy(..) is repeated, the value(U_t) I load changes. > U_t, the value that should be divided, becomes 0 and an error occurs. > What becomes zero? U_t or exp(initial_EnergyDensity)? Where is initial_EnergyDensity updated?
Best, Bruno -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/d507c2d9-9b9e-4df2-822b-d57d166c974d%40googlegroups.com.
