OptUtils.C

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#include "OptUtils.h"
#include "libmesh/id_types.h"
#include "libmesh/petsc_vector.h"
#include "libmesh/petsc_matrix.h"

// static functions to copy Petsc vectors to and from vectors for a Reporter

namespace OptUtils
{

void
copyReporterIntoPetscVector(const std::vector<std::vector<Real> *> reporterVectors,
                            libMesh::PetscVector<Number> & x)
{
  dof_id_type n = 0;
  for (const auto & data : reporterVectors)
    for (const auto & val : *data)
      x.set(n++, val);

  x.close();
}

void
copyPetscVectorIntoReporter(const libMesh::PetscVector<Number> & x,
                            std::vector<std::vector<Real> *> reporterVectors)
{
  dof_id_type j = 0;
  for (auto & data : reporterVectors)
    for (auto & val : *data)
      val = x(j++);
}
void
copyReporterIntoPetscMatrix(const std::vector<std::vector<Real> *> reporterVectors,
                            libMesh::PetscMatrix<Number> & x)
{
  for (const auto i : index_range(reporterVectors))
    for (const auto j : index_range(*reporterVectors[i]))
      x.set(i, j, (*reporterVectors[i])[j]);

  x.close();
}

void
copyPetscMatrixIntoReporter(const libMesh::PetscMatrix<Number> & x,
                            std::vector<std::vector<Real> *> reporterVectors)
{
  for (const auto i : index_range(reporterVectors))
    for (const auto j : index_range(*reporterVectors[i]))
      (*reporterVectors[i])[j] = x(i, j);
}
}