CFEASolver.cpp File Reference

Main subroutines for solving direct FEM elasticity problems. More...

#include "../../include/solvers/CFEASolver.hpp"
#include "../../include/variables/CFEABoundVariable.hpp"
#include "../../include/numerics/elasticity/CFEAElasticity.hpp"
#include "../../../Common/include/toolboxes/printing_toolbox.hpp"
#include "../../../Common/include/toolboxes/geometry_toolbox.hpp"
#include "../../include/solvers/CHeatSolver.hpp"
#include <algorithm>

Macros
#define	SUBTRACT_PROJECTION(READ, WRITE)

Functions
template<class T , class U , su2enable_if< is_same< T, U >::value > = 0>
CSysVector< T >	computeLinearResidual (const CSysMatrix< T > &A, const CSysVector< U > &x, const CSysVector< U > &b)

Detailed Description

Main subroutines for solving direct FEM elasticity problems.

Author

R. Sanchez

Version

8.5.0 "Harrier"

SU2 Project Website: https://su2code.github.io

The SU2 Project is maintained by the SU2 Foundation (http://su2foundation.org)

Copyright 2012-2026, SU2 Contributors (cf. AUTHORS.md)

SU2 is free software; you can redistribute it and/or modify it under the terms of the GNU Lesser General Public License as published by the Free Software Foundation; either version 2.1 of the License, or (at your option) any later version.

SU2 is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License for more details.

You should have received a copy of the GNU Lesser General Public License along with SU2. If not, see http://www.gnu.org/licenses/.

Macro Definition Documentation

SUBTRACT_PROJECTION

#define SUBTRACT_PROJECTION	(		READ,
			WRITE
	)

Value:

    SubtractProjection(nDim, normal, [&](unsigned short iDim) { return nodes->READ(iPoint, iDim); },    \
                       [&](unsigned short iDim, const su2double& x) { nodes->WRITE(iPoint, iDim, x); });

Function Documentation

computeLinearResidual()

template<class T , class U , su2enable_if< is_same< T, U >::value > = 0>

CSysVector< T > computeLinearResidual	(	const CSysMatrix< T > &	A,
		const CSysVector< U > &	x,
		const CSysVector< U > &	b
	)