diff --git a/FunctionSpace/TetEdgeBasis.cpp b/FunctionSpace/TetEdgeBasis.cpp index ddcfefd9b9a2f080c8e05fca9cf94bdf7409d2e6..f90eeeaba0a71a160db89a7d320a89829c57b2eb 100644 --- a/FunctionSpace/TetEdgeBasis.cpp +++ b/FunctionSpace/TetEdgeBasis.cpp @@ -172,22 +172,16 @@ TetEdgeBasis::TetEdgeBasis(size_t order){ // Type 1 - basis[s][i] = - new vector<Polynomial>((u * v).gradient()); - + basis[s][i] = new vector<Polynomial>((u * v).gradient()); i++; // Type 2 - basis[s][i] = - new vector<Polynomial>(subGradUV); - + basis[s][i] = new vector<Polynomial>(subGradUV); i++; // Type 3 if(l1 == 1){ - basis[s][i] = - new vector<Polynomial>(subGradL1L2V); - + basis[s][i] = new vector<Polynomial>(subGradL1L2V); i++; } } diff --git a/FunctionSpace/TriEdgeBasis.cpp b/FunctionSpace/TriEdgeBasis.cpp index 516c055bc0719a9fc5b5c8d3858f0497ff4e9882..6066dcd7e14398f708d8a91abe4024d0ea1f83d1 100644 --- a/FunctionSpace/TriEdgeBasis.cpp +++ b/FunctionSpace/TriEdgeBasis.cpp @@ -51,6 +51,9 @@ TriEdgeBasis::TriEdgeBasis(size_t order){ Polynomial(Polynomial(1, 0, 1, 0)) }; + // One // + Polynomial one(1, 0, 0, 0); + // Basis // basis = new vector<Polynomial>**[nOrientation]; @@ -118,7 +121,7 @@ TriEdgeBasis::TriEdgeBasis(size_t order){ lagrange[faceIdx[s][0][1]]); Polynomial v = lagrange[faceIdx[s][0][2]] * - legendre[l2].compose(lagrange[faceIdx[s][0][2]] * 2); + legendre[l2].compose(lagrange[faceIdx[s][0][2]] * 2 - one); // Preliminary Type 2 vector<Polynomial> gradU = u.gradient(); @@ -164,22 +167,16 @@ TriEdgeBasis::TriEdgeBasis(size_t order){ // Type 1 - basis[s][i] = - new vector<Polynomial>((u * v).gradient()); - + basis[s][i] = new vector<Polynomial>((u * v).gradient()); i++; // Type 2 - basis[s][i] = - new vector<Polynomial>(subGradUV); - + basis[s][i] = new vector<Polynomial>(subGradUV); i++; // Type 3 if(l1 == 1){ - basis[s][i] = - new vector<Polynomial>(subGradL1L2V); - + basis[s][i] = new vector<Polynomial>(subGradL1L2V); i++; } }