Closed DavidSagan closed 5 months ago
Currently my checksymp returns 1-norm of transpose(M)JM - J. I'm thinking about instead making checksymp just return the matrix transpose(M)JM - J , and then the user could calculate the norm if it is to be a quantified as just a single scalar number. How does that sound?
Yes that sounds more modular. Do it.
added in latest commit 6a87b1f
This is helpful since looking at the significantly nonzero terms in this matrix gives clues as to where the non-symplecticity is happening.