151 std::vector<int> tags( tagSize * this->getCardinality() );
154 for (
int i=0;i<this->getCardinality();i++) {
156 tags[tagSize*i+1] = 0;
157 tags[tagSize*i+2] = i;
158 tags[tagSize*i+3] = this->getCardinality();
161 Basis<Scalar,ArrayScalar> &xBasis_ = *this->bases_[0][0];
162 Basis<Scalar,ArrayScalar> &yBasis_ = *this->bases_[0][1];
163 Basis<Scalar,ArrayScalar> &zBasis_ = *this->bases_[0][2];
168 const std::vector<std::vector<int> >& xdoftags = xBasis_.getAllDofTags();
169 const std::vector<std::vector<int> >& ydoftags = yBasis_.getAllDofTags();
170 const std::vector<std::vector<int> >& zdoftags = zBasis_.getAllDofTags();
172 std::map<int,std::map<int,int> > total_dof_per_entity;
173 std::map<int,std::map<int,int> > current_dof_per_entity;
176 for (
int i=0;i<8;i++) {
177 total_dof_per_entity[0][i] = 0;
178 current_dof_per_entity[0][i] = 0;
181 for (
int i=0;i<12;i++) {
182 total_dof_per_entity[1][i] = 0;
183 current_dof_per_entity[1][i] = 0;
186 for (
int i=0;i<6;i++) {
187 total_dof_per_entity[2][i] = 0;
188 current_dof_per_entity[2][i] = 0;
191 total_dof_per_entity[3][0] = 0;
195 for (
int k=0;k<zBasis_.getCardinality();k++) {
196 const int zdim = zdoftags[k][0];
197 const int zent = zdoftags[k][1];
198 for (
int j=0;j<yBasis_.getCardinality();j++) {
199 const int ydim = ydoftags[j][0];
200 const int yent = ydoftags[j][1];
201 for (
int i=0;i<xBasis_.getCardinality();i++) {
202 const int xdim = xdoftags[i][0];
203 const int xent = xdoftags[i][1];
207 total_dof_per_entity[dofdim][dofent] += 1;
214 for (
int k=0;k<zBasis_.getCardinality();k++) {
215 const int zdim = zdoftags[k][0];
216 const int zent = zdoftags[k][1];
217 for (
int j=0;j<yBasis_.getCardinality();j++) {
218 const int ydim = ydoftags[j][0];
219 const int yent = ydoftags[j][1];
220 for (
int i=0;i<xBasis_.getCardinality();i++) {
221 const int xdim = xdoftags[i][0];
222 const int xent = xdoftags[i][1];
226 tags[4*tagcur] = dofdim;
227 tags[4*tagcur+1] = dofent;
228 tags[4*tagcur+2] = current_dof_per_entity[dofdim][dofent];
229 current_dof_per_entity[dofdim][dofent]++;
230 tags[4*tagcur+3] = total_dof_per_entity[dofdim][dofent];
237 this -> ordinalToTag_,
239 this -> basisCardinality_,
249 const ArrayScalar &inputPoints ,
250 const EOperator operatorType )
const
252#ifdef HAVE_INTREPID_DEBUG
253 Intrepid::getValues_HGRAD_Args<Scalar, ArrayScalar>(outputValues,
256 this -> getBaseCellTopology(),
257 this -> getCardinality() );
260 Basis<Scalar,ArrayScalar> &xBasis_ = *this->bases_[0][0];
261 Basis<Scalar,ArrayScalar> &yBasis_ = *this->bases_[0][1];
262 Basis<Scalar,ArrayScalar> &zBasis_ = *this->bases_[0][2];
265 FieldContainer<Scalar> xInputPoints(inputPoints.dimension(0),1);
266 FieldContainer<Scalar> yInputPoints(inputPoints.dimension(0),1);
267 FieldContainer<Scalar> zInputPoints(inputPoints.dimension(0),1);
269 for (
int i=0;i<inputPoints.dimension(0);i++) {
270 xInputPoints(i,0) = inputPoints(i,0);
271 yInputPoints(i,0) = inputPoints(i,1);
272 zInputPoints(i,0) = inputPoints(i,2);
275 switch (operatorType) {
278 FieldContainer<Scalar> xBasisValues(xBasis_.getCardinality(),xInputPoints.dimension(0));
279 FieldContainer<Scalar> yBasisValues(yBasis_.getCardinality(),yInputPoints.dimension(0));
280 FieldContainer<Scalar> zBasisValues(zBasis_.getCardinality(),zInputPoints.dimension(0));
282 xBasis_.getValues(xBasisValues,xInputPoints,OPERATOR_VALUE);
283 yBasis_.getValues(yBasisValues,yInputPoints,OPERATOR_VALUE);
284 zBasis_.getValues(zBasisValues,zInputPoints,OPERATOR_VALUE);
287 for (
int k=0;k<zBasis_.getCardinality();k++) {
288 for (
int j=0;j<yBasis_.getCardinality();j++) {
289 for (
int i=0;i<xBasis_.getCardinality();i++) {
290 for (
int l=0;l<inputPoints.dimension(0);l++) {
291 outputValues(bfcur,l) = xBasisValues(i,l) * yBasisValues(j,l) * zBasisValues(k,l);
302 FieldContainer<Scalar> xBasisValues(xBasis_.getCardinality(),xInputPoints.dimension(0));
303 FieldContainer<Scalar> yBasisValues(yBasis_.getCardinality(),yInputPoints.dimension(0));
304 FieldContainer<Scalar> zBasisValues(zBasis_.getCardinality(),zInputPoints.dimension(0));
305 FieldContainer<Scalar> xBasisDerivs(xBasis_.getCardinality(),xInputPoints.dimension(0),1);
306 FieldContainer<Scalar> yBasisDerivs(yBasis_.getCardinality(),yInputPoints.dimension(0),1);
307 FieldContainer<Scalar> zBasisDerivs(zBasis_.getCardinality(),zInputPoints.dimension(0),1);
309 xBasis_.getValues(xBasisValues,xInputPoints,OPERATOR_VALUE);
310 yBasis_.getValues(yBasisValues,yInputPoints,OPERATOR_VALUE);
311 zBasis_.getValues(zBasisValues,zInputPoints,OPERATOR_VALUE);
312 xBasis_.getValues(xBasisDerivs,xInputPoints,OPERATOR_D1);
313 yBasis_.getValues(yBasisDerivs,yInputPoints,OPERATOR_D1);
314 zBasis_.getValues(zBasisDerivs,zInputPoints,OPERATOR_D1);
317 for (
int k=0;k<zBasis_.getCardinality();k++) {
318 for (
int j=0;j<yBasis_.getCardinality();j++) {
319 for (
int i=0;i<xBasis_.getCardinality();i++) {
320 for (
int l=0;l<inputPoints.dimension(0);l++) {
321 outputValues(bfcur,l,0) = xBasisDerivs(i,l,0) * yBasisValues(j,l) * zBasisValues(k,l);
322 outputValues(bfcur,l,1) = xBasisValues(i,l) * yBasisDerivs(j,l,0) * zBasisValues(k,l);
323 outputValues(bfcur,l,2) = xBasisValues(i,l) * yBasisValues(j,l) * zBasisDerivs(k,l,0);
341 FieldContainer<Scalar> xBasisValues(xBasis_.getCardinality(),xInputPoints.dimension(0));
342 FieldContainer<Scalar> yBasisValues(yBasis_.getCardinality(),yInputPoints.dimension(0));
343 FieldContainer<Scalar> zBasisValues(yBasis_.getCardinality(),zInputPoints.dimension(0));
345 Teuchos::Array<int> partialMult;
349 if (partialMult[0] == 0) {
350 xBasisValues.resize(xBasis_.getCardinality(),xInputPoints.dimension(0));
351 xBasis_.getValues( xBasisValues , xInputPoints, OPERATOR_VALUE );
354 xBasisValues.resize(xBasis_.getCardinality(),xInputPoints.dimension(0),1);
355 EOperator xop = (EOperator) ( (
int) OPERATOR_D1 + partialMult[0] - 1 );
356 xBasis_.getValues( xBasisValues , xInputPoints, xop );
357 xBasisValues.resize(xBasis_.getCardinality(),xInputPoints.dimension(0));
359 if (partialMult[1] == 0) {
360 yBasisValues.resize(yBasis_.getCardinality(),yInputPoints.dimension(0));
361 yBasis_.getValues( yBasisValues , yInputPoints, OPERATOR_VALUE );
364 yBasisValues.resize(yBasis_.getCardinality(),yInputPoints.dimension(0),1);
365 EOperator yop = (EOperator) ( (
int) OPERATOR_D1 + partialMult[1] - 1 );
366 yBasis_.getValues( yBasisValues , yInputPoints, yop );
367 yBasisValues.resize(yBasis_.getCardinality(),yInputPoints.dimension(0));
369 if (partialMult[2] == 0) {
370 zBasisValues.resize(zBasis_.getCardinality(),zInputPoints.dimension(0));
371 zBasis_.getValues( zBasisValues , zInputPoints, OPERATOR_VALUE );
374 zBasisValues.resize(zBasis_.getCardinality(),zInputPoints.dimension(0),1);
375 EOperator zop = (EOperator) ( (
int) OPERATOR_D1 + partialMult[2] - 1 );
376 zBasis_.getValues( zBasisValues , zInputPoints, zop );
377 zBasisValues.resize(zBasis_.getCardinality(),zInputPoints.dimension(0));
382 for (
int k=0;k<zBasis_.getCardinality();k++) {
383 for (
int j=0;j<yBasis_.getCardinality();j++) {
384 for (
int i=0;i<xBasis_.getCardinality();i++) {
385 for (
int l=0;l<inputPoints.dimension(0);l++) {
386 outputValues(bfcur,l,d) = xBasisValues(i,l) * yBasisValues(j,l) * zBasisValues(k,l);
396 TEUCHOS_TEST_FOR_EXCEPTION(
true , std::invalid_argument,
397 ">>> ERROR (Basis_HGRAD_HEX_Cn_FEM): Operator type not implemented");