Commit 91718fb9 authored by Thomas Witkowski's avatar Thomas Witkowski

* Bugfix for file writing of higher order elements

parent 6c978c46
......@@ -351,9 +351,9 @@ namespace AMDiS {
basisFcts_->getLocalIndices(elInfo->getElement(), localAdmin_, localDOFs_);
::std::vector<int> elemInterpPoints(nBasisFcts_);
std::vector<int> elemInterpPoints(0);
for (int i = mesh_->getGeo(VERTEX); i < nBasisFcts_; i++) {
elemInterpPoints[i] = (*interpPointInd_)[localDOFs_[i]];
elemInterpPoints.push_back((*interpPointInd_)[localDOFs_[i]]);
}
interpPoints_.push_back(elemInterpPoints);
......
......@@ -132,10 +132,15 @@ namespace AMDiS {
void VtkWriter::writeVertexValues(::std::ofstream &file, int componentNo)
{
DOFVector<int> *interpPointInd;
DOFVector<double> *values;
DOFVector< ::std::list<WorldVector<double> > > *dofCoords;
// DOFVector<int> *interpPointInd;
// DOFVector<double> *values;
// DOFVector< ::std::list<WorldVector<double> > > *dofCoords;
DOFVector<int> *interpPointInd = (*dc_)[componentNo]->getInterpPointInd();
DOFVector<double> *values = (*dc_)[componentNo]->getValues();
DOFVector< ::std::list<WorldVector<double> > > *dofCoords = (*dc_)[componentNo]->getDofCoords();
/*
#ifdef _OPENMP
#pragma omp critical
#endif
......@@ -144,7 +149,7 @@ namespace AMDiS {
values = (*dc_)[componentNo]->getValues();
dofCoords = (*dc_)[componentNo]->getDofCoords();
}
*/
DOFVector<int>::Iterator intPointIt(interpPointInd, USED_DOFS);
DOFVector<double>::Iterator valueIt(values, USED_DOFS);
DOFVector< ::std::list<WorldVector<double> > >::Iterator coordIt(dofCoords, USED_DOFS);
......
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