https://gitlab.opengeosys.org/ogs/ogs.git
Tip revision: 2074dab4c68311146369efc8639b882c03d864ce authored by Lars Bilke on 15 December 2020, 07:52:41 UTC
[CMake] Fixed definition USE_PETSC.
[CMake] Fixed definition USE_PETSC.
Tip revision: 2074dab
CoupledSolutionsForStaggeredScheme.h
/**
* \copyright
* Copyright (c) 2012-2020, OpenGeoSys Community (http://www.opengeosys.org)
* Distributed under a Modified BSD License.
* See accompanying file LICENSE.txt or
* http://www.opengeosys.org/project/license
*
* \file
*
* Created on November 7, 2016, 12:14 PM
*/
#pragma once
#include <utility>
#include <vector>
#include "MathLib/LinAlg/GlobalMatrixVectorTypes.h"
namespace ProcessLib
{
/**
* A struct to keep the references of the current solutions of the equations of
* the coupled processes.
*
* During staggered coupling iteration, an instance of this struct is created
* and passed through interfaces to global and local assemblers for each
* process.
*/
struct CoupledSolutionsForStaggeredScheme
{
explicit CoupledSolutionsForStaggeredScheme(
std::vector<GlobalVector*> const& coupled_xs_);
/// References to the current solutions of the coupled processes.
std::vector<GlobalVector*> const& coupled_xs;
/// Pointers to the vector of the solutions of the previous time step.
std::vector<GlobalVector*> coupled_xs_t0;
};
/**
* A struct to keep the references to the local element solutions of the
* current and previous time step solutions of the equations of the coupled
* processes.
*
* During the global assembly loop, an instance of this struct is created for
* each element and it is then passed to local assemblers.
*/
struct LocalCoupledSolutions
{
explicit LocalCoupledSolutions(std::vector<double>&& local_coupled_xs0_)
: local_coupled_xs0(std::move(local_coupled_xs0_))
{
}
/// Local solutions of the previous time step.
std::vector<double> const local_coupled_xs0;
};
/**
* Fetch the nodal solutions of all coupled processes from the given vector of
* global solutions for each process into a flat vector.
*/
std::vector<double> getCoupledLocalSolutions(
std::vector<GlobalVector*> const& global_solutions,
std::vector<std::vector<GlobalIndexType>> const& indices);
} // namespace ProcessLib