fdm_mg_linear_system3.h

Go to the documentation of this file.

// Copyright (c) 2018 Doyub Kim
//
// I am making my contributions/submissions to this project solely in my
// personal capacity and am not conveying any rights to any intellectual
// property of any third parties.
 
#ifndef INCLUDE_JET_FDM_MG_LINEAR_SYSTEM3_H_
#define INCLUDE_JET_FDM_MG_LINEAR_SYSTEM3_H_
 
#include <jet/face_centered_grid3.h>
#include <jet/fdm_linear_system3.h>
#include <jet/mg.h>
 
namespace jet {
 
typedef MgMatrix<FdmBlas3> FdmMgMatrix3;
 
typedef MgVector<FdmBlas3> FdmMgVector3;
 
struct FdmMgLinearSystem3 {
    FdmMgMatrix3 A;
 
    FdmMgVector3 x;
 
    FdmMgVector3 b;
 
    void clear();
 
    size_t numberOfLevels() const;
 
    void resizeWithCoarsest(const Size3 &coarsestResolution,
                            size_t numberOfLevels);
 
    void resizeWithFinest(const Size3 &finestResolution,
                          size_t maxNumberOfLevels);
};
 
class FdmMgUtils3 {
 public:
    static void restrict(const FdmVector3 &finer, FdmVector3 *coarser);
 
    static void correct(const FdmVector3 &coarser, FdmVector3 *finer);
 
    template <typename T>
    static void resizeArrayWithCoarsest(const Size3 &coarsestResolution,
                                        size_t numberOfLevels,
                                        std::vector<Array3<T>> *levels);
 
    template <typename T>
    static void resizeArrayWithFinest(const Size3 &finestResolution,
                                      size_t maxNumberOfLevels,
                                      std::vector<Array3<T>> *levels);
};
 
}  // namespace jet
 
#include "detail/fdm_mg_linear_system3-inl.h"
 
#endif  // INCLUDE_JET_FDM_MG_LINEAR_SYSTEM3_H_