svMultiPhysics
Loading...
Searching...
No Matches
cmm.h
1/* Copyright (c) Stanford University, The Regents of the University of California, and others.
2 *
3 * All Rights Reserved.
4 *
5 * See Copyright-SimVascular.txt for additional details.
6 *
7 * Permission is hereby granted, free of charge, to any person obtaining
8 * a copy of this software and associated documentation files (the
9 * "Software"), to deal in the Software without restriction, including
10 * without limitation the rights to use, copy, modify, merge, publish,
11 * distribute, sublicense, and/or sell copies of the Software, and to
12 * permit persons to whom the Software is furnished to do so, subject
13 * to the following conditions:
14 *
15 * The above copyright notice and this permission notice shall be included
16 * in all copies or substantial portions of the Software.
17 *
18 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS
19 * IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
20 * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A
21 * PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER
22 * OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
23 * EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
24 * PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
25 * PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
26 * LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
27 * NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
28 * SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
29 */
30
31#ifndef CMM_H
32#define CMM_H
33
34#include "ComMod.h"
35#include "Simulation.h"
36
37/// @brief These subroutines implement the Coupled Momentum Method (CMM).
38namespace cmm {
39
40void cmm_3d(ComMod& com_mod, const int eNoN, const double w, const Vector<double>& N, const Array<double>& Nx,
41 const Array<double>& al, const Array<double>& yl, const Array<double>& bfl, const Array<double>& Kxi,
42 Array<double>& lR, Array3<double>& lK);
43
44void cmm_b(ComMod& com_mod, const faceType& lFa, const int e, const Array<double>& al, const Array<double>& dl,
45 const Array<double>& xl, const Array<double>& bfl, const Vector<double>& pS0l, const Vector<double>& vwp,
46 const Vector<int>& ptr);
47
48void bcmmi(ComMod& com_mod, const int eNoN, const int idof, const double w, const Vector<double>& N, const Array<double>& Nxi,
49 const Array<double>& xl, const Array<double>& tfl, Array<double>& lR);
50
51void cmmi(ComMod& com_mod, const mshType& lM, const Array<double>& al, const Array<double>& dl, const Array<double>& xl,
52 const Array<double>& bfl, const Array<double>& pS0l, const Vector<double>& vwp, const Vector<int>& ptr);
53
54void cmm_mass(ComMod& com_mod, const double w, const Vector<double>& N, const Array<double>& al,
55 const Array<double>& bfl, const Vector<double>& vwp, Array<double>& lR, Array3<double>& lK);
56
57void cmm_stiffness(ComMod& com_mod, const Array<double>& Nxi, const Array<double>& xl, const Array<double>& dl,
58 const Vector<double>& pS0l, const Vector<double>& vwp, Vector<double>& pSl, Array<double>& lR, Array3<double>& lK);
59
60void construct_cmm(ComMod& com_mod, const mshType& lM, const Array<double>& Ag, const Array<double>& Yg, const Array<double>& Dg);
61
62};
63
64#endif
65
The Array3 template class implements a simple interface to 3D arrays.
Definition Array3.h:52
The ComMod class duplicates the data structures in the Fortran COMMOD module defined in MOD....
Definition ComMod.h:1354
The Vector template class is used for storing int and double data.
Definition Vector.h:50
The face type containing mesh at boundary.
Definition ComMod.h:521
This is the container for a mesh or NURBS patch, those specific to NURBS are noted.
Definition ComMod.h:832
These subroutines implement the Coupled Momentum Method (CMM).
Definition cmm.cpp:43
void cmmi(ComMod &com_mod, const mshType &lM, const Array< double > &al, const Array< double > &dl, const Array< double > &xl, const Array< double > &bfl, const Vector< double > &pS0l, const Vector< double > &vwp, const Vector< int > &ptr)
CMM initialization (interior).
Definition cmm.cpp:378
void construct_cmm(ComMod &com_mod, const mshType &lM, const Array< double > &Ag, const Array< double > &Yg, const Array< double > &Dg)
Reproduces Fortran 'CONSTRUCT_CMM'.
Definition cmm.cpp:785