Test for Disentangled Spatiotemporal SSM

Testing/OptimizeTests/OptimizeTests.cpp

@@ -696,10 +696,38 @@ TEST(OptimizeTests, vtk_output) {
   for (int i = 0; i < values.size(); i++) {
     std::cerr << "Eigenvalue " << i << " : " << values[i] << "\n";
   }
-
   // check the first mode of variation.
   // If Procrustes scaling is working, this should be small.
   // Otherwise it is quite large (>4000).
   double value = values[values.size() - 1];
   ASSERT_LT(value, 100);
 }
+
+//---------------------------------------------------------------------------
+TEST(OptimizeTests, disentangled_spatiotemporal_test) {
+  prep_temp("/optimize/disentangled_spatiotemporal", "disentangled_spatiotemporal_test");
+
+  std::remove("test1_particles/ID01_TIME01_world.particles");
+
+  Optimize app;
+  ProjectHandle project = std::make_shared<Project>();
+  ASSERT_TRUE(project->load("test1.xlsx"));
+  OptimizeParameters params(project);
+  ASSERT_TRUE(params.set_up_optimize(&app));
+  bool runStatus = app.Run();
+
+  ASSERT_TRUE(runStatus);
+
+  // compute and test stats
+  ParticleShapeStatistics stats;
+  stats.ReadPointFiles("test1_analyze.xml");
+  stats.ComputeModes();
+  stats.PrincipalComponentProjections();
+
+  auto values = stats.Eigenvalues();
+  for (int i = 0; i < values.size(); i++) {
+    std::cerr << "Eigenvalue " << i << " : " << values[i] << "\n";
+  }
+  ASSERT_GT(values[values.size() - 1], 200000);
+  ASSERT_LT(values[values.size() - 2], 120000);
+}

Testing/data/optimize/disentangled_spatiotemporal/test1_analyze.xml

@@ -0,0 +1,66 @@
+<point_files>
+  test1_particles/ID01_TIME01_world.particles
+  test1_particles/ID01_TIME02_world.particles
+  test1_particles/ID01_TIME03_world.particles
+  test1_particles/ID01_TIME04_world.particles
+  test1_particles/ID01_TIME05_world.particles
+  test1_particles/ID01_TIME06_world.particles
+  test1_particles/ID01_TIME07_world.particles
+  test1_particles/ID01_TIME08_world.particles
+  test1_particles/ID02_TIME01_world.particles
+  test1_particles/ID02_TIME02_world.particles
+  test1_particles/ID02_TIME03_world.particles
+  test1_particles/ID02_TIME04_world.particles
+  test1_particles/ID02_TIME05_world.particles
+  test1_particles/ID02_TIME06_world.particles
+  test1_particles/ID02_TIME07_world.particles
+  test1_particles/ID02_TIME08_world.particles
+  test1_particles/ID03_TIME01_world.particles
+  test1_particles/ID03_TIME02_world.particles
+  test1_particles/ID03_TIME03_world.particles
+  test1_particles/ID03_TIME04_world.particles
+  test1_particles/ID03_TIME05_world.particles
+  test1_particles/ID03_TIME06_world.particles
+  test1_particles/ID03_TIME07_world.particles
+  test1_particles/ID03_TIME08_world.particles
+  test1_particles/ID04_TIME01_world.particles
+  test1_particles/ID04_TIME02_world.particles
+  test1_particles/ID04_TIME03_world.particles
+  test1_particles/ID04_TIME04_world.particles
+  test1_particles/ID04_TIME05_world.particles
+  test1_particles/ID04_TIME06_world.particles
+  test1_particles/ID04_TIME07_world.particles
+  test1_particles/ID04_TIME08_world.particles
+  test1_particles/ID05_TIME01_world.particles
+  test1_particles/ID05_TIME02_world.particles
+  test1_particles/ID05_TIME03_world.particles
+  test1_particles/ID05_TIME04_world.particles
+  test1_particles/ID05_TIME05_world.particles
+  test1_particles/ID05_TIME06_world.particles
+  test1_particles/ID05_TIME07_world.particles
+  test1_particles/ID05_TIME08_world.particles
+  test1_particles/ID06_TIME01_world.particles
+  test1_particles/ID06_TIME02_world.particles
+  test1_particles/ID06_TIME03_world.particles
+  test1_particles/ID06_TIME04_world.particles
+  test1_particles/ID06_TIME05_world.particles
+  test1_particles/ID06_TIME06_world.particles
+  test1_particles/ID06_TIME07_world.particles
+  test1_particles/ID06_TIME08_world.particles
+  test1_particles/ID07_TIME01_world.particles
+  test1_particles/ID07_TIME02_world.particles
+  test1_particles/ID07_TIME03_world.particles
+  test1_particles/ID07_TIME04_world.particles
+  test1_particles/ID07_TIME05_world.particles
+  test1_particles/ID07_TIME06_world.particles
+  test1_particles/ID07_TIME07_world.particles
+  test1_particles/ID07_TIME08_world.particles
+  test1_particles/ID08_TIME01_world.particles
+  test1_particles/ID08_TIME02_world.particles
+  test1_particles/ID08_TIME03_world.particles
+  test1_particles/ID08_TIME04_world.particles
+  test1_particles/ID08_TIME05_world.particles
+  test1_particles/ID08_TIME06_world.particles
+  test1_particles/ID08_TIME07_world.particles
+  test1_particles/ID08_TIME08_world.particles
+</point_files>