Dear fipy users, I'm trying to solve a problem in an anisotropic system. Referring to 12.4 examples.phase.anisotropy, I use two different methods to define the orientation as the following: theta=arctan2(phi.getFaceGrad()[1], phi.getFaceGrad()[0]) theta=arctan(phi.getFaceGrad()[1] / phi.getFaceGrad()[0]) However, I got a slight different result from above. Can someone tell me what's the difference between these two methods? Which one is more correct to define the theta? And, when I put the theta into the equation of psi=1+0.5cos(4*theta) to calcalute the my anisotropy system , it is expected there is a fourfold symmetry. Unforturnately, only threefold symmetry is appeared. If I try 1+0.5cos(5*theta), only fourfold symmetry is shown. I'm not sure what's problem in my calculation, in the definition of theta or the governing equation. Hope someone can help me to clarify these problem. Thanks a lot in advance.
Regards, Kundar _______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
