ComputeResiduum.h
/**
* \file
* \copyright
* Copyright (c) 2012-2021, OpenGeoSys Community (http://www.opengeosys.org)
* Distributed under a Modified BSD License.
* See accompanying file LICENSE.txt or
* http://www.opengeosys.org/project/license
*/
#pragma once
#include "MathLib/LinAlg/GlobalMatrixVectorTypes.h"
namespace ProcessLib
{
/// Computes the residuum r = -M*x_dot - K*x + b. Negation for consistency with
/// the Newton scheme.
GlobalVector computeResiduum(GlobalVector const& x, GlobalVector const& xdot,
GlobalMatrix const& M, GlobalMatrix const& K,
GlobalVector const& b);
} // namespace ProcessLib