I feel this should be an elementary question, but I can't seem to figure
out how to answer it. I am solving a simple linear elleptic-ish equation
with

eq =
(DiffusionTerm(var=Psi,coeff=DiffCoeff)+ExponentialConvectionTerm(var=Psi,coeff=convCoeff))
eq.solve(var=Psi)


This works fine; the solution matches what I would expect.  I have two
questions:

First, how can I obtain the value of the individual terms of the equation,
evaluated with the solution in Psi?

Second, is there any way to define my own new Psi (which is not an exact
solution to the equation), and easily evaluate the DiffusionTerm and
ExponentialConvectionTerm for that Psi?

I am trying to illustrate how various parts of the solution evolve over the
solution space, and how various approximations to the full solution differ
from the full solution.

Thanks to all who glance at this
James Pringle
University of New Hampshire
_______________________________________________
fipy mailing list
[email protected]
http://www.ctcms.nist.gov/fipy
  [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]

Reply via email to