Memosa-FVM  0.2
Public Types | Public Member Functions | Private Attributes | List of all members
CometMatrix< T > Class Template Reference

#include <CometMatrix.h>

Inheritance diagram for CometMatrix< T >:
Inheritance graph
[legend]
Collaboration diagram for CometMatrix< T >:
Collaboration graph
[legend]

Public Types

typedef Array< T > TArray
 
typedef NumTypeTraits< T >
::T_Scalar 		T_Scalar
 
- Public Types inherited from MatrixJML< T >
typedef Array< T > TArray
 

Public Member Functions

 CometMatrix (const int order)
 
T & getElement (const int i, const int j)
 
void printElement (const int &i, const int &j)
 
void Solve (TArray &bVec)
 
void zero ()
 
T getTraceAbs ()
 
- Public Member Functions inherited from MatrixJML< T >
 MatrixJML ()
 
virtual ~MatrixJML ()
 
 MatrixJML ()
 
virtual ~MatrixJML ()
 
virtual void multiply (const TArray &x, TArray &b)=0
 

Private Attributes

int _elements
 
TArray _values
 
int _order
 
SquareMatrixESBGK< T > _AMat
 
TArray _vel
 

Detailed Description

template<class T>
class CometMatrix< T >

Definition at line 16 of file CometMatrix.h.

Member Typedef Documentation

template<class T >
typedef NumTypeTraits<T>::T_Scalar CometMatrix< T >::T_Scalar

Definition at line 20 of file CometMatrix.h.

template<class T >
typedef Array<T> CometMatrix< T >::TArray

Definition at line 19 of file CometMatrix.h.

Constructor & Destructor Documentation

template<class T >
CometMatrix< T >::CometMatrix ( const int  order)
inline

Definition at line 22 of file CometMatrix.h.

References CometMatrix< T >::_AMat, CometMatrix< T >::_values, and CometMatrix< T >::_vel.

22  :
23  _elements(7*order-18),
24  _values(_elements),
25  _order(order),
26  _AMat(3),
27  _vel(3)
28  {
29  _values=0.;
30  _AMat.zero();
31  _vel=0;
32  }
CometMatrix::_AMat
SquareMatrixESBGK< T > _AMat
Definition: CometMatrix.h:198
CometMatrix::_values
TArray _values
Definition: CometMatrix.h:196
CometMatrix::_vel
TArray _vel
Definition: CometMatrix.h:199

Member Function Documentation

template<class T >
T& CometMatrix< T >::getElement ( const int  i,
const int  j 
)
inlinevirtual

Implements MatrixJML< T >.

Definition at line 34 of file CometMatrix.h.

References CometMatrix< T >::_order, and CometMatrix< T >::_values.

Referenced by CometMatrix< T >::Solve().

35  {
36  int index;
37  if(i==j)
38  {
39  index=(i-1);
40  return _values[index];
41  }
42  else if(j==_order)
43  {
44  index=_order+3*i-1;
45  return _values[index];
46  }
47  else if(j==_order-1)
48  {
49  index=_order+3*i-2;
50  return _values[index];
51  }
52  else if(j==_order-2)
53  {
54  index=_order+3*i-3;
55  return _values[index];
56  }
57  else if(i==_order)
58  {
59  index=6*_order+j-16;
60  return _values[index];
61  }
62  else if(i==_order-1)
63  {
64  index=5*_order+j-13;
65  return _values[index];
66  }
67  else if(i==_order-2)
68  {
69  index=4*_order+j-10;
70  return _values[index];
71  }
72  else
73  {
74  throw CException("Invalid index for Comet matrix");
75  return _values[0];
76  }
77  }
CometMatrix::_values
TArray _values
Definition: CometMatrix.h:196
template<class T >
T CometMatrix< T >::getTraceAbs ( )
inline

Definition at line 186 of file CometMatrix.h.

References CometMatrix< T >::_order, CometMatrix< T >::_values, and fabs().

187  {
188  T trace=0.;
189  for(int i=0;i<_order;i++)
190  trace+=fabs(_values[i]);
191  return trace;
192  }
CometMatrix::_values
TArray _values
Definition: CometMatrix.h:196
fabs
Tangent fabs(const Tangent &a)
Definition: Tangent.h:312
template<class T >
void CometMatrix< T >::printElement ( const int &  i,
const int &  j 
)
inline

Definition at line 79 of file CometMatrix.h.

References CometMatrix< T >::_order, and CometMatrix< T >::_values.

80  {
81  int index;
82  if(i==j)
83  {
84  index=(i-1);
85  cout<<_values[index]<<endl;
86  }
87  else if(j==_order)
88  {
89  index=i-1+_order;
90  cout<<_values[index]<<endl;
91  }
92  else if(i==_order)
93  {
94  index=2*_order+j-2;
95  cout<<_values[index]<<endl;
96  }
97  else
98  {
99  throw CException("Invalid index for Comet matrix");
100  }
101  }
CometMatrix::_values
TArray _values
Definition: CometMatrix.h:196
template<class T >
void CometMatrix< T >::Solve ( TArray &  bVec)
inline

Definition at line 103 of file CometMatrix.h.

References CometMatrix< T >::_AMat, CometMatrix< T >::_order, CometMatrix< T >::_vel, and CometMatrix< T >::getElement().

104  {
105 
106  //Replaces bVec with solution vector.
107 
108  T an1i;
109  T an2i;
110  T an3i;
111  T ain1;
112  T ain2;
113  T ain3;
114  T aii;
115  T alpha0;
116  T alpha1;
117  T alpha2;
118 
119  alpha0=0.;
120  alpha1=0.;
121  alpha2=0.;
122 
123  //#pragma omp parallel for default(shared) private(i,an1i,an2i,an3i,aii,ain1,ain2,ain3)
124  {
125  for(int i=1;i<_order-2;i++)
126  {
127  an1i=getElement(_order-2,i);
128  an2i=getElement(_order-1,i);
129  an3i=getElement(_order,i);
130  aii=getElement(i,i);
131  ain1=getElement(i,_order-2);
132  ain2=getElement(i,_order-1);
133  ain3=getElement(i,_order);
134 
135  _AMat.getElement(1,1)-=an1i*ain1/aii;
136  _AMat.getElement(1,2)-=an1i*ain2/aii;
137  _AMat.getElement(1,3)-=an1i*ain3/aii;
138 
139  _AMat.getElement(2,1)-=an2i*ain1/aii;
140  _AMat.getElement(2,2)-=an2i*ain2/aii;
141  _AMat.getElement(2,3)-=an2i*ain3/aii;
142 
143  _AMat.getElement(3,1)-=an3i*ain1/aii;
144  _AMat.getElement(3,2)-=an3i*ain2/aii;
145  _AMat.getElement(3,3)-=an3i*ain3/aii;
146 
147  alpha0+=an1i*bVec[i-1]/aii;
148  alpha1+=an2i*bVec[i-1]/aii;
149  alpha2+=an3i*bVec[i-1]/aii;
150  }
151  }
152 
153  _AMat.getElement(1,1)+=getElement(_order-2,_order-2);
154  _AMat.getElement(2,2)+=getElement(_order-1,_order-1);
155  _AMat.getElement(3,3)+=getElement(_order,_order);
156 
157  _vel[0]=(bVec[_order-3]-alpha0);
158  _vel[1]=(bVec[_order-2]-alpha1);
159  _vel[2]=(bVec[_order-1]-alpha2);
160 
161  _AMat.Solve(_vel);
162 
163  bVec[_order-3]=_vel[0];
164  bVec[_order-2]=_vel[1];
165  bVec[_order-1]=_vel[2];
166 
167  //#pragma omp parallel for default(shared) private(i,aii,ain1,ain2,ain3)
168  {
169  for(int i=1;i<_order-2;i++)
170  {
171  ain1=getElement(i,_order-2);
172  ain2=getElement(i,_order-1);
173  ain3=getElement(i,_order);
174  aii=getElement(i,i);
175  bVec[i-1]=(bVec[i-1]-ain1*bVec[_order-3]-ain2*bVec[_order-2]-ain3*bVec[_order-1])/aii;
176  }
177  }
178  }
CometMatrix::_AMat
SquareMatrixESBGK< T > _AMat
Definition: CometMatrix.h:198
CometMatrix::_vel
TArray _vel
Definition: CometMatrix.h:199
CometMatrix::getElement
T & getElement(const int i, const int j)
Definition: CometMatrix.h:34
template<class T >
void CometMatrix< T >::zero ( )
inlinevirtual

Implements MatrixJML< T >.

Definition at line 180 of file CometMatrix.h.

References CometMatrix< T >::_elements, and CometMatrix< T >::_values.

181  {
182  for(int i=0;i<_elements;i++)
183  _values[i]=T_Scalar(0);
184  }
CometMatrix::_values
TArray _values
Definition: CometMatrix.h:196
CometMatrix::T_Scalar
NumTypeTraits< T >::T_Scalar T_Scalar
Definition: CometMatrix.h:20

Member Data Documentation

template<class T >
SquareMatrixESBGK<T> CometMatrix< T >::_AMat
private

Definition at line 198 of file CometMatrix.h.

Referenced by CometMatrix< T >::CometMatrix(), and CometMatrix< T >::Solve().

template<class T >
int CometMatrix< T >::_elements
private

Definition at line 195 of file CometMatrix.h.

Referenced by CometMatrix< T >::zero().

template<class T >
int CometMatrix< T >::_order
private

Definition at line 197 of file CometMatrix.h.

Referenced by CometMatrix< T >::getElement(), CometMatrix< T >::getTraceAbs(), CometMatrix< T >::printElement(), and CometMatrix< T >::Solve().

template<class T >
TArray CometMatrix< T >::_values
private

Definition at line 196 of file CometMatrix.h.

Referenced by CometMatrix< T >::CometMatrix(), CometMatrix< T >::getElement(), CometMatrix< T >::getTraceAbs(), CometMatrix< T >::printElement(), and CometMatrix< T >::zero().

template<class T >
TArray CometMatrix< T >::_vel
private

Definition at line 199 of file CometMatrix.h.

Referenced by CometMatrix< T >::CometMatrix(), and CometMatrix< T >::Solve().

The documentation for this class was generated from the following file: