tao_null_space.cpp File Reference

Test null-space TAO optimization of a scaled epigraph constraint. More...

#include <MoFEM.hpp>
#include <limits>

Go to the source code of this file.

Functions
int	main (int argc, char *argv[])

Variables
static char	help []

Detailed Description

Test null-space TAO optimization of a scaled epigraph constraint.

Definition in file tao_null_space.cpp.

Function Documentation

main()

int main	(	int	argc,
		char *	argv[]
	)

Definition at line 775 of file tao_null_space.cpp.

                                 {
  MoFEM::Core::Initialize(&argc, &argv, nullptr, help);
  try {
    CHKERR TaoRegisterNullSpace();
    PetscBool activation = PETSC_FALSE, domain = PETSC_FALSE,
              domain_budget = PETSC_FALSE, flow = PETSC_FALSE,
              flow_budget = PETSC_FALSE;
    CHKERR PetscOptionsGetBool(nullptr, nullptr, "-test_nullspace_activation",
                               &activation, nullptr);
    CHKERR PetscOptionsGetBool(nullptr, nullptr, "-test_nullspace_domain",
                               &domain, nullptr);
    CHKERR PetscOptionsGetBool(nullptr, nullptr,
                               "-test_nullspace_domain_budget", &domain_budget,
                               nullptr);
    CHKERR PetscOptionsGetBool(nullptr, nullptr, "-test_nullspace_flow", &flow,
                               nullptr);
    CHKERR PetscOptionsGetBool(nullptr, nullptr, "-test_nullspace_flow_budget",
                               &flow_budget, nullptr);
    if (activation) {
      CHKERR checkCandidateActivation();
    } else if (domain_budget) {
      CHKERR checkObjectiveDomainLineSearch(true);
    } else if (domain) {
      CHKERR checkObjectiveDomainLineSearch(false);
    } else if (flow || flow_budget) {
      CHKERR checkAlgorithmOneStep(flow_budget);
    } else {
      const std::array<std::array<PetscScalar, 3>, 3> starts = {
          {{1, 1, 3}, {3, 4, 2}, {1, 1, 6}}};
      for (const auto scale : {1.0, 0.1, 10.0})
        for (size_t i = 0; i != starts.size(); ++i)
          CHKERR solveProblem(starts[i], scale, i != 1);
    }
  }
  CATCH_ERRORS;
  CHKERR MoFEM::Core::Finalize();
  return 0;
}

Variable Documentation

help

char help[]

static

Initial value:

=
    "Test the null-space TAO solver on an epigraph problem.\n\n"

Definition at line 10 of file tao_null_space.cpp.