Hello, I'm using FiPy to solve a Convection Diffusion Reaction Equation (one spacial dimension, plus time). I want my source to:
da^2(H(a, t)) * P(a, t) Where da^2(H(a,t)) is the second special derivative of a function H (a known function and NOT the function I'm solving for) and P the function I'm solving for. Both are represented as CellVariable. I have tried to write: H.getGrad().getGrad() * P but the ranks don't match (H has rank-0 and H.getGrad().getGrad() has rank-3). I don't get why the rank is increased when I take the gradient of a function that is defined on a 1-d mesh. Any ideas? Thanks in advance, Matej _______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
