Hello, On Wed, Jan 28, 2009 at 4:50 PM, Jonathan Guyer <[email protected]> wrote:
> In FiPy, you would write them both the same way, simply changing from a 1D > to a 2D Mesh. Yes, I know that the difference is only the Mesh. > Please see the Cahn-Hilliard examples in the manual. I saw the example you mentioned but I am not sure if I did correct in my case. Just I wanted to compare with your version, that's why I asked for little help. Here is what I have right now: [..skip...] faceVar = var.getArithmeticFaceValue() fDerivative = ..... from fipy.terms.implicitDiffusionTerm import ImplicitDiffusionTerm from fipy.terms.transientTerm import TransientTerm diffTerm4 = - ImplicitDiffusionTerm(coeff = (h*h*h, -fDerivative)) eqch = TransientTerm() + diffTerm4 from fipy.solvers.linearPCGSolver import LinearPCGSolver solver = LinearPCGSolver(tolerance = 1e-15,steps = 1000) [..skip...] Please let me know what you think, is it correct or I forgot something to add? Thanks, Ionut
