Skip to content
Snippets Groups Projects
Select Git revision
  • 3f147e49852e19b79b873f3f88edaff2314ae196
  • master default protected
  • patches-4.14
  • steplayer
  • bl
  • pluginMeshQuality
  • fixBugsAmaury
  • hierarchical-basis
  • alphashapes
  • relaying
  • new_export_boris
  • oras_vs_osm
  • reassign_partitions
  • distributed_fwi
  • rename-classes
  • fix/fortran-api-example-t4
  • robust_partitions
  • reducing_files
  • fix_overlaps
  • 3115-issue-fix
  • 3023-Fillet2D-Update
  • gmsh_4_14_0
  • gmsh_4_13_1
  • gmsh_4_13_0
  • gmsh_4_12_2
  • gmsh_4_12_1
  • gmsh_4_12_0
  • gmsh_4_11_1
  • gmsh_4_11_0
  • gmsh_4_10_5
  • gmsh_4_10_4
  • gmsh_4_10_3
  • gmsh_4_10_2
  • gmsh_4_10_1
  • gmsh_4_10_0
  • gmsh_4_9_5
  • gmsh_4_9_4
  • gmsh_4_9_3
  • gmsh_4_9_2
  • gmsh_4_9_1
  • gmsh_4_9_0
41 results

FuncHeaviside.h

Blame
  • FuncHeaviside.h 1.12 KiB
    // Gmsh - Copyright (C) 1997-2019 C. Geuzaine, J.-F. Remacle
    //
    // See the LICENSE.txt file for license information. Please report all
    // issues on https://gitlab.onelab.info/gmsh/gmsh/issues.
    //
    // Contributor(s):
    //   Boris Sedji
    //
    
    #ifndef FUNCHEAVISIDE_H
    #define FUNCHEAVISIDE_H
    
    #include "simpleFunction.h"
    #include "gmshLevelset.h"
    
    class FuncHeaviside : public simpleFunctionOnElement<double> {
    private:
      gLevelset *_ls;
    
    public:
      FuncHeaviside(gLevelset *ls) : _ls(ls) {}
      virtual double operator()(double x, double y, double z) const
      {
        if(_ls->isInsideDomain(x, y, z))
          return 1;
        else
          return -1;
      }
      virtual double operator()(double x, double y, double z, MElement *e) const
      {
        if(_ls->isInsideDomain(x, y, z))
          return 1;
        else
          return -1;
      }
      virtual void gradient(double x, double y, double z, double &dfdx,
                            double &dfdy, double &dfdz) const
      {
        dfdx = dfdy = dfdz = 0.0;
      }
      virtual void gradient(double x, double y, double z, double &dfdx,
                            double &dfdy, double &dfdz, MElement *e) const
      {
        dfdx = dfdy = dfdz = 0.0;
      }
    };
    
    #endif