fix(main): setup

added moment tensor and COM calculation to setup
2026-02-18 07:37:15 -05:00
parent 7d047ead16
commit 423914bc59
1 changed files with 3 additions and 0 deletions
--- 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");
  }