Commit f86fc6c8 authored by Giovanni La Mura's avatar Giovanni La Mura
Browse files

Implement runtime inversion selector for CLUSTER

parent 70f1e6a8

src/cluster/cluster.cpp

+13 −0
Original line number Diff line number Diff line
@@ -234,6 +234,19 @@ void cluster(const string& config_file, const string& data_file, const string& o 
      sprintf(virtual_line, "%.5lg m.\n", sconf->get_particle_radius(gconf));

      message = "INFO: particle radius is " + (string)virtual_line;

      logger->log(message);

      if (gconf->invert_mode == RuntimeSettings::INV_MODE_LU) {

	message = "INFO: using LU factorization for inversion.\n";

	logger->log(message);

      } else if(gconf->invert_mode == RuntimeSettings::INV_MODE_GESV) {

	message = "INFO: using GESV system solver for inversion.\n";

	logger->log(message);

      } else if(gconf->invert_mode == RuntimeSettings::INV_MODE_RBT) {

	message = "INFO: using RBT for inversion.\n";

	logger->log(message);

      } else if(gconf->invert_mode == RuntimeSettings::INV_MODE_SVD) {

	message = "INFO: using SVD for inversion.\n";

	logger->log(message);

      }

      // Overlapping spheres test

      double tolerance = gconf->tolerance;

      if (tolerance < 0.0) {