Commit d46c6604 authored by Tomáš Oberhuber's avatar Tomáš Oberhuber
Browse files

Fixing mean curvature test.

parent 4a647ea1

examples/CMakeLists.txt

+1 −1
Original line number Diff line number Diff line 
add_subdirectory( heat-equation )

add_subdirectory( navier-stokes )

add_subdirectory( mean-curvature-flow )
 
#add_subdirectory( mean-curvature-flow )

examples/mean-curvature-flow/tnl-mean-curvature-flow.h

+3 −4
Original line number Diff line number Diff line
@@ -25,10 +25,9 @@ 
#include <operators/tnlNeumannBoundaryConditions.h>

#include <functions/tnlConstantFunction.h>

#include <problems/tnlMeanCurvatureFlowProblem.h>

#include <operators/diffusion/tnlNonlinearDiffusion.h>

#include <operators/diffusion/tnlOneSidedNonlinearDiffusion.h>

#include <operators/operator-Q/tnlOneSideDiffOperatorQ.h>

#include <operators/operator-Q/tnlFiniteVolumeOperatorQ.h>

#include <operators/diffusion/nonlinear-diffusion-operators/tnlOneSideDiffNonlinearOperator.h>

#include <operators/diffusion/nonlinear-diffusion-operators/tnlFiniteVolumeNonlinearOperator.h>

#include <functions/tnlMeshFunction.h>
@@ -85,7 +84,7 @@ class meanCurvatureFlowSetter 
      if( numericalScheme == "fdm" )

      {

         typedef tnlOneSideDiffOperatorQ<MeshType, Real, Index > QOperator;

         typedef tnlOneSideDiffNonlinearOperator<MeshType, QOperator, Real, Index > NonlinearOperator;         

         typedef tnlOneSideNonlinearDiffusion<MeshType, QOperator, Real, Index > NonlinearOperator;         

         return setBoundaryConditions< NonlinearOperator, QOperator >( parameters );

      }

      if( numericalScheme == "fvm" )
@@ -101,7 +100,7 @@ class meanCurvatureFlowSetter 
             typename QOperator >

   static bool setBoundaryConditions( const tnlParameterContainer& parameters )

   {

      typedef tnlNonlinearDiffusion< MeshType, NonlinearOperator, Real, Index > ApproximateOperator;

      typedef tnlOneSidedNonlinearDiffusion< MeshType, NonlinearOperator, Real, Index > ApproximateOperator;

      typedef tnlConstantFunction< Dimensions, Real > RightHandSide;

      typedef tnlStaticVector< MeshType::meshDimensions, Real > Vertex;

src/functions/CMakeLists.txt

+2 −2
Original line number Diff line number Diff line
@@ -15,9 +15,9 @@ SET( headers tnlConstantFunction.h 
             tnlMeshFunctionEvaluator_impl.h

             tnlMeshFunctionGnuplotWriter.h

             tnlMeshFunctionGnuplotWriter_impl.h

             tnlMeshFunctionNormGetter.h

             tnlMeshFunctionVTKWriter.h

             tnlMeshFunctionVTKWriter_impl.h             

             tnlMeshFunctionNormGetter.h

             tnlOperatorFunction.h

             tnlSinBumpsFunction.h

             tnlSinBumpsFunction_impl.h

src/functions/tnlConstantFunction.h

+1 −1
Original line number Diff line number Diff line
@@ -60,7 +60,7 @@ class tnlConstantFunction : public tnlDomain< dimensions, NonspaceDomain > 
      RealType operator()( const VertexType& v,

                           const Real& time = 0.0 ) const

      {

         return getValue< 0, 0, 0 >( v, time );

         return constant;

      }

      

       __cuda_callable__ inline

src/functions/tnlMeshFunctionEvaluator.h

+1 −1
Original line number Diff line number Diff line
@@ -49,7 +49,7 @@ class tnlMeshFunctionEvaluator : public tnlDomain< OutMeshFunction::getEntitiesD 
      typedef tnlMeshFunctionEvaluatorTraverserUserData< OutMeshFunction, InFunction, RealType > TraverserUserData;



      

      const static int meshEntityDimensions = OutMeshFunction::entityDimensions;

      const static int meshEntityDimensions = OutMeshFunction::getEntitiesDimensions();

      

      static_assert( MeshType::meshDimensions == InFunction::getDimensions(), 

         "Input function and the mesh of the mesh function have both different number of dimensions." );