/usr/include/InsightToolkit/Numerics/FEM/itkFEMSolverHyperbolic.h is in libinsighttoolkit3-dev 3.20.1-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 | /*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: itkFEMSolverHyperbolic.h
Language: C++
Date: $Date$
Version: $Revision$
Copyright (c) Insight Software Consortium. All rights reserved.
See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.
This software is distributed WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
PURPOSE. See the above copyright notices for more information.
=========================================================================*/
#ifndef __itkFEMSolverHyperbolic_h
#define __itkFEMSolverHyperbolic_h
#include "itkFEMSolver.h"
namespace itk {
namespace fem {
/**
* \class SolverHyperbolic
* \brief Solver class suitable for hyperbolic problems.
*
* M*ddu + C*du + K*u=F
*
*/
class SolverHyperbolic : public Solver
{
public:
/**
* Default constructor
*/
SolverHyperbolic();
/**
* Initialize the linear system wrapper.
*/
virtual void InitializeLinearSystemWrapper(void);
/**
* When assembling the element matrix into master matrix, we
* need to assemble the mass matrix too.
*/
virtual void AssembleElementMatrix(Element::Pointer e);
/**
* Initializes the storasge for all master matrices.
*/
virtual void InitializeMatrixForAssembly(unsigned int N);
/**
* Combines the M, C and K matrices into one big system of linear
* equations.
*/
virtual void FinalizeMatrixAfterAssembly( void );
/**
* Solves the system for the next time step.
*/
virtual void Solve( void );
virtual Float GetTimeStep( void ) const { return m_deltaT; }
virtual void SetTimeStep(Float dt) { this->m_deltaT=dt; }
/**
* Constants that specify, where matrices are strored.
*/
enum { matrix_K=1, matrix_M=2, matrix_C=3, matrix_tmp=4 };
/**
* Constants that specify, where vectors are strored.
*/
enum { solution_d=0, solution_v=1, solution_a=2};
enum { vector_dhat=2, vector_vhat=3, vector_ahat=4, vector_tmp=5 };
Float m_gamma;
Float m_beta;
Float m_deltaT;
};
}} // end namespace itk::fem
#endif // #ifndef __itkFEMSolverHyperbolic_h
|