thank you for the response - I will try your suggestion on the anisotropic diffusion coeff.
is it possible the conversion to python 3 is the issue with the equivalency in the case of a constant anisotropic diffusion coeff. ? My feeling is that [[1., 5.]] shoud be equivalent to [[[1., 0.], [0., 5.]]]. I'm not sure why that isn't working correctly. On 4/18/2016 8:54 AM, Daniel Wheeler wrote: > On Fri, Apr 15, 2016 at 9:29 AM, Phil <[email protected]> wrote: >> But when written like >> >> TransientTerm() == DiffusionTerm([[1.,5.]]), which I think the >> manual said was equivalent ( under freq asked questions) >> I get an assertion error, which I can post but before posting should >> these two diffusion terms be equivalent? > My feeling is that [[1., 5.]] shoud be equivalent to [[[1., 0.], [0., > 5.]]]. I'm not sure why that isn't working correctly. > >> Sticking with the complete matrix form of the diffusion Coefficient >> for Dx and Dz ( both similar functions of the cellVariable C) >> >> when I use this form >> >> TransientTerm() == DiffusionTerm([[[Dx,0],[0,Dz]]]) >> >> I get the following traceback error > That's not surprising. Use a single FaceVariable or CellVariable of > rank 2 for the diffusion coefficient. It's safer to do that for > anisotropy. > >> I see questions posted on this list, Gist and stack overflow - is there a >> preference for posting ? > I prefer Stackoverflow, but depends on the nature of the question. > Stackoverflow isn't good for very open ended questions or > conversations. > _______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
