On Wed, Aug 3, 2016 at 9:27 AM, Nils Becker <nils.bec...@bioquant.uni-heidelberg.de> wrote: >> Could you possible post your code so I can play around with it? I >> can't think offhand what could be the issue or exactly what the error >> is that you're calculating, it seems like the normalized residual or >> something. > > hi daniel, thanks for your reply. the way i was calculating the residual > is to export the solution on a grid, import that into mathematica, > interpolate in space and time and have mathematica calculate > ||lhs-rhs|| / ||lhs|| where ||...|| is the L2 norm. what i saw was that > the relative error was low in the beginning but increased over time, > especially in places where the solution starts to grow exponentially.
I'm not sure if using ||lhs|| is a good way to normalize the residual. Surely that could be very small. Is ||lhs - rhs|| getting smaller without normalization? -- Daniel Wheeler _______________________________________________ fipy mailing list fipy@nist.gov http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]