Closed Pbellive closed 7 years ago
I wonder if we should drop redundant output argument(s) in functions with exclamation mark which modify the input argument(s). The function concerned here is DiagTimesMTimesDiag!
. I find it confusing to see, e.g., CURL = DiagTimesMTimesDiag!(FSZi, T, ESZ)
where CURL
and T
will point to the same memory location. What's the official Julia rule on this?
I think the Julia best practice is to avoid the redundant outputs. We should do this.
Go back to using old method of constructing edge mass matrices for isotropic models. Provides massive speedups (20X and higher ) on mass matrix and derivative of mass matrix construction.