The following runs for me with a few changes. On Fri, May 6, 2016 at 1:25 PM, Mohammad Kazemi <kazem...@gmail.com> wrote: > > Pk.constrain([Pd], mesh.facesLeft) > Pk.constrain([Pu], mesh.facesRight)
Pk.constrain(Pd.faceValue, mesh.facesLeft) Pk.constrain(Pu.faceValue, mesh.facesRight) > Pf.constrain([Pd], mesh.facesLeft) > Pf.constrain([Pu], mesh.facesRight) Pf.constrain(Pd.faceValue, mesh.facesLeft) Pf.constrain(Pu.faceValue, mesh.facesRight) This makes sense since the Pf and Pk are scalars and [Pd] and [Pu] are effectively vectors. Also the face value is required since the constraint is on the faces. I hope this helps. -- Daniel Wheeler _______________________________________________ fipy mailing list fipy@nist.gov http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
