On Sep 30, 2013, at 4:38 AM, Damian Kösters <koesters.dam...@gmail.com> wrote:
> I try now to work with a Constraint object. This doesn't give me an error > message but I also do not see the rectangle in the plot. > >>>x, y = mesh.cellCenters[0], mesh.cellCenters[1] > >>>X, Y = mesh.faceCenters[0], mesh.faceCenters[1] > >>> > >>>rect2 = (x>9)&(x<22)&(y>50)&(y<100) > >>> > >>>crect2 = Constraint(10,rect2 ) > >>> > >>>rect2.constrain(crect2) Plot of what? This code constrains the boolean mask rect2 to have value of 10 where that mask is true (I think). I think you want phi.constrain(crect2). _______________________________________________ fipy mailing list fipy@nist.gov http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]