diff --git a/sdpm/src/sdpmSolver.cpp b/sdpm/src/sdpmSolver.cpp
index 56679fabd739b6f1ce63fbd8098a5a895fe84f57..1ee540b4c8a336dacb2bc75eaf6652cae50a5bf0 100644
--- a/sdpm/src/sdpmSolver.cpp
+++ b/sdpm/src/sdpmSolver.cpp
@@ -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);
}
}
diff --git a/sdpm/src/sdpmSolver.h b/sdpm/src/sdpmSolver.h
index a92ca0c2eae911cbbbbf72054f6b008acc54a37f..313f019818401410c44ba5d07f8f623c99a8c277 100644
--- a/sdpm/src/sdpmSolver.h
+++ b/sdpm/src/sdpmSolver.h
@@ -46,9 +46,8 @@ protected:
private:
// Solution
- std::vector<sdpmDouble> _phi; ///< perturbation potential
std::vector<sdpmDouble> _tau; ///< source strength
- std::vector<sdpmDouble> _mu; ///< doublet strength
+ std::vector<sdpmDouble> _mu; ///< doublet strength (perturbation potential)
std::vector<sdpmVector3d> _u; ///< steady velocity
std::vector<sdpmDouble> _cp; ///< steady pressure coefficient
std::vector<std::vector<std::vector<sdpmVector3d>>> _u1Re; ///< real part of unsteady first harmonic velocity