Skip to content
Snippets Groups Projects

Version 1.2

Merged Version 1.2
adri into master
Merged Adrien Crovato requested to merge adri into master
Viewing commit e476b7a3
Prev Next
Show latest version
2 files
+ 2
− 7
Compare changes
  • Side-by-side
  • Inline
Files
2
sdpm/src/sdpmSolver.cpp
+ 1
− 5
@@ -63,7 +63,6 @@ Solver::Solver(Problem &pbl) : _pbl(pbl), _cl(0), _cd(0), _cs(0), _cm(0)
_B.resize(nf, sdpmMatrixXcd::Zero(_nP, _nP));
// Set up variables
_phi.resize(ne, 0.);
_tau.resize(ne, 0.);
_mu.resize(ne, 0.);
_u.resize(ne, sdpmVector3d::Zero());
@@ -326,7 +325,6 @@ void Solver::save(GmshExport const &writer, std::string const &suffix)
// set up results
Results results;
// steady
results.scalarsAtElems["phi"] = &_phi;
results.scalarsAtElems["mu"] = &_mu;
results.scalarsAtElems["tau"] = &_tau;
results.vectorsAtElems["u"] = &_u;
@@ -407,8 +405,7 @@ void Solver::save(GmshExport const &writer, std::string const &suffix)
*/
void Solver::post(sdpmVectorXd const &etau, sdpmVectorXd const &emu)
{
// Store source and doublet strength and compute perturbation potential on elements
sdpmVectorXd ephi = -_A0 * emu - _B0 * etau;
// Store source and doublet strength on elements
std::vector<Body *> bodies = _pbl.getBodies();
for (auto body : bodies)
{
@@ -418,7 +415,6 @@ void Solver::post(sdpmVectorXd const &etau, sdpmVectorXd const &emu)
int i = _rows[e];
_tau[id] = etau(i);
_mu[id] = emu(i);
_phi[id] = ephi(i);
}
}