dbpp_LegacyNumericalAlgo.h

Go to the documentation of this file.

#pragma once
 
//C++ include
#include <vector>
// SfxBase19 include
#include "include/Sfx_DefineTypes.h"
 
namespace dbpp 
{
    void TraitementTermeP(std::vector<float64>& PF2, std::vector<float64>& P2,
        const std::vector<float64>& U1, const int NbSections, const float64 B = 1);
    void CalculS0(std::vector<float64>& S0, const std::vector<float64>& Z, const float64 dx, const int NbSections);
    void TraitementTermeSource(std::vector<float64>& S, const std::vector<float64>& Q,
        const std::vector<float64>& A, const std::vector<float64>& Z,
        const std::vector<float64>& n, const float64 dx, const int NbSections, const float64 B = 1);
    void TraitementTermeSource2(std::vector<float64>& S, const std::vector<float64>& Q,
        const std::vector<float64>& A, const std::vector<float64>& H,
        const std::vector<float64>& n, const float64 dx, const int NbSections = 101, const float64 B = 1.);
 
    // use Rvalue reference to modify some variable by adding b.c. in the computation phase
    // in our programming we impose b.c. by adding a ghost node, but we don't want to modify
    // the original which is a scalar field (bind temporary to rvalue reference)
    // Usage: TreatmentSource2(w_s, aQ.values().asStdVector(), aA.values().asStdVector(), ...)
    //
    void TraitementTermeSource2(std::vector<double>& aS, // return value
        std::vector<double>&& aQ, std::vector<double>&& aA, std::vector<double>&& aH, //rvalue reference
        const std::vector<double>& n, double aDx /*= 10.*/, unsigned aDim = 100, const float64 B = 1.);
 
    // deprecated
    //void CalculFFHLL(std::vector<float64>& FF1, std::vector<float64>& FF2, const std::vector<float64>& U1,
    //  const std::vector<float64>& U2, std::vector<float64>& dU1, std::vector<float64>& dU2,
    //  const int NbSections, const float64 B = 1.);
 
    void CalculFFHLL(std::vector<float64>& FF1, std::vector<float64>& FF2, const std::vector<float64>& U1,
        const std::vector<float64>& U2, const float64 B = 1.);
} // End of namespace