>> 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?
if you don't see any obvious errors in building the solution in fipy, i could continue poking around if i knew how to evaluate the various terms of the solution and their errors at different times. i see methods eq_diffprolif.justErrorvector, but i don't quite understand what exactly this returns. is it lhs-rhs? what is justResidualVector? can i get those for the individual terms on the rhs as well? cheers, n.
0xF0A5C638.asc
Description: application/pgp-keys
_______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
