Code: utility

Under construction

• Operator computation: A, B, C, D, E (mass matrix)
• Frequency response computation: H(jw) = C*inv(sE-A)^{-1}*B(+D) with real linear algebra
• Export utils (spy, save)
• Debug utils (export subdomains)
• Mpi utils (= shortcuts)
• Use it as backend in an optimization tool (as in Jussiau, W., Demourant, F., Leclercq, C., & Apkarian, P. (2025). Control of a Class of High-Dimensional Nonlinear Oscillators: Application to Flow Stabilization. IEEE Transactions on Control Systems Technology.).