backend.cc

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) NumeriCor GmbH
// SPDX-License-Identifier: Apache-2.0
 
#include "backend.h"
#include "basics.h"
 
#include <cstdlib>
#include <cstring>
#include <string>
 
namespace opencarp {
 
namespace {
 
bool petsc_running = false;
bool mpi_started_here = false;
 
void stop_if_flavor_is_missing()
{
  const std::string flavor = param_globals::flavor != nullptr ? param_globals::flavor : "";
#ifndef WITH_GINKGO
  if (flavor == "ginkgo") {
    fprintf(stderr, "\n*** The flavor \"ginkgo\" is not available, openCARP was built without the Ginkgo backend.\n\n");
    exit(EXIT_FAILURE);
  }
#endif
#ifndef WITH_PETSC
  if (!flavor.empty() && flavor != "ginkgo") {
    fprintf(stderr, "\n*** The flavor \"%s\" is not available, openCARP was built with only the Ginkgo backend.\n\n",
            flavor.c_str());
    exit(EXIT_FAILURE);
  }
#endif
}
 
}  // namespace
 
bool ginkgo_selected()
{
#ifdef WITH_PETSC
  return param_globals::flavor != nullptr && std::string{param_globals::flavor} == "ginkgo";
#else
  return true;
#endif
}
 
bool petsc_active()
{
  return petsc_running;
}
 
void initialize_backend(int *argc, char **argv, char *rc, char *help_msg)
{
  stop_if_flavor_is_missing();
 
  if (!ginkgo_selected()) {
#ifdef WITH_PETSC
    initialize_PETSc(argc, argv, rc, help_msg);
#endif
    return;
  }
 
  int delimiter = 1;
  while (delimiter < *argc && strcmp(argv[delimiter], DELIMIT_OPTS) != 0) delimiter++;
  const bool has_petsc_options = delimiter + 1 < *argc;
  *argc = delimiter;
 
  int initialized = 0;
  MPI_Initialized(&initialized);
  if (!initialized) {
    // the level PETSc asks for: the OpenMP kernels are threaded, MPI is only used outside of them
    int provided = 0;
    MPI_Init_thread(nullptr, nullptr, MPI_THREAD_FUNNELED, &provided);
    mpi_started_here = true;
  }
 
  if (has_petsc_options)
    log_msg(NULL, 3, 0, "The options after \"%s\" are PETSc options, which the Ginkgo backend ignores.", DELIMIT_OPTS);
}
 
void finalize_backend()
{
#ifdef WITH_PETSC
  if (petsc_running) {
    finalize_PETSc();
    return;
  }
#endif
  if (mpi_started_here) {
    mpi_started_here = false;
    MPI_Finalize();
  }
}
 
void finalize_backend_and_exit(int code)
{
  finalize_backend();
  exit(code);
}
 
#ifdef WITH_PETSC
void set_petsc_active(bool active)
{
  petsc_running = active;
}
#endif
 
}  // namespace opencarp