From 423914bc59668f925f3bfc0be3c232673e8c38f5 Mon Sep 17 00:00:00 2001 From: Emily Boudreaux Date: Wed, 18 Feb 2026 07:37:15 -0500 Subject: [PATCH] fix(main): setup added moment tensor and COM calculation to setup --- main.cpp | 3 +++ 1 file changed, 3 insertions(+) diff --git a/main.cpp b/main.cpp index ae8ed68..2fa5b9e 100644 --- a/main.cpp +++ b/main.cpp @@ -339,6 +339,9 @@ FEM setup_fem(const std::string& filename, bool verbose) { project_scalar_function(*fem_setup.rho_gf, initial_density); fem_setup.rho_gf = conserve_mass(fem_setup, *fem_setup.rho_gf, MASS); + fem_setup.com = get_com(fem_setup, *fem_setup.rho_gf); + fem_setup.Q = compute_quadrupole_moment_tensor(fem_setup, *fem_setup.rho_gf, fem_setup.com); + if (verbose) { std::println("Setup {}", fem_setup.okay() ? "OK" : "FAIL"); }