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

Fixing the heat equation solver to work with single precision.

parent 00c903f5

src/TNL/Functions/MeshFunction_impl.h

+1 −1
Original line number Diff line number Diff line
@@ -207,7 +207,7 @@ bind( const MeshPointer& meshPointer, 
{

   TNL_ASSERT_GE( data->getSize(), offset + meshPointer->template getEntitiesCount< typename MeshType::template EntityType< MeshEntityDimension > >(),

                  "The input vector is not large enough for binding to the mesh function." );



   static_assert( std::is_same< typename Vector::RealType, RealType >::value, "Cannot bind Vector with different Real type." );

   setupSynchronizer(meshPointer->getDistributedMesh());

   this->meshPointer=meshPointer;

   this->data.bind( *data, offset, getMesh().template getEntitiesCount< typename Mesh::template EntityType< MeshEntityDimension > >() );

src/TNL/Problems/HeatEquationProblem.h

+4 −4
Original line number Diff line number Diff line
@@ -38,15 +38,15 @@ template< typename Mesh, 
                                                              typename BoundaryCondition::RealType > >

class HeatEquationProblem : public PDEProblem< Mesh,

                                               Communicator,

                                               typename DifferentialOperator::RealType,

                                               typename Mesh::RealType,

                                               typename Mesh::DeviceType,

                                               typename DifferentialOperator::IndexType  >

                                               typename Mesh::IndexType  >

{

   public:



      typedef typename DifferentialOperator::RealType RealType;

      typedef typename Mesh::RealType RealType;

      typedef typename Mesh::DeviceType DeviceType;

      typedef typename DifferentialOperator::IndexType IndexType;

      typedef typename Mesh::IndexType IndexType;

      typedef Functions::MeshFunction< Mesh > MeshFunctionType;

      typedef Pointers::SharedPointer< MeshFunctionType, DeviceType > MeshFunctionPointer;

      typedef PDEProblem< Mesh, Communicator, RealType, DeviceType, IndexType > BaseType;