queso-0.53.0
Public Member Functions | Public Attributes | List of all members
QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M > Class Template Reference

#include <GcmJointTildeInfo.h>

Collaboration diagram for QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >:
Collaboration graph
[legend]

Public Member Functions

 GcmJointTildeInfo (const GpmsaComputerModelOptions &gcmOptionsObj, const GcmExperimentInfo< S_V, S_M, D_V, D_M, P_V, P_M > &e, const GcmJointInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M > &jj)
 
 ~GcmJointTildeInfo ()
 

Public Attributes

const BaseEnvironment & m_env
 
D_M m_Bmat_tilde
 
unsigned int m_Bmat_tilde_rank
 
VectorSpace< D_V, D_M > m_vu_tilde_space
 
D_M m_Lbmat
 
D_M m_Btildet_Wy_Btilde
 
D_M m_Btildet_Wy_Btilde_inv
 
D_V m_Zvec_tilde_hat_vu
 
unsigned int m_a_y_modifier_tilde
 
unsigned int m_b_y_modifier_tilde
 

Detailed Description

template<class S_V = GslVector, class S_M = GslMatrix, class D_V = GslVector, class D_M = GslMatrix, class P_V = GslVector, class P_M = GslMatrix, class Q_V = GslVector, class Q_M = GslMatrix>
class QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >

Definition at line 40 of file GcmJointTildeInfo.h.

Constructor & Destructor Documentation

template<class S_V , class S_M , class D_V , class D_M , class P_V , class P_M , class Q_V , class Q_M >
QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::GcmJointTildeInfo ( const GpmsaComputerModelOptions &  gcmOptionsObj,
const GcmExperimentInfo< S_V, S_M, D_V, D_M, P_V, P_M > &  e,
const GcmJointInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M > &  jj 
)

Definition at line 32 of file GcmJointTildeInfo.C.

References QUESO::BaseEnvironment::checkingLevel(), QUESO::VectorSpace< V, M >::dimGlobal(), QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::m_a_y_modifier_tilde, QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::m_b_y_modifier_tilde, QUESO::GcmJointInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::m_Bmat_rank, QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::m_Bmat_tilde, QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::m_Bmat_tilde_rank, QUESO::GcmJointInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::m_Bmat_with_permut, QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::m_Btildet_Wy_Btilde, QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::m_Btildet_Wy_Btilde_inv, QUESO::GcmOptionsValues::m_dataOutputAllowedSet, QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::m_env, QUESO::GcmExperimentInfo< S_V, S_M, D_V, D_M, P_V, P_M >::m_experimentStorage, QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::m_Lbmat, QUESO::GpmsaComputerModelOptions::m_ov, QUESO::GcmExperimentInfo< S_V, S_M, D_V, D_M, P_V, P_M >::m_paper_n_y, QUESO::GcmJointInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::m_vu_size, QUESO::GcmJointInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::m_vu_space, QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::m_vu_tilde_space, QUESO::GcmExperimentInfo< S_V, S_M, D_V, D_M, P_V, P_M >::m_Wy, QUESO::GcmExperimentInfo< S_V, S_M, D_V, D_M, P_V, P_M >::m_y_space, QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::m_Zvec_tilde_hat_vu, QUESO::VectorSpace< V, M >::map(), queso_require_equal_to_msg, QUESO::scalarProduct(), QUESO::BaseEnvironment::subDisplayFile(), QUESO::BaseEnvironment::subId(), and QUESO::VectorSpace< V, M >::zeroVector().

36  :
37  m_env (jj.m_env),
38  m_Bmat_tilde (m_env,e.m_y_space.map(),jj.m_Bmat_rank),
39  m_Bmat_tilde_rank (0),
40  m_vu_tilde_space (m_env, "vu_tilde_", jj.m_Bmat_rank, NULL), // rr0 check
41  m_Lbmat (m_env,m_vu_tilde_space.map(),jj.m_Bmat_with_permut->numCols()), // rr0 check
42  m_Btildet_Wy_Btilde (m_vu_tilde_space.zeroVector()),
43  m_Btildet_Wy_Btilde_inv(m_vu_tilde_space.zeroVector()),
44  m_Zvec_tilde_hat_vu (m_vu_tilde_space.zeroVector()),
45  m_a_y_modifier_tilde (0),
46  m_b_y_modifier_tilde (0)
47 {
48  std::set<unsigned int> tmpSet;
49  tmpSet.insert(m_env.subId());
50  //******************************************************************************
51  // Tilde situation: form 'm_Bmat_tilde'
52  // Tilde situation: form 'm_vu_tilde_space'
53  // Tilde situation: form 'm_Lbmat'
54  //******************************************************************************
55  if (jj.m_Bmat_with_permut->numRowsGlobal() >= jj.m_Bmat_with_permut->numCols()) {
56  D_M matbU(m_env,e.m_y_space.map(),jj.m_vu_size); // same as m_Bmat
57  matbU = jj.m_Bmat_with_permut->svdMatU();
58  unsigned int buMatRank = matbU.rank(0.,1.e-8 ); // todo: should be an option
59  unsigned int buMatRank14 = matbU.rank(0.,1.e-14);
60  if (m_env.subDisplayFile()) {
61  *m_env.subDisplayFile() << "In GcmJointTildeInfo<S_V,S_M,D_V,D_M,P_V,P_M,Q_V,Q_M>::constructor()"
62  << ": matbU.numRowsLocal() = " << matbU.numRowsLocal()
63  << ", matbU.numCols() = " << matbU.numCols()
64  << ", matbU.rank(0.,1.e-8) = " << buMatRank
65  << ", matbU.rank(0.,1.e-14) = " << buMatRank14
66  << std::endl;
67  }
68 
69  if (m_env.checkingLevel() >= 1) {
70  D_M matbUcheck(jj.m_vu_space.zeroVector());
71  D_V vecI(e.m_y_space.zeroVector());
72  D_V vecJ(e.m_y_space.zeroVector());
73  for (unsigned int i = 0; i < matbU.numCols(); ++i) {
74  matbU.getColumn(i,vecI);
75  for (unsigned int j = i; j < matbU.numCols(); ++j) {
76  matbU.getColumn(j,vecJ);
77  matbUcheck(i,j) = scalarProduct(vecI,vecJ);
78  }
79  }
80  matbUcheck.setPrintHorizontally(false);
81  if (m_env.subDisplayFile()) {
82  *m_env.subDisplayFile() << "In GcmJointTildeInfo<S_V,S_M,D_V,D_M,P_V,P_M,Q_V,Q_M>::constructor()"
83  << ": m_Bmat_with_permut->numRowsLocal() = " << jj.m_Bmat_with_permut->numRowsLocal()
84  << ", m_Bmat_with_permut->numCols() = " << jj.m_Bmat_with_permut->numCols()
85  << ", m_Bmat_rank = " << jj.m_Bmat_rank
86  << ", matbU.numRowsLocal() = " << matbU.numRowsLocal()
87  << ", matbU.numCols() = " << matbU.numCols()
88  << ", buMatrank(0.,1.e-8) = " << buMatRank
89  << ", buMatrank(0.,1.e-14) = " << buMatRank14
90  << ", matbUcheck.numRowsLocal() = " << matbUcheck.numRowsLocal()
91  << ", matbUcheck.numCols() = " << matbUcheck.numCols()
92  << ", matbUcheck =\n" << matbUcheck
93  << std::endl;
94  }
95  }
96 
97  D_V vecbJ(e.m_y_space.zeroVector());
98  for (unsigned int j = 0; j < jj.m_Bmat_rank; ++j) {
99  matbU.getColumn(j,vecbJ);
100  m_Bmat_tilde.setColumn(j,vecbJ);
101  }
102  if (gcmOptionsObj.m_ov.m_dataOutputAllowedSet.find(m_env.subId()) != gcmOptionsObj.m_ov.m_dataOutputAllowedSet.end()) {
103  m_Bmat_tilde.subWriteContents("Btilde",
104  "Btilde1",
105  "m",
106  tmpSet);
107  }
108  if (m_env.subDisplayFile()) {
109  *m_env.subDisplayFile() << "In GcmJointTildeInfo<S_V,S_M,D_V,D_M,P_V,P_M,Q_V,Q_M>::constructor()"
110  << ": m_Bmat_tilde computed (1)"
111  << std::endl;
112  }
113  }
114  else {
115  D_M Bmat_t(m_env,jj.m_vu_space.map(),e.m_paper_n_y); // same as m_Bmat^T
116  Bmat_t.fillWithTranspose(0,0,*jj.m_Bmat_with_permut,true,true);
117 
118  D_M matbV(e.m_y_space.zeroVector());
119  matbV = Bmat_t.svdMatV();
120  unsigned int bvMatRank = matbV.rank(0.,1.e-8); // todo: should be an option
121  unsigned int bvMatRank14 = matbV.rank(0.,1.e-14);
122  if (m_env.subDisplayFile()) {
123  *m_env.subDisplayFile() << "In GcmJointTildeInfo<S_V,S_M,D_V,D_M,P_V,P_M,Q_V,Q_M>::constructor()"
124  << ": matbV.numRowsLocal() = " << matbV.numRowsLocal()
125  << ", matbV.numCols() = " << matbV.numCols()
126  << ", matbV.rank(0.,1.e-8) = " << bvMatRank
127  << ", matbV.rank(0.,1.e-14) = " << bvMatRank14
128  << std::endl;
129  }
130 
131  if (m_env.checkingLevel() >= 1) {
132  D_M matbVcheck(e.m_y_space.zeroVector());
133  D_V vecI(e.m_y_space.zeroVector());
134  D_V vecJ(e.m_y_space.zeroVector());
135  for (unsigned int i = 0; i < matbV.numCols(); ++i) {
136  matbV.getColumn(i,vecI);
137  for (unsigned int j = i; j < matbV.numCols(); ++j) {
138  matbV.getColumn(j,vecJ);
139  matbVcheck(i,j) = scalarProduct(vecI,vecJ);
140  }
141  }
142  matbVcheck.setPrintHorizontally(false);
143  if (m_env.subDisplayFile()) {
144  *m_env.subDisplayFile() << "In GcmJointTildeInfo<S_V,S_M,D_V,D_M,P_V,P_M,Q_V,Q_M>::constructor()"
145  << ": m_Bmat_with_permut->numRowsLocal() = " << jj.m_Bmat_with_permut->numRowsLocal()
146  << ", m_Bmat_with_permut->numCols() = " << jj.m_Bmat_with_permut->numCols()
147  << ", m_Bmat_rank = " << jj.m_Bmat_rank
148  << ", matbV.numRowsLocal() = " << matbV.numRowsLocal()
149  << ", matbV.numCols() = " << matbV.numCols()
150  << ", bvMatrank(0.,1.e-8) = " << bvMatRank
151  << ", bvMatrank(0.,1.e-14) = " << bvMatRank14
152  << ", matbVcheck.numRowsLocal() = " << matbVcheck.numRowsLocal()
153  << ", matbVcheck.numCols() = " << matbVcheck.numCols()
154  << ", matbVcheck =\n" << matbVcheck
155  << std::endl;
156  }
157  }
158 
159  D_V vecbJ(e.m_y_space.zeroVector());
160  for (unsigned int j = 0; j < jj.m_Bmat_rank; ++j) {
161  matbV.getColumn(j,vecbJ);
162  m_Bmat_tilde.setColumn(j,vecbJ);
163  }
164  if (gcmOptionsObj.m_ov.m_dataOutputAllowedSet.find(m_env.subId()) != gcmOptionsObj.m_ov.m_dataOutputAllowedSet.end()) {
165  m_Bmat_tilde.subWriteContents("Btilde",
166  "Btilde2",
167  "m",
168  tmpSet);
169  }
170  if (m_env.subDisplayFile()) {
171  *m_env.subDisplayFile() << "In GcmJointTildeInfo<S_V,S_M,D_V,D_M,P_V,P_M,Q_V,Q_M>::constructor()"
172  << ": m_Bmat_tilde computed (2)"
173  << std::endl;
174  }
175  }
176 
177  if (m_env.subDisplayFile()) {
178  *m_env.subDisplayFile() << "In GcmJointTildeInfo<S_V,S_M,D_V,D_M,P_V,P_M,Q_V,Q_M>::constructor()"
179  << ": finished forming 'm_Bmat_tilde'"
180  << ", m_Bmat_tilde.numRowsLocal() = " << m_Bmat_tilde.numRowsLocal()
181  << ", m_Bmat_tilde.numCols() = " << m_Bmat_tilde.numCols()
182  << std::endl;
183  }
184 
185  m_Bmat_tilde.svdSolve(*jj.m_Bmat_with_permut,m_Lbmat);
186  if (gcmOptionsObj.m_ov.m_dataOutputAllowedSet.find(m_env.subId()) != gcmOptionsObj.m_ov.m_dataOutputAllowedSet.end()) {
187  m_Lbmat.subWriteContents("Lbmat",
188  "Lbmat",
189  "m",
190  tmpSet);
191  }
192  if (m_env.subDisplayFile()) {
193  *m_env.subDisplayFile() << "In GcmJointTildeInfo<S_V,S_M,D_V,D_M,P_V,P_M,Q_V,Q_M>::constructor()"
194  << ": m_Lbmat_tilde computed"
195  << std::endl;
196  }
197 
198  //********************************************************************************
199  // Form 'Btilde^T' matrix
200  //********************************************************************************
201  D_M Btildet(m_env,m_vu_tilde_space.map(),e.m_paper_n_y);
202  Btildet.fillWithTranspose(0,0,m_Bmat_tilde,true,true);
203 
204  if (m_env.checkingLevel() >= 1) {
205  // Check transpose operation
206  D_M Btildett(m_env,e.m_y_space.map(),m_vu_tilde_space.dimGlobal());
207  Btildett.fillWithTranspose(0,0,Btildet,true,true);
208  Btildett -= m_Bmat_tilde;
209  double btDiffNorm = Btildett.normFrob();
210  if (m_env.subDisplayFile()) {
211  *m_env.subDisplayFile() << "In GcmJointTildeInfo<S_V,S_M,D_V,D_M,P_V,P_M,Q_V,Q_M>::constructor()"
212  << ": ||Btildett - Btilde||_2 = " << btDiffNorm
213  << std::endl;
214  }
215  }
216 
217  //********************************************************************************
218  // Compute 'Btilde' rank
219  //********************************************************************************
220  double bTildeRank14 = 0.;
221  if (m_Bmat_tilde.numRowsGlobal() >= m_Bmat_tilde.numCols()) {
222  m_Bmat_tilde_rank = m_Bmat_tilde.rank(0.,1.e-8 ); // todo: should be an option
223  bTildeRank14 = m_Bmat_tilde.rank(0.,1.e-14);
224  }
225  else {
226  m_Bmat_tilde_rank = Btildet.rank(0.,1.e-8 ); // todo: should be an option
227  bTildeRank14 = Btildet.rank(0.,1.e-14);
228  }
229 
230  if (m_env.subDisplayFile()) {
231  *m_env.subDisplayFile() << "In GcmJointTildeInfo<S_V,S_M,D_V,D_M,P_V,P_M,Q_V,Q_M>::constructor()"
232  << ": m_Bmat_tilde.numRowsLocal() = " << m_Bmat_tilde.numRowsLocal()
233  << ", m_Bmat_tilde.numCols() = " << m_Bmat_tilde.numCols()
234  << ", m_Bmat_tilde.rank(0.,1.e-8) = " << m_Bmat_tilde_rank
235  << ", m_Bmat_tilde.rank(0.,1.e-14) = " << bTildeRank14
236  << std::endl;
237  }
238  queso_require_equal_to_msg(m_Bmat_tilde_rank, std::min(m_Bmat_tilde.numRowsGlobal(),m_Bmat_tilde.numCols()), "'m_Bmat_tilde' does not have a proper rank");
239 
240  //******************************************************************************
241  // Tilde situation: compute 'Btilde^T W_y Btilde' matrix, and its inverse
242  //******************************************************************************
243  m_Btildet_Wy_Btilde = Btildet * (*e.m_Wy * m_Bmat_tilde); // todo: add 1.e-4 to diagonal
244  if (m_env.subDisplayFile()) {
245  *m_env.subDisplayFile() << "In GcmJointTildeInfo<S_V,S_M,D_V,D_M,P_V,P_M,Q_V,Q_M>::constructor()"
246  << ": finished computing 'm_Btildet_Wy_Btilde'"
247  << std::endl;
248  }
249  if (gcmOptionsObj.m_ov.m_dataOutputAllowedSet.find(m_env.subId()) != gcmOptionsObj.m_ov.m_dataOutputAllowedSet.end()) {
250  m_Btildet_Wy_Btilde.subWriteContents("Btildet_Wy_Btilde",
251  "Btildet_Wy_Btilde",
252  "m",
253  tmpSet);
254  }
255  if (m_env.subDisplayFile()) {
256  *m_env.subDisplayFile() << "In GcmJointTildeInfo<S_V,S_M,D_V,D_M,P_V,P_M,Q_V,Q_M>::constructor()"
257  << ": m_Btildet_Wy_Btilde computed"
258  << std::endl;
259  }
260 
261  double btildetWyBtildeLnDeterminant = m_Btildet_Wy_Btilde.lnDeterminant();
262  unsigned int btildetWyBtildeRank = m_Btildet_Wy_Btilde.rank(0.,1.e-8 ); // todo: should be an option
263  unsigned int btildetWyBtildeRank14 = m_Btildet_Wy_Btilde.rank(0.,1.e-14);
264  if (m_env.subDisplayFile()) {
265  *m_env.subDisplayFile() << "In GcmJointTildeInfo<S_V,S_M,D_V,D_M,P_V,P_M,Q_V,Q_M>::constructor()"
266  << ": m_Btildet_Wy_Btilde.numRowsLocal() = " << m_Btildet_Wy_Btilde.numRowsLocal()
267  << ", m_Btildet_Wy_Btilde.numCols() = " << m_Btildet_Wy_Btilde.numCols()
268  << ", m_Btildet_Wy_Btilde.lnDeterminant() = " << btildetWyBtildeLnDeterminant
269  << ": m_Btildet_Wy_Btilde.rank(0.,1.e-8) = " << btildetWyBtildeRank
270  << ": m_Btildet_Wy_Btilde.rank(0.,1.e-14) = " << btildetWyBtildeRank14
271  << std::endl;
272  }
273 
274  m_Btildet_Wy_Btilde_inv = m_Btildet_Wy_Btilde.inverse(); // todo: add 1.e-6 to diagonal
275  if (m_env.subDisplayFile()) {
276  *m_env.subDisplayFile() << "In GcmJointTildeInfo<S_V,S_M,D_V,D_M,P_V,P_M,Q_V,Q_M>::constructor()"
277  << ": finished computing 'm_Btildet_Wy_Btilde_inv'"
278  << ", m_Btildet_Wy_Btilde_inv.lnDeterminant() = " << m_Btildet_Wy_Btilde_inv.lnDeterminant()
279  << std::endl;
280  }
281  if (gcmOptionsObj.m_ov.m_dataOutputAllowedSet.find(m_env.subId()) != gcmOptionsObj.m_ov.m_dataOutputAllowedSet.end()) {
282  m_Btildet_Wy_Btilde_inv.subWriteContents("Btildet_Wy_Btilde_inv",
283  "Btildet_Wy_Btilde_inv",
284  "m",
285  tmpSet);
286  }
287  if (m_env.subDisplayFile()) {
288  *m_env.subDisplayFile() << "In GcmJointTildeInfo<S_V,S_M,D_V,D_M,P_V,P_M,Q_V,Q_M>::constructor()"
289  << ": m_Btildet_Wy_Btilde_inv computed"
290  << std::endl;
291  }
292 
293  double btildetWyBtildeInvLnDeterminant = m_Btildet_Wy_Btilde_inv.lnDeterminant();
294  unsigned int btildetWyBtildeInvRank = m_Btildet_Wy_Btilde_inv.rank(0.,1.e-8 ); // todo: should be an option
295  unsigned int btildetWyBtildeInvRank14 = m_Btildet_Wy_Btilde_inv.rank(0.,1.e-14);
296  if (m_env.subDisplayFile()) {
297  *m_env.subDisplayFile() << "In GcmJointTildeInfo<S_V,S_M,D_V,D_M,P_V,P_M,Q_V,Q_M>::constructor()"
298  << ": m_Btildet_Wy_Btilde_inv.numRowsLocal() = " << m_Btildet_Wy_Btilde_inv.numRowsLocal()
299  << ", m_Btildet_Wy_Btilde_inv.numCols() = " << m_Btildet_Wy_Btilde_inv.numCols()
300  << ": m_Btildet_Wy_Btilde_inv.lnDeterminant() = " << btildetWyBtildeInvLnDeterminant
301  << ": m_Btildet_Wy_Btilde_inv.rank(0.,1.e-8) = " << btildetWyBtildeInvRank
302  << ": m_Btildet_Wy_Btilde_inv.rank(0.,1.e-14) = " << btildetWyBtildeInvRank14
303  << std::endl;
304  }
305 
306  //********************************************************************************
307  // Compute 'tilde' exponent modifiers
308  //********************************************************************************
309  m_a_y_modifier_tilde = ((double) (e.m_paper_n_y - m_Bmat_tilde_rank)) / 2.;
310  if (m_env.subDisplayFile()) {
311  *m_env.subDisplayFile() << "In GcmJointTildeInfo<S_V,S_M,D_V,D_M,P_V,P_M,Q_V,Q_M>::constructor()"
312  << ": m_a_y_modifier_tilde = " << m_a_y_modifier_tilde
313  << std::endl;
314  }
315 
316  D_V yVec_transformed(e.m_experimentStorage.yVec_transformed());
317  if (m_env.subDisplayFile()) {
318  *m_env.subDisplayFile() << "In GcmJointTildeInfo<S_V,S_M,D_V,D_M,P_V,P_M,Q_V,Q_M>::constructor()"
319  << ": m_Zvec_tilde_hat_vu.sizeLocal() = " << m_Zvec_tilde_hat_vu.sizeLocal()
320  << ", m_Btildet_Wy_Btilde_inv.numRowsLocal() = " << m_Btildet_Wy_Btilde_inv.numRowsLocal()
321  << ", m_Btildet_Wy_Btilde_inv.numCols() = " << m_Btildet_Wy_Btilde_inv.numCols()
322  << ", Btildet.numRowsLocal() = " << Btildet.numRowsLocal()
323  << ", Btildet.numCols() = " << Btildet.numCols()
324  << ", m_Wy->numRowsLocal() = " << e.m_Wy->numRowsLocal()
325  << ", m_Wy->numCols() = " << e.m_Wy->numCols()
326  << ", yVec_transformed.sizeLocal() = " << yVec_transformed.sizeLocal()
327  << std::endl;
328  }
329  m_Zvec_tilde_hat_vu = m_Btildet_Wy_Btilde_inv * (Btildet * (*e.m_Wy * yVec_transformed));
330  D_V tmpVec2(yVec_transformed - (m_Bmat_tilde * m_Zvec_tilde_hat_vu));
331  tmpVec2 = *e.m_Wy * tmpVec2;
332  m_b_y_modifier_tilde = scalarProduct(yVec_transformed,tmpVec2) / 2.;
333 
334  if (m_env.subDisplayFile()) {
335  *m_env.subDisplayFile() << "In GcmJointTildeInfo<S_V,S_M,D_V,D_M,P_V,P_M,Q_V,Q_M>::constructor()"
336  << ": m_b_y_modifier_tilde = " << m_b_y_modifier_tilde
337  << std::endl;
338  }
339 
340  if (m_env.subDisplayFile()) {
341  *m_env.subDisplayFile() << "In GcmJointTildeInfo<S_V,S_M,D_V,D_M,P_V,P_M,Q_V,Q_M>::constructor()"
342  << ": finished computing 'tilde' exponent modifiers"
343  << std::endl;
344  }
345 }
QUESO::GcmJointTildeInfo::m_env
const BaseEnvironment & m_env
Definition: GcmJointTildeInfo.h:48
QUESO::VectorSpace::map
const Map & map() const
Map.
Definition: VectorSpace.C:157
QUESO::GcmJointTildeInfo::m_vu_tilde_space
VectorSpace< D_V, D_M > m_vu_tilde_space
Definition: GcmJointTildeInfo.h:51
QUESO::VectorSpace::zeroVector
const V & zeroVector() const
Returns a vector filled with zeros.
Definition: VectorSpace.C:189
QUESO::VectorSpace::dimGlobal
unsigned int dimGlobal() const
Definition: VectorSpace.C:176
queso_require_equal_to_msg
#define queso_require_equal_to_msg(expr1, expr2, msg)
Definition: asserts.h:85
QUESO::scalarProduct
double scalarProduct(const GslVector &x, const GslVector &y)
Definition: GslVector.C:1150
QUESO::BaseEnvironment::subDisplayFile
std::ofstream * subDisplayFile() const
Access function for m_subDisplayFile (displays file on stream).
Definition: Environment.C:274
QUESO::BaseEnvironment::subId
unsigned int subId() const
Access function to the number of each sub-environment Id: m_subId.
Definition: Environment.C:295
QUESO::BaseEnvironment::checkingLevel
unsigned int checkingLevel() const
Access function to private attribute m_checkingLevel.
Definition: Environment.C:410
template<class S_V , class S_M , class D_V , class D_M , class P_V , class P_M , class Q_V , class Q_M >
QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::~GcmJointTildeInfo ( )

Definition at line 348 of file GcmJointTildeInfo.C.

349 {
350 }

Member Data Documentation

template<class S_V = GslVector, class S_M = GslMatrix, class D_V = GslVector, class D_M = GslMatrix, class P_V = GslVector, class P_M = GslMatrix, class Q_V = GslVector, class Q_M = GslMatrix>
unsigned int QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::m_a_y_modifier_tilde

Definition at line 56 of file GcmJointTildeInfo.h.

Referenced by QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::GcmJointTildeInfo().

template<class S_V = GslVector, class S_M = GslMatrix, class D_V = GslVector, class D_M = GslMatrix, class P_V = GslVector, class P_M = GslMatrix, class Q_V = GslVector, class Q_M = GslMatrix>
unsigned int QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::m_b_y_modifier_tilde

Definition at line 57 of file GcmJointTildeInfo.h.

Referenced by QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::GcmJointTildeInfo().

template<class S_V = GslVector, class S_M = GslMatrix, class D_V = GslVector, class D_M = GslMatrix, class P_V = GslVector, class P_M = GslMatrix, class Q_V = GslVector, class Q_M = GslMatrix>
D_M QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::m_Bmat_tilde

Definition at line 49 of file GcmJointTildeInfo.h.

Referenced by QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::GcmJointTildeInfo(), and QUESO::GcmZTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::GcmZTildeInfo().

template<class S_V = GslVector, class S_M = GslMatrix, class D_V = GslVector, class D_M = GslMatrix, class P_V = GslVector, class P_M = GslMatrix, class Q_V = GslVector, class Q_M = GslMatrix>
unsigned int QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::m_Bmat_tilde_rank

Definition at line 50 of file GcmJointTildeInfo.h.

Referenced by QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::GcmJointTildeInfo().

template<class S_V = GslVector, class S_M = GslMatrix, class D_V = GslVector, class D_M = GslMatrix, class P_V = GslVector, class P_M = GslMatrix, class Q_V = GslVector, class Q_M = GslMatrix>
D_M QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::m_Btildet_Wy_Btilde

Definition at line 53 of file GcmJointTildeInfo.h.

Referenced by QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::GcmJointTildeInfo().

template<class S_V = GslVector, class S_M = GslMatrix, class D_V = GslVector, class D_M = GslMatrix, class P_V = GslVector, class P_M = GslMatrix, class Q_V = GslVector, class Q_M = GslMatrix>
D_M QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::m_Btildet_Wy_Btilde_inv

Definition at line 54 of file GcmJointTildeInfo.h.

Referenced by QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::GcmJointTildeInfo().

template<class S_V = GslVector, class S_M = GslMatrix, class D_V = GslVector, class D_M = GslMatrix, class P_V = GslVector, class P_M = GslMatrix, class Q_V = GslVector, class Q_M = GslMatrix>
const BaseEnvironment& QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::m_env

Definition at line 48 of file GcmJointTildeInfo.h.

Referenced by QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::GcmJointTildeInfo().

template<class S_V = GslVector, class S_M = GslMatrix, class D_V = GslVector, class D_M = GslMatrix, class P_V = GslVector, class P_M = GslMatrix, class Q_V = GslVector, class Q_M = GslMatrix>
D_M QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::m_Lbmat

Definition at line 52 of file GcmJointTildeInfo.h.

Referenced by QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::GcmJointTildeInfo(), and QUESO::GcmZTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::GcmZTildeInfo().

template<class S_V = GslVector, class S_M = GslMatrix, class D_V = GslVector, class D_M = GslMatrix, class P_V = GslVector, class P_M = GslMatrix, class Q_V = GslVector, class Q_M = GslMatrix>
VectorSpace<D_V,D_M> QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::m_vu_tilde_space

Definition at line 51 of file GcmJointTildeInfo.h.

Referenced by QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::GcmJointTildeInfo().

template<class S_V = GslVector, class S_M = GslMatrix, class D_V = GslVector, class D_M = GslMatrix, class P_V = GslVector, class P_M = GslMatrix, class Q_V = GslVector, class Q_M = GslMatrix>
D_V QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::m_Zvec_tilde_hat_vu

Definition at line 55 of file GcmJointTildeInfo.h.

Referenced by QUESO::GcmJointTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::GcmJointTildeInfo(), and QUESO::GcmZTildeInfo< S_V, S_M, D_V, D_M, P_V, P_M, Q_V, Q_M >::GcmZTildeInfo().

The documentation for this class was generated from the following files:
Generated on Thu Jun 11 2015 13:52:34 for queso-0.53.0 by  doxygen 1.8.5