Skip to content
Snippets Groups Projects
Select Git revision
  • 6d2114aac14613bceb69d488811c078870daf192
  • master default protected
  • overlaps_tags_and_distributed_export
  • 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
  • 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

HierarchicalBasis.h

Blame
  • FuncHeaviside.h 1.07 KiB
    //
    // Description : Heaviside function based on level set discontinuity description
    //
    //
    // Author:  <Boris Sedji>,  01/2010
    //
    // Copyright: See COPYING file that comes with this distribution
    //
    //
    
    #ifndef _FUNCHEAVISIDE_H_
    #define _FUNCHEAVISIDE_H_
    
    #include "simpleFunction.h"
    #include "DILevelset.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