Skip to content
Snippets Groups Projects
Commit b89c9209 authored by Nicolas Marsic's avatar Nicolas Marsic
Browse files

Fix illegal corruption in Basis::getDerivative for 0 and 1 forms

parent 57e7cb5e
No related branches found
No related tags found
No related merge requests found
......@@ -82,21 +82,10 @@ getFunctions(fullMatrix<double>& retValues,
void BasisHierarchical0Form::getDerivative(fullMatrix<double>& retValues,
const MElement& element,
double u, double v, double w) const{
const size_t nOrientation = ReferenceSpaceManager::getNOrientation(getType());
// Get Grad //
if(!hasGrad)
getGrad();
// Delete if older //
if(preEvaluatedGrad){
for(size_t i = 0; i < nOrientation; i++)
delete preEvaluatedGradFunction[i];
delete[] preEvaluatedGradFunction;
}
// Define Orientation //
const size_t orientation = ReferenceSpaceManager::getOrientation(element);
......
......@@ -92,21 +92,10 @@ getFunctions(fullMatrix<double>& retValues,
void BasisHierarchical1Form::getDerivative(fullMatrix<double>& retValues,
const MElement& element,
double u, double v, double w) const{
const size_t nOrientation = ReferenceSpaceManager::getNOrientation(getType());
// Build Curl //
if(!hasCurl)
getCurl();
// Delete if older //
if(preEvaluatedCurl){
for(size_t i = 0; i < nOrientation; i++)
delete preEvaluatedCurlFunction[i];
delete[] preEvaluatedCurlFunction;
}
// Define Orientation //
const size_t orientation = ReferenceSpaceManager::getOrientation(element);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment