Function opensn::MaxEigenvalueOnStiffnessRange

Defined in File boundary_closure.h

Function Documentation

double opensn::MaxEigenvalueOnStiffnessRange(const DenseMatrix<double> &N, const DenseMatrix<double> &K)

Largest generalized eigenvalue \(\lambda\) of \(N v = \lambda K v\) over the range of K, with K the cell stiffness matrix \(\int_K \nabla b_i \cdot \nabla b_j\) (symmetric positive semi-definite, null space the constants) and N the face matrix \(\int_f \partial_n b_i \partial_n b_j\), which also vanishes on constants.

Times \(h_f\) this is the trace constant \(C_{grad}\) of the basis on the face.