vtkAMRGaussianPulseSource.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
// SPDX-License-Identifier: BSD-3-Clause
#ifndef vtkAMRGaussianPulseSource_h
#define vtkAMRGaussianPulseSource_h
 
#include "vtkFiltersAMRModule.h" // For export macro
#include "vtkOverlappingAMRAlgorithm.h"
 
#include <cmath> // For std::exp
 
VTK_ABI_NAMESPACE_BEGIN
class vtkOverlappingAMR;
class vtkUniformGrid;
class vtkInformation;
class vtkInformationVector;
 
class VTKFILTERSAMR_EXPORT vtkAMRGaussianPulseSource : public vtkOverlappingAMRAlgorithm
{
public:
  static vtkAMRGaussianPulseSource* New();
  vtkTypeMacro(vtkAMRGaussianPulseSource, vtkOverlappingAMRAlgorithm);
  void PrintSelf(ostream& os, vtkIndent indent) override;
 
 
  vtkSetMacro(Dimension, int);
 
 
  vtkSetMacro(NumberOfLevels, int);
 
  void SetRefinementRatio(int r)
  {
    this->RefinmentRatio = r;
    this->Modified();
  }
 
 
  void SetRootSpacing(double h0)
  {
    this->RootSpacing[0] = this->RootSpacing[1] = this->RootSpacing[2] = h0;
    this->Modified();
  }
 
 
  vtkSetVector3Macro(PulseOrigin, double);
  vtkGetVector3Macro(PulseOrigin, double);
  void SetXPulseOrigin(double f)
  {
    this->PulseOrigin[0] = f;
    this->Modified();
  }
  void SetYPulseOrigin(double f)
  {
    this->PulseOrigin[1] = f;
    this->Modified();
  }
  void SetZPulseOrigin(double f)
  {
    this->PulseOrigin[2] = f;
    this->Modified();
  }
 
 
  vtkSetVector3Macro(PulseWidth, double);
  vtkGetVector3Macro(PulseWidth, double);
  void SetXPulseWidth(double f)
  {
    this->PulseWidth[0] = f;
    this->Modified();
  }
  void SetYPulseWidth(double f)
  {
    this->PulseWidth[1] = f;
    this->Modified();
  }
  void SetZPulseWidth(double f)
  {
    this->PulseWidth[2] = f;
    this->Modified();
  }
 
 
  vtkSetMacro(PulseAmplitude, double);
  vtkGetMacro(PulseAmplitude, double);
 
protected:
  vtkAMRGaussianPulseSource();
  ~vtkAMRGaussianPulseSource() override;
 
  int RequestData(vtkInformation* request, vtkInformationVector** inputVector,
    vtkInformationVector* outputVector) override;
 
 
  double ComputePulseAt(const double x, const double y, const double z)
  {
    double xyz[3];
    xyz[0] = x;
    xyz[1] = y;
    xyz[2] = z;
    return (this->ComputePulseAt(xyz));
  }
  double ComputePulseAt(double pt[3])
  {
    double pulse = 0.0;
    double r = 0.0;
    for (int i = 0; i < this->Dimension; ++i)
    {
      double d = pt[i] - this->PulseOrigin[i];
      double d2 = d * d;
      double L2 = this->PulseWidth[i] * this->PulseWidth[i];
      r += d2 / L2;
    }
    pulse = this->PulseAmplitude * std::exp(-r);
    return (pulse);
  }
 
  void ComputeCellCenter(vtkUniformGrid* grid, vtkIdType cellIdx, double centroid[3]);
 
  void GeneratePulseField(vtkUniformGrid* grid);
 
  vtkUniformGrid* GetGrid(double origin[3], double h[3], int ndim[3]);
 
  vtkUniformGrid* RefinePatch(vtkUniformGrid* parent, int patchExtent[6]);
 
 
  void Generate2DDataSet(vtkOverlappingAMR* amr);
  void Generate3DDataSet(vtkOverlappingAMR* amr);
 
  double RootSpacing[3];
  double PulseOrigin[3];
  double PulseWidth[3];
  double PulseAmplitude;
  int RefinmentRatio;
  int Dimension;
  int NumberOfLevels;
 
private:
  vtkAMRGaussianPulseSource(const vtkAMRGaussianPulseSource&) = delete;
  void operator=(const vtkAMRGaussianPulseSource&) = delete;
};
 
VTK_ABI_NAMESPACE_END
#endif /* vtkAMRGaussianPulseSource_h */