Dear Wolfgang,
Thank you for your kind reply.
I have managed to fix the problem, I was only updating the hessians in the
residuals, but not in my LHS matrix.
Regards,
Alex
El vie., 16 oct. 2020 a las 8:52, Wolfgang Bangerth ()
escribió:
>
> > I am currently implementing the entropy
I was able to diagnose and fix the problem. Thank you. I still have one
remaining error since changing all my matrices to FullMatrix. I have
diagnosed that the issue is the very last line of code below, when
attempting to update my system. I have printed the line in red.
for (unsigned int m =
I am currently implementing the entropy viscosity model (see Guermond et al.,
J. Sci. Comput. 49:35–50 (2011)), which is similar to the artificial viscosity
model introduced in step-31 for the transient advection-diffusion equation for
temperature. However, in this case, the viscosity is for
Nikki,
I can't tell from these pieces of code what dimension is wrong. But you will
find it very useful to run your program in a debugger so that you can see
which matrix access it is that triggers the error, and why. There are a number
of video lectures about running a program in a
Dear all,
I hope everyone is safe and healthy.
I am currently implementing the entropy viscosity model (see Guermond et
al., J. Sci. Comput. 49:35–50 (2011)), which is similar to the artificial
viscosity model introduced in step-31 for the transient advection-diffusion
equation for
Dear all,
I hope you are well :).
We are currently facing a weird bottle neck. There are some operation which
require us to set a particle property to zero and I think we are doing
something which is not optimal at all. Right now, we loop over the
particles, then we get the properties array