DGFunction.cpp 3.12 KB
Newer Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
/*
 * Copyright (c) 2017-2019 Boerge Struempfel, Daniel Drzisga, Dominik Thoennes, Nils Kohl.
 *
 * This file is part of HyTeG
 * (see https://i10git.cs.fau.de/hyteg/hyteg).
 *
 * This program is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program. If not, see <http://www.gnu.org/licenses/>.
 */

#include "hyteg/dgfunctionspace/DGFunction.hpp"

#include "hyteg/dgfunctionspace/DGFunctionMacroEdge.hpp"
#include "hyteg/dgfunctionspace/DGFunctionMacroFace.hpp"
#include "hyteg/dgfunctionspace/DGFunctionMacroVertex.hpp"

namespace hyteg {

template < typename ValueType >
void DGFunction< ValueType >::projectP1( P1Function< ValueType >& src, uint_t level, DoFType flag, UpdateType updateType )
{
   this->startTiming( "projectP1" );

   src.template startCommunication< Edge, Vertex >( level );
   src.template startCommunication< Face, Edge >( level );
   src.template endCommunication< Edge, Vertex >( level );

   for ( auto& it : this->getStorage()->getVertices() )
   {
      Vertex& vertex = *it.second;

      const DoFType vertexBC = this->getBoundaryCondition().getBoundaryType( vertex.getMeshBoundaryFlag() );
      if ( testFlag( vertexBC, flag ) )
      {
         dgfunction::macrovertex::projectP1< ValueType >(
             level, vertex, this->getStorage(), src.getVertexDataID(), this->getVertexDataID(), updateType );
      }
   }

   this->template startCommunication< Vertex, Edge >( level );

   src.template endCommunication< Face, Edge >( level );

   for ( auto& it : this->getStorage()->getEdges() )
   {
      Edge& edge = *it.second;

      const DoFType edgeBC = this->getBoundaryCondition().getBoundaryType( edge.getMeshBoundaryFlag() );

      if ( testFlag( edgeBC, flag ) )
      {
         dgfunction::macroedge::projectP1< ValueType >(
             level, edge, this->getStorage(), src.getEdgeDataID(), this->getEdgeDataID(), updateType );
      }
   }

   this->template endCommunication< Vertex, Edge >( level );

   this->template startCommunication< Edge, Face >( level );

   for ( auto& it : this->getStorage()->getFaces() )
   {
      Face& face = *it.second;

      const DoFType faceBC = this->getBoundaryCondition().getBoundaryType( face.getMeshBoundaryFlag() );
      if ( testFlag( faceBC, flag ) )
      {
         dgfunction::macroface::projectP1< ValueType >(
             level, face, this->getStorage(), src.getFaceDataID(), this->getFaceDataID(), updateType );
      }
   }

   this->template endCommunication< Edge, Face >( level );

   this->stopTiming( "projectP1" );
}

template class DGFunction< real_t >;
template class DGFunction< int32_t >;
template class DGFunction< int64_t >;

} // namespace hyteg