lsVelocityField.hpp

Go to the documentation of this file.

#pragma once
 
#include <vcVectorType.hpp>
 
namespace viennals {
 
using namespace viennacore;

template <class T> class VelocityField {
public:
  VelocityField() = default;

  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 Vec3D<T>{0, 0, 0};
  }

  virtual T getDissipationAlpha(int /*direction*/, int /*material*/,
                                const Vec3D<T> & /*centralDifferences*/) {
    return 0;
  }
 
  virtual ~VelocityField() = default;
};
 
} // namespace viennals