Commit 4ce5b65b authored by Thomas Witkowski's avatar Thomas Witkowski

* removed axpy from ProblemVec

parent 8f4f2d4d
......@@ -95,4 +95,6 @@
#include "Triangle.h"
#include "ValueWriter.h"
#include "demangle.h"
#include "parareal/ProblemBase.h"
#include "parareal/AdaptParaReal.h"
#endif
......@@ -1014,32 +1014,5 @@ namespace AMDiS {
solution_->deserialize(in);
}
void axpy(double a,
const Matrix<DOFMatrix*> &x,
Matrix<DOFMatrix*> &y)
{
FUNCNAME("ProblemVec::axpy()");
int numRows = x.getNumRows(), numCols = x.getNumCols();
int i, j;
const DOFMatrix *xMatrix;
DOFMatrix *yMatrix;
for(i = 0; i < numRows; i++) {
for(j = 0; j < numCols; j++) {
xMatrix = x[i][j];
yMatrix = y[i][j];
if(!xMatrix) continue;
if(!yMatrix) {
yMatrix = NEW DOFMatrix(xMatrix->getRowFESpace(),
xMatrix->getColFESpace(), "");
y[i][j] = yMatrix;
}
yMatrix->axpy(a, *xMatrix, *yMatrix);
}
}
}
}
......@@ -641,12 +641,6 @@ namespace AMDiS {
*/
static ProblemVec *traversePtr_;
};
void axpy(double a,
const Matrix<DOFMatrix*> &x,
Matrix<DOFMatrix*> &y);
}
#endif
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