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");
   }