36 |
|
* [1] Meineke, et al., J. Comp. Chem. 26, 252-271 (2005). |
37 |
|
* [2] Fennell & Gezelter, J. Chem. Phys. 124, 234104 (2006). |
38 |
|
* [3] Sun, Lin & Gezelter, J. Chem. Phys. 128, 24107 (2008). |
39 |
< |
* [4] Vardeman & Gezelter, in progress (2009). |
39 |
> |
* [4] Kuang & Gezelter, J. Chem. Phys. 133, 164101 (2010). |
40 |
> |
* [5] Vardeman, Stocker & Gezelter, J. Chem. Theory Comput. 7, 834 (2011). |
41 |
|
*/ |
42 |
|
|
43 |
|
/** |
133 |
|
for (unsigned int i = 0; i < Dim ; i++) |
134 |
|
tmp += this->data_[i][i]; |
135 |
|
|
136 |
+ |
return tmp; |
137 |
+ |
} |
138 |
+ |
|
139 |
+ |
/** |
140 |
+ |
* Returns the tensor contraction (double dot product) of two rank 2 |
141 |
+ |
* tensors (or Matrices) |
142 |
+ |
* @param t1 first tensor |
143 |
+ |
* @param t2 second tensor |
144 |
+ |
* @return the tensor contraction (double dot product) of t1 and t2 |
145 |
+ |
*/ |
146 |
+ |
Real doubleDot( const SquareMatrix<Real, Dim>& t1, const SquareMatrix<Real, Dim>& t2 ) { |
147 |
+ |
Real tmp; |
148 |
+ |
tmp = 0; |
149 |
+ |
|
150 |
+ |
for (unsigned int i = 0; i < Dim; i++) |
151 |
+ |
for (unsigned int j =0; j < Dim; j++) |
152 |
+ |
tmp += t1[i][j] * t2[i][j]; |
153 |
+ |
|
154 |
|
return tmp; |
155 |
|
} |
156 |
|
|
181 |
|
return false; |
182 |
|
|
183 |
|
return true; |
184 |
+ |
} |
185 |
+ |
|
186 |
+ |
/** |
187 |
+ |
* Returns a column vector that contains the elements from the |
188 |
+ |
* diagonal of m in the order R(0) = m(0,0), R(1) = m(1,1), and so |
189 |
+ |
* on. |
190 |
+ |
*/ |
191 |
+ |
Vector<Real, Dim> diagonals() const { |
192 |
+ |
Vector<Real, Dim> result; |
193 |
+ |
for (unsigned int i = 0; i < Dim; i++) { |
194 |
+ |
result(i) = this->data_[i][i]; |
195 |
+ |
} |
196 |
+ |
return result; |
197 |
|
} |
198 |
|
|
199 |
|
/** Tests if this matrix is the unit matrix. */ |