Remove dedicated variable for the perturbation potential

a311999e · Adrien Crovato · cfada2f9 · a311999e · a311999e
Commit a311999e authored 1 year ago by Adrien Crovato
--- 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());
@@ -327,7 +326,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;
@@ -408,8 +406,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)
    {
@@ -419,7 +416,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);
        }
    }


--- a/sdpm/src/sdpmSolver.h
+++ b/sdpm/src/sdpmSolver.h
@@ -43,9 +43,8 @@ class SDPM_API Solver : public sdpmObject
    std::vector<sdpmMatrixXcd> _A;  ///< unsteady body-body doublet matrix
    std::vector<sdpmMatrixXcd> _B;  ///< unsteady body-body source matrix
    // 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