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)).)
Zebo
2014-06-06 9:28 GMT-05:00 Daniel Wheeler :
> On Thu, Jun 5, 2014 at 12:23 PM, Zebo
On Thu, Jun 5, 2014 at 12:23 PM, Zebo LI wrote:
> Hi,
>
> I currently working on a diffusion problem using Fipy. And I need to
> calculate the convection velocity via the product of a matrix and a vector.
> I used numeric.array and numeric.NUMERIX.dot, but it does not work well, as
> shown in
