ceed-ref-basis.c (dc16b3712a706ab37fa09198d6046649aa573ff7) ceed-ref-basis.c (fd364f385fcf8d7fc7b4a70eff195ce997d4a2a4)
1// Copyright (c) 2017-2018, Lawrence Livermore National Security, LLC.
2// Produced at the Lawrence Livermore National Laboratory. LLNL-CODE-734707.
3// All Rights reserved. See files LICENSE and NOTICE for details.
4//
5// This file is part of CEED, a collection of benchmarks, miniapps, software
6// libraries and APIs for efficient high-order finite element and spectral
7// element discretizations for exascale applications. For more information and
8// source code availability see http://github.com/ceed.

--- 38 unchanged lines hidden (view full) ---

47
48 // Clear v if operating in transpose
49 if (tmode == CEED_TRANSPOSE) {
50 const CeedInt vsize = nelem*ncomp*nnodes;
51 for (CeedInt i = 0; i < vsize; i++)
52 v[i] = (CeedScalar) 0.0;
53 }
54 bool tensorbasis;
1// Copyright (c) 2017-2018, Lawrence Livermore National Security, LLC.
2// Produced at the Lawrence Livermore National Laboratory. LLNL-CODE-734707.
3// All Rights reserved. See files LICENSE and NOTICE for details.
4//
5// This file is part of CEED, a collection of benchmarks, miniapps, software
6// libraries and APIs for efficient high-order finite element and spectral
7// element discretizations for exascale applications. For more information and
8// source code availability see http://github.com/ceed.

--- 38 unchanged lines hidden (view full) ---

47
48 // Clear v if operating in transpose
49 if (tmode == CEED_TRANSPOSE) {
50 const CeedInt vsize = nelem*ncomp*nnodes;
51 for (CeedInt i = 0; i < vsize; i++)
52 v[i] = (CeedScalar) 0.0;
53 }
54 bool tensorbasis;
55 ierr = CeedBasisGetTensorStatus(basis, &tensorbasis); CeedChk(ierr);
55 ierr = CeedBasisIsTensor(basis, &tensorbasis); CeedChk(ierr);
56 // Tensor basis
57 if (tensorbasis) {
58 CeedInt P1d, Q1d;
59 ierr = CeedBasisGetNumNodes1D(basis, &P1d); CeedChk(ierr);
60 ierr = CeedBasisGetNumQuadraturePoints1D(basis, &Q1d); CeedChk(ierr);
61 switch (emode) {
62 // Interpolate to/from quadrature points
63 case CEED_EVAL_INTERP: {

--- 341 unchanged lines hidden --- 		56 // Tensor basis
57 if (tensorbasis) {
58 CeedInt P1d, Q1d;
59 ierr = CeedBasisGetNumNodes1D(basis, &P1d); CeedChk(ierr);
60 ierr = CeedBasisGetNumQuadraturePoints1D(basis, &Q1d); CeedChk(ierr);
61 switch (emode) {
62 // Interpolate to/from quadrature points
63 case CEED_EVAL_INTERP: {

--- 341 unchanged lines hidden ---