My second week I spent on introducing Lame coefficients into
Unfortunately, our work is restricted by SymPy structure so, we don’t have too much freedom
in our implementation. Hopefully, with my mentor, Francesco, we found some solution how to
achieve our goals without destroying
vector module. This week shows that I have lack in my
knowledge about object-oriented programming and SymPy.
Having access to Lame coefficient I was able to modified
Del operator (in mathematics nabla operators) to handle spherical
and cylindrical coordinate system. It was rather easy task. I had to only modified i.e. generalized
cross methods to any kind of orthogonal curvilinear system simply by multiplying by Lame factors. At that moment I
stick only to this coordinate system. We can easily expand it to another in the future.
Below is example how we can calculate
divergence in spherical coordinate system:
To sum up, This week I end up with two PR, which I think will be merged soon. We also started discussion about next step. I’d say that the hardest difficulties are now in front of us.