定义球面坐标中的单位向量以与Eigen3一起使用

Question

我在这里发帖是因为我有点迷茫

我正在尝试实现一个solution to Maxwells equations (p47 2-2)

​

，

它是以C++的球面坐标给出的，因此它可以用于更大的建模项目中。我使用Eigen3作为线性代数的基础，据我所知，它并不显式地支持球面坐标(我对其他选择持开放态度)。

为了实现解决方案，我需要(或者至少我认为我需要)将球面单位向量定义为球面坐标，就像笛卡尔坐标中的they're not constants一样，我不知道如何做到这一点。

我很犹豫是否将解决方案转换为笛卡尔坐标，因为我认为我不理解这样做的含义(它是否有效？)

任何和所有的意见和建议都是值得感谢的

Answer 1

解决方案是将Spherical Unit Vector Identities实现为3个函数(每个单位向量一个)，将r、Theta和Phi作为参数并返回一个向量。