lsVelocityField.hpp

Go to the documentation of this file.

#pragma once
 
#include <array>
 
#include <vcVectorUtil.hpp>
 
namespace viennals {
 
using namespace viennacore;
 
template <class T> class VelocityField {
public:
  VelocityField() {}
 
  virtual T getScalarVelocity(const Vec3D<T> & /*coordinate*/, int /*material*/,
                              const Vec3D<T> & /*normalVector*/,
                              unsigned long /*pointId*/) {
    return 0;
  }
 
  virtual Vec3D<T> getVectorVelocity(const Vec3D<T> & /*coordinate*/,
                                     int /*material*/,
                                     const Vec3D<T> & /*normalVector*/,
                                     unsigned long /*pointId*/) {
    return {0, 0, 0};
  }
 
  virtual T getDissipationAlpha(int /*direction*/, int /*material*/,
                                const Vec3D<T> & /*centralDifferences*/) {
    return 0;
  }
 
  virtual ~VelocityField() {}
};
 
} // namespace viennals