Hi Ingo, I'm not aware of any, but I had a quick google to find the formula (I know vaguely what divergence and curl are having watched a 3blue1brown video about it some time ago). I couldn't find any implementations in Python or Fortran.
This (https://openstax.org/books/calculus-volume-3/pages/6-5-divergence-and-curl formula 6.17) indicates the formula for curl is: If F=⟨P,Q,R⟩ is a vector field in R3, and Px, Py, Pz, Qy, Qx, Qz, Rz, Rx, and Ry all exist, then the curl of F is defined by curl F = (Ry−Qz)i + (Pz−Rx)j + (Qx−Py)k = (∂R/∂y − ∂Q/∂z)i + (∂P/∂z − ∂R/∂x)j + (∂Q/∂x − ∂P/∂y)k This suggests to me that for a 3D problem, you need a coordinates ndarray dim (3,x,y,z) and a vector field ndarray with the same dims. You can do the partial differentiation numerically by some mv-ing and then diffover (https://metacpan.org/pod/PDL::Ufunc#diffover), then the notional i,j,k above obviously indicate the final components of curl vector at each point. More to follow after I've figured out how to do the partial stuff! Best regards, Ed ________________________________ From: Ingo Schmid via pdl-devel <pdl-devel@lists.sourceforge.net> Sent: 19 November 2024 6:06 PM To: perldl <pdl-gene...@lists.sourceforge.net>; pdl-devel <pdl-devel@lists.sourceforge.net> Subject: [Pdl-devel] curl of vector Hi, is there any implementations of calculating the curl of a vector field around? Best wishes Ingo
_______________________________________________ pdl-devel mailing list pdl-devel@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/pdl-devel
