dufort_frankel.hpp

Go to the documentation of this file.

#pragma once
/**
 * @file dufort_frankel.hpp
 * @brief DuFort–Frankel explicit method (3-level scheme).
 *
 * Update:
 *   T_i^{n+1} = [ (1 - 2r) T_i^{n-1} + 2r (T_{i+1}^n + T_{i-1}^n) ] / (1 + 2r),
 *   with r = D dt / dx^2.
 */
#include "method.hpp"
 
class DuFortFrankel final : public Method
{
  public:
    std::string name() const override
    {
        return "DuFort-Frankel";
    }
    void step(const Grid& g,
              double D,
              double dt,
              const std::vector<double>& Tprev,
              const std::vector<double>& Tcurr,
              std::vector<double>& Tnext) const override;
    bool uses_previous_step() const noexcept override
    {
        return true;
    }
};