On Fri, Jun 6, 2014 at 11:36 AM, Zebo LI <[email protected]> wrote:
> Sorry,
>
> Here is the code. In this case, if we just set velo=epsilon.grad, then
> the code works well. But if we use velo = DOT(D,epsilon.grad), it doesn't.
> (Actually, D is the unity Matrix ((1,0),(0,1)).)
Change
velo = DOT(D,epsilon.grad)
to
velo = epsilon.grad.dot(D)
and it should work for you. Not sure why "DOT" isn't returning a
variable object.
--
Daniel Wheeler
_______________________________________________
fipy mailing list
[email protected]
http://www.ctcms.nist.gov/fipy
[ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
