Commit 3e65261d authored by Bernd Flemisch's avatar Bernd Flemisch
Browse files

removed argument u from PDELabBoxAssembler::assemble

git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@4504 2fb0f335-1f38-0410-981e-8018bf24f1b0
parent 4a6a83dd
......@@ -165,7 +165,7 @@ public:
bool inJacobianAssemble() const
{ return inJacobianAssemble_; }
void assemble(SolutionVector &u)
void assemble()
{
// assemble the global jacobian matrix
if (!reuseMatrix_) {
......@@ -182,14 +182,14 @@ public:
// we actually need to reassemle!
resetMatrix_();
inJacobianAssemble_ = true;
gridOperatorSpace_->jacobian(u, *matrix_);
gridOperatorSpace_->jacobian(model_().curSol(), *matrix_);
inJacobianAssemble_ = false;
}
reuseMatrix_ = false;
// calculate the global residual
residual_ = 0;
gridOperatorSpace_->residual(u, residual_);
gridOperatorSpace_->residual(model_().curSol(), residual_);
typedef typename Matrix::block_type BlockType;
// set the entries for the ghost nodes
......@@ -203,7 +203,7 @@ public:
(*matrix_)[globI] = 0;
(*matrix_)[globI][globI] = Id;
residual_[globI] = 0;
u[globI] = 0;
model_().curSol()[globI] = 0;
}
}
......
......@@ -137,7 +137,7 @@ protected:
std::cout.flush();
}
// linearize the problem at the current solution
jacobianAsm.assemble(u);
jacobianAsm.assemble();
// solve the resultuing linear equation system
if (ctl.verbose()) {
......
Supports Markdown
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