136 std::vector<int> tags( tagSize * this->getCardinality() );
139 for (
int i=0;i<this->getCardinality();i++) {
141 tags[tagSize*i+1] = 0;
142 tags[tagSize*i+2] = i;
143 tags[tagSize*i+3] = this->getCardinality();
146 Basis<Scalar,ArrayScalar> &xBasis_ = *this->getBases()[0][0];
147 Basis<Scalar,ArrayScalar> &yBasis_ = *this->getBases()[0][1];
151 const std::vector<std::vector<int> >& xdoftags = xBasis_.getAllDofTags();
152 const std::vector<std::vector<int> >& ydoftags = yBasis_.getAllDofTags();
154 std::map<int,std::map<int,int> > total_dof_per_entity;
155 std::map<int,std::map<int,int> > current_dof_per_entity;
157 for (
int i=0;i<4;i++) {
158 total_dof_per_entity[0][i] = 0;
159 current_dof_per_entity[0][i] = 0;
161 for (
int i=0;i<4;i++) {
162 total_dof_per_entity[1][i] = 0;
163 current_dof_per_entity[1][i] = 0;
165 total_dof_per_entity[2][0] = 0;
166 current_dof_per_entity[2][0] = 0;
169 for (
int j=0;j<yBasis_.getCardinality();j++) {
170 const int ydim = ydoftags[j][0];
171 const int yent = ydoftags[j][1];
172 for (
int i=0;i<xBasis_.getCardinality();i++) {
173 const int xdim = xdoftags[i][0];
174 const int xent = xdoftags[i][1];
178 total_dof_per_entity[dofdim][dofent] += 1;
183 for (
int j=0;j<yBasis_.getCardinality();j++) {
184 const int ydim = ydoftags[j][0];
185 const int yent = ydoftags[j][1];
186 for (
int i=0;i<xBasis_.getCardinality();i++) {
187 const int xdim = xdoftags[i][0];
188 const int xent = xdoftags[i][1];
192 tags[4*tagcur] = dofdim;
193 tags[4*tagcur+1] = dofent;
194 tags[4*tagcur+2] = current_dof_per_entity[dofdim][dofent];
195 current_dof_per_entity[dofdim][dofent]++;
196 tags[4*tagcur+3] = total_dof_per_entity[dofdim][dofent];
202 this -> ordinalToTag_,
204 this -> basisCardinality_,
214 const ArrayScalar &inputPoints ,
215 const EOperator operatorType )
const
217#ifdef HAVE_INTREPID_DEBUG
218 getValues_HGRAD_Args<Scalar, ArrayScalar>(outputValues,
221 this -> getBaseCellTopology(),
222 this -> getCardinality() );
225 FieldContainer<Scalar> xInputPoints(inputPoints.dimension(0),1);
226 FieldContainer<Scalar> yInputPoints(inputPoints.dimension(0),1);
228 const Basis<Scalar,ArrayScalar> &xBasis_ = *this->bases_[0][0];
229 const Basis<Scalar,ArrayScalar> &yBasis_ = *this->bases_[0][1];
231 for (
int i=0;i<inputPoints.dimension(0);i++) {
232 xInputPoints(i,0) = inputPoints(i,0);
233 yInputPoints(i,0) = inputPoints(i,1);
236 switch (operatorType) {
239 FieldContainer<Scalar> xBasisValues(xBasis_.getCardinality(),xInputPoints.dimension(0));
240 FieldContainer<Scalar> yBasisValues(yBasis_.getCardinality(),yInputPoints.dimension(0));
242 xBasis_.getValues(xBasisValues,xInputPoints,OPERATOR_VALUE);
243 yBasis_.getValues(yBasisValues,yInputPoints,OPERATOR_VALUE);
246 for (
int j=0;j<yBasis_.getCardinality();j++) {
247 for (
int i=0;i<xBasis_.getCardinality();i++) {
248 for (
int k=0;k<inputPoints.dimension(0);k++) {
249 outputValues(bfcur,k) = xBasisValues(i,k) * yBasisValues(j,k);
259 FieldContainer<Scalar> xBasisValues(xBasis_.getCardinality(),xInputPoints.dimension(0));
260 FieldContainer<Scalar> yBasisValues(yBasis_.getCardinality(),yInputPoints.dimension(0));
261 FieldContainer<Scalar> xBasisDerivs(xBasis_.getCardinality(),xInputPoints.dimension(0),1);
262 FieldContainer<Scalar> yBasisDerivs(yBasis_.getCardinality(),yInputPoints.dimension(0),1);
264 xBasis_.getValues(xBasisValues,xInputPoints,OPERATOR_VALUE);
265 yBasis_.getValues(yBasisValues,yInputPoints,OPERATOR_VALUE);
266 xBasis_.getValues(xBasisDerivs,xInputPoints,OPERATOR_D1);
267 yBasis_.getValues(yBasisDerivs,yInputPoints,OPERATOR_D1);
272 for (
int j=0;j<yBasis_.getCardinality();j++) {
273 for (
int i=0;i<xBasis_.getCardinality();i++) {
274 for (
int k=0;k<inputPoints.dimension(0);k++) {
275 outputValues(bfcur,k,0) = xBasisDerivs(i,k,0) * yBasisValues(j,k);
276 outputValues(bfcur,k,1) = xBasisValues(i,k) * yBasisDerivs(j,k,0);
293 FieldContainer<Scalar> xBasisValues(xBasis_.getCardinality(),xInputPoints.dimension(0));
294 FieldContainer<Scalar> yBasisValues(yBasis_.getCardinality(),yInputPoints.dimension(0));
296 Teuchos::Array<int> partialMult;
300 if (partialMult[0] == 0) {
301 xBasisValues.resize(xBasis_.getCardinality(),xInputPoints.dimension(0));
302 xBasis_.getValues( xBasisValues , xInputPoints, OPERATOR_VALUE );
305 xBasisValues.resize(xBasis_.getCardinality(),xInputPoints.dimension(0),1);
306 EOperator xop = (EOperator) ( (
int) OPERATOR_D1 + partialMult[0] - 1 );
307 xBasis_.getValues( xBasisValues , xInputPoints, xop );
308 xBasisValues.resize(xBasis_.getCardinality(),xInputPoints.dimension(0));
310 if (partialMult[1] == 0) {
311 yBasisValues.resize(yBasis_.getCardinality(),yInputPoints.dimension(0));
312 yBasis_.getValues( yBasisValues , yInputPoints, OPERATOR_VALUE );
315 yBasisValues.resize(yBasis_.getCardinality(),yInputPoints.dimension(0),1);
316 EOperator yop = (EOperator) ( (
int) OPERATOR_D1 + partialMult[1] - 1 );
317 yBasis_.getValues( yBasisValues , yInputPoints, yop );
318 yBasisValues.resize(yBasis_.getCardinality(),yInputPoints.dimension(0));
323 for (
int j=0;j<yBasis_.getCardinality();j++) {
324 for (
int i=0;i<xBasis_.getCardinality();i++) {
325 for (
int k=0;k<inputPoints.dimension(0);k++) {
326 outputValues(bfcur,k,d) = xBasisValues(i,k) * yBasisValues(j,k);
336 FieldContainer<Scalar> xBasisValues(xBasis_.getCardinality(),xInputPoints.dimension(0));
337 FieldContainer<Scalar> yBasisValues(yBasis_.getCardinality(),yInputPoints.dimension(0));
338 FieldContainer<Scalar> xBasisDerivs(xBasis_.getCardinality(),xInputPoints.dimension(0),1);
339 FieldContainer<Scalar> yBasisDerivs(yBasis_.getCardinality(),yInputPoints.dimension(0),1);
341 xBasis_.getValues(xBasisValues,xInputPoints,OPERATOR_VALUE);
342 yBasis_.getValues(yBasisValues,yInputPoints,OPERATOR_VALUE);
343 xBasis_.getValues(xBasisDerivs,xInputPoints,OPERATOR_D1);
344 yBasis_.getValues(yBasisDerivs,yInputPoints,OPERATOR_D1);
349 for (
int j=0;j<yBasis_.getCardinality();j++) {
350 for (
int i=0;i<xBasis_.getCardinality();i++) {
351 for (
int k=0;k<inputPoints.dimension(0);k++) {
352 outputValues(bfcur,k,0) = xBasisValues(i,k) * yBasisDerivs(j,k,0);
353 outputValues(bfcur,k,1) = -xBasisDerivs(i,k,0) * yBasisValues(j,k);
361 TEUCHOS_TEST_FOR_EXCEPTION(
true , std::invalid_argument,
362 ">>> ERROR (Basis_HGRAD_QUAD_Cn_FEM): Operator type not implemented");