search

DedalusProject / dedalus

A flexible framework for solving PDEs with modern spectral methods.
http://dedalus-project.org/
GNU General Public License v3.0
489 stars 115 forks source link

CartesianTrace with rank > 2 #292

Closed lecoanet closed 3 months ago

lecoanet commented 3 months ago

The matrix form of the CartesianTrace operator assumed the tensor had rank 2. I've modified it to work for rank >=2. I've also added tests for the rank=3 case.