Commit 59160f27 authored by Christophe Geuzaine's avatar Christophe Geuzaine

_mm_malloc is not portable - maybe hxt should use hxtMalloc ? reverting to standard malloc for now

parent 4f5e071b
Pipeline #1387 failed with stage
in 39 minutes and 32 seconds
......@@ -279,7 +279,7 @@ HXTStatus hxtLinearSystemLUCreate(HXTLinearSystemLU **pSystem, int nElements, in
}
free(nodeRowStart);
free(nodeRowEnd);
system->M = _mm_malloc(sizeof(double)*totalSize, PADDING*8);
system->M = malloc(sizeof(double)*totalSize);
system->rows = malloc(sizeof(double*)*system->n);
for (int i = 0; i < totalSize; ++i)
system->M[i] = 0;
......@@ -291,7 +291,7 @@ HXTStatus hxtLinearSystemLUCreate(HXTLinearSystemLU **pSystem, int nElements, in
totalSize += system->rowEnd[i]-system->rowStart[i]+(paddedStart-start);
system->rows[i] = system->M + paddedStart;
}
system->x = _mm_malloc(sizeof(double)*system->n, PADDING*8);
system->x = malloc(sizeof(double)*system->n);
return HXT_STATUS_OK;
}
......@@ -384,7 +384,7 @@ HXTStatus hxtLinearSystemLUAddMatrixEntry(HXTLinearSystemLU *system, int node0,
HXT_ERROR_MSG(HXT_STATUS_FAILED, "node %i or %i not in the domain", node0, node1);
int row0 = system->nodeMap[node0]*system->nFields + field0;
int col1 = system->nodeMap[node1]*system->nFields + field1;
system->rows[row0][col1] += v;
return HXT_STATUS_OK;
}
......@@ -403,7 +403,7 @@ HXTStatus hxtLinearSystemLUSolve(HXTLinearSystemLU *system, double *rhs, double
LUPDecompose(system);
system->flaglu=1;
}
LUPSolve(system, rhs);
for (int i = 0; i < system->nNodes; ++i){
int ii = system->nodeMap[i];
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment