│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
106 EntityShiftTable() =
delete;
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
109 static constexpr int nextOffset(
int& offset,
int codim)
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
116 offset += subEnt<dim>(dim, codim-1);
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
120 template<std::size_t... I>
│ │ │ │ +
121 static constexpr std::array<int,
sizeof...(I)> computeOffsets(std::index_sequence<I...>)
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
124 return {{ (nextOffset(offset, I))... }};
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
128 static constexpr unsigned char nextValue(
int& codim,
int& i)
│ │ │ │ +
│ │ │ │ +
130 const auto result = F::evaluate(i, codim);
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
133 if (i >= subEnt<dim>(dim, codim)) {
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
141 template<std::size_t... I>
│ │ │ │ +
142 static constexpr std::array<
unsigned char,
sizeof...(I)> computeValues(std::index_sequence<I...>)
│ │ │ │ +
│ │ │ │ +
144 int codim = 0, i = 0;
│ │ │ │ +
145 return {{ ((void)I, nextValue(codim, i))... }};
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
148 static constexpr std::array<int,dim+1> _offsets = computeOffsets(std::make_index_sequence<dim+1>{});
│ │ │ │ +
149 static constexpr std::array<
unsigned char,Dune::power(3,dim)> _values = computeValues(std::make_index_sequence<Dune::power(3,dim)>{});
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
155 struct calculate_entity_shift
│ │ │ │ +
│ │ │ │ +
157 static constexpr unsigned long long evaluate(
int index,
int cc)
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
160 for (
int d = dim; d>0; d--)
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
164 if (index < subEnt<dim>(d-1,cc))
│ │ │ │ +
165 result |= 1ull << (d-1);
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
168 index = (index - subEnt<dim>(d-1, cc)) % subEnt<dim>(d-1,cc-1);
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
185 std::bitset<dim> entityShift(
int index,
int cc)
│ │ │ │ +
│ │ │ │ +
187 return EntityShiftTable<calculate_entity_shift<dim>,dim>::evaluate(index,cc);
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
192 struct calculate_entity_move
│ │ │ │ +
│ │ │ │ +
194 static constexpr unsigned long long evaluate(
int index,
int cc)
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
197 for (
int d = dim; d>0; d--)
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
202 result &= ~(1ull << (d-1));
│ │ │ │ +
203 result |= index & (1ull << (d-1));
│ │ │ │ +
│ │ │ │ +
205 index &= ~(1<<(d-1));
│ │ │ │ +
│ │ │ │ +
207 if (index >= subEnt<dim>(d-1,cc))
│ │ │ │ +
│ │ │ │ +
209 if ((index - subEnt<dim>(d-1,cc)) / subEnt<dim>(d-1,cc-1) == 1)
│ │ │ │ +
│ │ │ │ +
211 result |= 1ull << (d-1);
│ │ │ │ +
│ │ │ │ +
213 index = (index - subEnt<dim>(d-1, cc)) % subEnt<dim>(d-1,cc-1);
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
230 std::bitset<dim> entityMove(
int index,
int cc)
│ │ │ │ +
│ │ │ │ +
232 return EntityShiftTable<calculate_entity_move<dim>,dim>::evaluate(index,cc);
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
239 template<
int codim,
int dim,
class Gr
idImp>
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
244 template<
int, PartitionIteratorType,
typename>
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
248 typedef typename GridImp::ctype
ctype;
│ │ │ │ +
│ │ │ │ +
250 typedef typename GridImp::template Codim<codim>::Geometry
Geometry;
│ │ │ │ +
251 typedef typename GridImp::Traits::template Codim<codim>::GeometryImpl
GeometryImpl;
│ │ │ │ +
│ │ │ │ +
253 typedef typename GridImp::template Codim<codim>::EntitySeed
EntitySeed;
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
279 constexpr GeometryType
type ()
const
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
307 for (
int j=0; j<d; ++j)
│ │ │ │ -
308 dist += std::abs(
delta[j]);
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
290 return Dune::Yasp::subEnt<dim>(dim-codim,cc-codim);
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
296 if (
_g->interior[codim].inside(
_it.coord(),
_it.shift()))
│ │ │ │ +
│ │ │ │ +
298 if (
_g->interiorborder[codim].inside(
_it.coord(),
_it.shift()))
│ │ │ │ +
│ │ │ │ +
300 if (
_g->overlap[codim].inside(
_it.coord(),
_it.shift()))
│ │ │ │ +
│ │ │ │ +
302 if (
_g->overlapfront[codim].inside(
_it.coord(),
_it.shift()))
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
307 typedef typename GridImp::YGridLevelIterator
YGLI;
│ │ │ │ +
308 typedef typename GridImp::YGrid::Iterator
I;
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
333 typename std::deque<CommPartner>::const_iterator i;
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
336 std::array<int,dim> size;
│ │ │ │ +
│ │ │ │ +
338 for (
int i=0; i<dim; i++)
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
341 size[i] =
_g->mg->levelSize(
_g->level(), i);
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
354 for (
int i=dim-1; i>=0; i--)
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
365 task.buffer = buffer;
│ │ │ │ -
│ │ │ │ -
367 if (
rank!=_comm.rank())
│ │ │ │ -
368 _sendrequests.push_back(task);
│ │ │ │ -
│ │ │ │ -
370 _localsendrequests.push_back(task);
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
378 task.buffer = buffer;
│ │ │ │ -
│ │ │ │ -
380 if (
rank!=_comm.rank())
│ │ │ │ -
381 _recvrequests.push_back(task);
│ │ │ │ -
│ │ │ │ -
383 _localrecvrequests.push_back(task);
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
390 if (_localsendrequests.size()!=_localrecvrequests.size())
│ │ │ │ -
│ │ │ │ -
392 std::cout <<
"[" <<
rank() <<
"]: ERROR: local sends/receives do not match in exchange!" << std::endl;
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
395 for (
unsigned int i=0; i<_localsendrequests.size(); i++)
│ │ │ │ -
│ │ │ │ -
397 if (_localsendrequests[i].size!=_localrecvrequests[i].size)
│ │ │ │ -
│ │ │ │ -
399 std::cout <<
"[" <<
rank() <<
"]: ERROR: size in local sends/receive does not match in exchange!" << std::endl;
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
402 memcpy(_localrecvrequests[i].buffer,_localsendrequests[i].buffer,_localsendrequests[i].size);
│ │ │ │ -
│ │ │ │ -
404 _localsendrequests.clear();
│ │ │ │ -
405 _localrecvrequests.clear();
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
410 std::vector<MPI_Request> requests(_sendrequests.size() + _recvrequests.size());
│ │ │ │ -
411 MPI_Request* req = requests.data();
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
366 return _it.superindex();
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
374 std::bitset<dim-codim> subent_shift = Dune::Yasp::entityShift<dim-codim>(i,cc-codim);
│ │ │ │ +
375 std::bitset<dim-codim> subent_move = Dune::Yasp::entityMove<dim-codim>(i,cc-codim);
│ │ │ │ +
│ │ │ │ +
377 std::bitset<dim> shift =
_it.shift();
│ │ │ │ +
378 std::array<int, dim> coord =
_it.coord();
│ │ │ │ +
379 for (
int j=0, k=0; j<dim; j++)
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
384 coord[j] += subent_move[k];
│ │ │ │ +
385 shift[j] = subent_shift[k];
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
389 int which =
_g->overlapfront[cc].shiftmapping(shift);
│ │ │ │ +
390 return _g->overlapfront[cc].superindex(coord,which);
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
405 template<
int dim,
class Gr
idImp>
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
409 constexpr static int dimworld = GridImp::dimensionworld;
│ │ │ │ +
│ │ │ │ +
411 typedef typename GridImp::Traits::template Codim< 0 >::GeometryImpl
GeometryImpl;
│ │ │ │
│ │ │ │ -
│ │ │ │ -
414 for (
unsigned int i=0; i<_sendrequests.size(); i++)
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
419 MPI_Isend(_sendrequests[i].buffer, _sendrequests[i].size, MPI_BYTE,
│ │ │ │ -
420 _sendrequests[i].
rank, _tag, _comm, req++);
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
424 for (
unsigned int i=0; i<_recvrequests.size(); i++)
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
429 MPI_Irecv(_recvrequests[i].buffer, _recvrequests[i].size, MPI_BYTE,
│ │ │ │ -
430 _recvrequests[i].
rank, _tag, _comm, req++);
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
434 MPI_Waitall(requests.size(), requests.data(), MPI_STATUSES_IGNORE);
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
437 _sendrequests.clear();
│ │ │ │ -
438 _recvrequests.clear();
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
446 _comm.template allreduce<Dune::Max<double>,
double>(&x, &res, 1);
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
413 template<
int, PartitionIteratorType,
typename>
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
420 typedef typename GridImp::ctype
ctype;
│ │ │ │ +
│ │ │ │ +
422 typedef typename GridImp::YGridLevelIterator
YGLI;
│ │ │ │ +
423 typedef typename GridImp::YGrid::Iterator
I;
│ │ │ │ +
│ │ │ │ +
425 typedef typename GridImp::template Codim< 0 >::Geometry
Geometry;
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
435 typedef typename GridImp::template Codim<0>::EntitySeed
EntitySeed;
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
445 typedef typename GridImp::YGrid::iTupel
iTupel;
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
453 s <<
"[" <<
rank() <<
"]: Torus " <<
procs() <<
" processor(s) arranged as " <<
dims() << std::endl;
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
456 s <<
"[" <<
rank() <<
"]: send to "
│ │ │ │ -
457 <<
"rank=" << i.rank()
│ │ │ │ -
458 <<
" index=" << i.index()
│ │ │ │ -
459 <<
" delta=" << i.delta() <<
" dist=" << i.distance() << std::endl;
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
463 s <<
"[" <<
rank() <<
"]: recv from "
│ │ │ │ -
464 <<
"rank=" << i.rank()
│ │ │ │ -
465 <<
" index=" << i.index()
│ │ │ │ -
466 <<
" delta=" << i.delta() <<
" dist=" << i.distance() << std::endl;
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
478 std::fill(delta.begin(), delta.end(), -1);
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
487 for (
int i=0; i<d; i++)
│ │ │ │ -
488 nb[i] = ( me[i]+_dims[i]+delta[i] ) % _dims[i];
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
482 if (
_g->interior[0].inside(
_it.coord(),
_it.shift()))
│ │ │ │ +
│ │ │ │ +
484 if (
_g->overlap[0].inside(
_it.coord(),
_it.shift()))
│ │ │ │ +
│ │ │ │ +
486 DUNE_THROW(
GridError,
"Impossible GhostEntity");
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
494 for (
int i=0; i<d; i++)
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
500 _recvlist.push_back(cp);
│ │ │ │ -
501 cp.index = last-index;
│ │ │ │ -
502 _sendlist.push_front(cp);
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
509 for (
int i=0; i<d; i++)
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
493 auto ll =
_it.lowerleft();
│ │ │ │ +
494 auto ur =
_it.upperright();
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
497 for (
int i=0; i<dimworld; i++) {
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
501 auto size =
_g->mg->domainSize()[i];
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
505 auto size =
_g->mg->domainSize()[i];
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
512 GeometryImpl _geometry(ll,ur);
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
519 constexpr GeometryType
type ()
const
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
529 std::deque<CommPartner> _sendlist;
│ │ │ │ -
530 std::deque<CommPartner> _recvlist;
│ │ │ │ -
│ │ │ │ -
532 mutable std::vector<CommTask> _sendrequests;
│ │ │ │ -
533 mutable std::vector<CommTask> _recvrequests;
│ │ │ │ -
534 mutable std::vector<CommTask> _localsendrequests;
│ │ │ │ -
535 mutable std::vector<CommTask> _localrecvrequests;
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
540 template <
class Communication,
int d>
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
This file provides tools to partition YaspGrids. If you want to write your own partitioner,...
│ │ │ │ -
│ │ │ │ -
std::ostream & operator<<(std::ostream &out, const PartitionType &type)
write a PartitionType to a stream
Definition gridenums.hh:72
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
530 return Dune::Yasp::subEnt<dim>(dim,cc);
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
539 return Dune::Yasp::subEnt<dim>(dim,codim);
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
548 std::bitset<dim> move = Dune::Yasp::entityMove<dim>(i,cc);
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
552 for (
int j=0; j<dim; j++)
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
556 int which =
_g->overlapfront[cc].shiftmapping(Dune::Yasp::entityShift<dim>(i,cc));
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
565 DUNE_THROW(
GridError,
"tried to call father on level 0");
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
575 for (
int k=0; k<dim; k++) coord[k] = coord[k]/2;
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
583 return (
_g->level()>0);
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
591 FieldVector<ctype,dim> ll(0.0),ur(0.5);
│ │ │ │ +
│ │ │ │ +
593 for (
int k=0; k<dim; k++)
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
613 return (
_g->level() ==
yaspgrid()->maxLevel());
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
694 for (
int i=dim-1; i>=0; i--)
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
706 return _it.superindex();
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
713 std::bitset<dim> shift = Dune::Yasp::entityShift<dim>(i,cc);
│ │ │ │ +
714 std::bitset<dim> move = Dune::Yasp::entityMove<dim>(i,cc);
│ │ │ │ +
│ │ │ │ +
716 int trailing = (cc == dim) ? 1000 : 0;
│ │ │ │ +
│ │ │ │ +
718 std::array<int,dim> size =
_g->mg->levelSize(
_g->level());
│ │ │ │ +
719 std::array<int, dim> coord =
_it.coord();
│ │ │ │ +
720 for (
int j=0; j<dim; j++)
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
731 for (
int j=0; j<dim; j++)
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
737 for (
int k=0; k<
_g->level(); k++)
│ │ │ │ +
738 if (coord[j] & (1<<k))
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
742 trailing = std::min(trailing,zeroes);
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
754 for (
int j=dim-1; j>=0; j--)
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
767 std::bitset<dim> shift = Dune::Yasp::entityShift<dim>(i,cc);
│ │ │ │ +
768 std::bitset<dim> move = Dune::Yasp::entityMove<dim>(i,cc);
│ │ │ │ +
│ │ │ │ +
770 std::array<int, dim> coord =
_it.coord();
│ │ │ │ +
771 for (
int j=0; j<dim; j++)
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
774 int which =
_g->overlapfront[cc].shiftmapping(shift);
│ │ │ │ +
775 return _g->overlapfront[cc].superindex(coord,which);
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
784 template<
int dim,
class Gr
idImp>
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
788 constexpr static int dimworld = GridImp::dimensionworld;
│ │ │ │ +
│ │ │ │ +
790 template<
int, PartitionIteratorType,
typename>
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
793 typedef typename GridImp::Traits::template Codim<dim>::GeometryImpl
GeometryImpl;
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
796 typedef typename GridImp::ctype
ctype;
│ │ │ │ +
│ │ │ │ +
798 typedef typename GridImp::YGridLevelIterator
YGLI;
│ │ │ │ +
799 typedef typename GridImp::YGrid::Iterator
I;
│ │ │ │ +
│ │ │ │ +
801 typedef typename GridImp::template Codim<dim>::Geometry
Geometry;
│ │ │ │ +
│ │ │ │ +
803 typedef typename GridImp::template Codim<dim>::EntitySeed
EntitySeed;
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
809 typedef typename GridImp::YGrid::iTupel
iTupel;
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
845 return Dune::Yasp::subEnt<dim>(dim-dim,cc-dim);
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
850 GeometryImpl _geometry((
_it).lowerleft());
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
857 constexpr GeometryType
type ()
const
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
865 if (
_g->interior[dim].inside(
_it.coord(),
_it.shift()))
│ │ │ │ +
│ │ │ │ +
867 if (
_g->interiorborder[dim].inside(
_it.coord(),
_it.shift()))
│ │ │ │ +
│ │ │ │ +
869 if (
_g->overlap[dim].inside(
_it.coord(),
_it.shift()))
│ │ │ │ +
│ │ │ │ +
871 if (
_g->overlapfront[dim].inside(
_it.coord(),
_it.shift()))
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
892 iTupel size =
_g->mg->levelSize(
_g->level());
│ │ │ │ +
│ │ │ │ +
894 for (
int i=0; i<dim; i++)
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
902 for (
int i=0; i<dim; i++)
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
906 for (
int j=0; j<
_g->level(); j++)
│ │ │ │ +
907 if (
_it.coord(i)&(1<<j))
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
911 trailing = std::min(trailing,zeros);
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
915 int level =
_g->level()-trailing;
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
925 for (
int i=dim-1; i>=0; i--)
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
PartitionType
Attributes used in the generic overlap model.
Definition gridenums.hh:30
│ │ │ │ +
@ FrontEntity
on boundary between overlap and ghost
Definition gridenums.hh:34
│ │ │ │ +
@ InteriorEntity
all interior entities
Definition gridenums.hh:31
│ │ │ │ +
@ GhostEntity
ghost entities
Definition gridenums.hh:35
│ │ │ │ +
@ BorderEntity
on boundary between interior and overlap
Definition gridenums.hh:32
│ │ │ │ +
@ OverlapEntity
all entities lying in the overlap zone
Definition gridenums.hh:33
│ │ │ │ +
│ │ │ │
Include standard header files.
Definition agrid.hh:60
│ │ │ │ -
a base class for the yaspgrid partitioning strategy
Definition partitioning.hh:38
│ │ │ │ -
virtual void partition(const iTupel &, int, iTupel &, int) const =0
│ │ │ │ -
│ │ │ │ -
Torus()
constructor making uninitialized object
Definition torus.hh:67
│ │ │ │ -
int color(int rank) const
assign color to given rank
Definition torus.hh:197
│ │ │ │ -
double partition(int rank, iTupel origin_in, iTupel size_in, iTupel &origin_out, iTupel &size_out) const
partition the given grid onto the torus and return the piece of the process with given rank; returns ...
Definition torus.hh:239
│ │ │ │ -
int dims(int i) const
return dimensions of torus in direction i
Definition torus.hh:118
│ │ │ │ -
iTupel coord() const
return own coordinates
Definition torus.hh:100
│ │ │ │ -
int rank() const
return own rank
Definition torus.hh:94
│ │ │ │ -
const iTupel & dims() const
return dimensions of torus
Definition torus.hh:112
│ │ │ │ -
int rank_relative(int rank, int dir, int cnt) const
return rank of process where its coordinate in direction dir has offset cnt (handles periodic case)
Definition torus.hh:166
│ │ │ │ -
void recv(int rank, void *buffer, int size) const
store a receive request; buffers are received in order; handles also local requests with memcpy
Definition torus.hh:374
│ │ │ │ -
void send(int rank, void *buffer, int size) const
store a send request; buffers are sent in order; handles also local requests with memcpy
Definition torus.hh:361
│ │ │ │ -
Torus(Communication comm, int tag, iTupel size, int overlap, const Yasp::Partitioning< d > *partitioner)
make partitioner from communicator and coarse mesh size
Definition torus.hh:71
│ │ │ │ -
int neighbors() const
return the number of neighbors, which is
Definition torus.hh:203
│ │ │ │ -
void print(std::ostream &s) const
print contents of torus object
Definition torus.hh:451
│ │ │ │ -
double global_max(double x) const
global max
Definition torus.hh:443
│ │ │ │ -
Communication comm() const
return communicator
Definition torus.hh:124
│ │ │ │ -
int color(const iTupel &coord) const
assign color to given coordinate
Definition torus.hh:174
│ │ │ │ -
ProcListIterator recvend() const
last process in receive list
Definition torus.hh:355
│ │ │ │ -
bool is_neighbor(iTupel delta, std::bitset< d > periodic) const
return true if neighbor with given delta is a neighbor under the given periodicity
Definition torus.hh:212
│ │ │ │ -
ProcListIterator sendend() const
end of send list
Definition torus.hh:343
│ │ │ │ -
int procs() const
return number of processes
Definition torus.hh:106
│ │ │ │ -
iTupel rank_to_coord(int rank) const
map rank to coordinate in torus using lexicographic ordering
Definition torus.hh:144
│ │ │ │ -
std::array< int, d > iTupel
type used to pass tupels in and out
Definition torus.hh:49
│ │ │ │ -
ProcListIterator sendbegin() const
first process in send list
Definition torus.hh:337
│ │ │ │ -
void exchange() const
exchange messages stored in request buffers; clear request buffers afterwards
Definition torus.hh:387
│ │ │ │ -
ProcListIterator recvbegin() const
first process in receive list
Definition torus.hh:349
│ │ │ │ -
int coord_to_rank(iTupel coord) const
map coordinate in torus to rank using lexicographic ordering
Definition torus.hh:157
│ │ │ │ -
bool inside(iTupel c) const
return true if coordinate is inside torus
Definition torus.hh:136
│ │ │ │ -
int tag() const
return tag used by torus
Definition torus.hh:130
│ │ │ │ -
│ │ │ │ -
iTupel delta() const
return distance vector
Definition torus.hh:291
│ │ │ │ -
bool operator==(const ProcListIterator &iter) const
Return true when two iterators point to same member.
Definition torus.hh:313
│ │ │ │ -
bool operator!=(const ProcListIterator &iter) const
Return true when two iterators do not point to same member.
Definition torus.hh:320
│ │ │ │ -
int rank() const
return rank of neighboring process
Definition torus.hh:285
│ │ │ │ -
ProcListIterator(typename std::deque< CommPartner >::const_iterator iter)
make an iterator
Definition torus.hh:279
│ │ │ │ -
ProcListIterator & operator++()
Increment iterator to next cell.
Definition torus.hh:326
│ │ │ │ -
int index() const
return index in proclist
Definition torus.hh:297
│ │ │ │ -
int distance() const
return 1-norm of distance vector
Definition torus.hh:303
│ │ │ │ +
const int yaspgrid_level_bits
Definition yaspgrid.hh:48
│ │ │ │ +
const int yaspgrid_dim_bits
Definition yaspgrid.hh:47
│ │ │ │ +
Default Implementations for EntityImp.
Definition common/entity.hh:542
│ │ │ │ +
Base class for exceptions in Dune grid modules.
Definition exceptions.hh:20
│ │ │ │ +
static constexpr int mydimension
geometry dimension
Definition common/geometry.hh:94
│ │ │ │ +
The general version that handles all codimensions but 0 and dim.
Definition yaspgridgeometry.hh:31
│ │ │ │ +
Definition yaspgridentity.hh:242
│ │ │ │ +
int level() const
level of this element
Definition yaspgridentity.hh:256
│ │ │ │ +
Geometry geometry() const
geometry of this entity
Definition yaspgridentity.hh:270
│ │ │ │ +
GridImp::ctype ctype
Definition yaspgridentity.hh:248
│ │ │ │ +
constexpr GeometryType type() const
Return the name of the reference element. The type can be used to access the Dune::ReferenceElement.
Definition yaspgridentity.hh:279
│ │ │ │ +
PersistentIndexType persistentIndex() const
globally unique, persistent index
Definition yaspgridentity.hh:333
│ │ │ │ +
const YGLI & gridlevel() const
Definition yaspgridentity.hh:394
│ │ │ │ +
const I & transformingsubiterator() const
Definition yaspgridentity.hh:393
│ │ │ │ +
YGLI _g
Definition yaspgridentity.hh:400
│ │ │ │ +
bool equals(const YaspEntity &e) const
Return true when two iterators over the same grid are equal (!).
Definition yaspgridentity.hh:321
│ │ │ │ +
PartitionType partitionType() const
return partition type attribute
Definition yaspgridentity.hh:294
│ │ │ │ +
GridImp::Traits::template Codim< codim >::GeometryImpl GeometryImpl
Definition yaspgridentity.hh:251
│ │ │ │ +
YGLI & gridlevel()
Definition yaspgridentity.hh:396
│ │ │ │ +
int subCompressedIndex(int i, unsigned int cc) const
subentity compressed index
Definition yaspgridentity.hh:370
│ │ │ │ +
GridImp::YGridLevelIterator YGLI
Definition yaspgridentity.hh:307
│ │ │ │ +
GridImp::YGrid::Iterator I
Definition yaspgridentity.hh:308
│ │ │ │ +
const GridImp * yaspgrid() const
Definition yaspgridentity.hh:397
│ │ │ │ +
GridImp::template Codim< codim >::EntitySeed EntitySeed
Definition yaspgridentity.hh:253
│ │ │ │ +
YaspEntity(YGLI &&g, const I &&it)
Definition yaspgridentity.hh:316
│ │ │ │ +
EntitySeed seed() const
Return the entity seed which contains sufficient information to generate the entity again and uses as...
Definition yaspgridentity.hh:264
│ │ │ │ +
int compressedIndex() const
consecutive, codim-wise, level-wise index
Definition yaspgridentity.hh:364
│ │ │ │ +
YaspEntity(const YGLI &g, const I &it)
Definition yaspgridentity.hh:312
│ │ │ │ +
I _it
Definition yaspgridentity.hh:399
│ │ │ │ +
YaspEntity()
Definition yaspgridentity.hh:309
│ │ │ │ +
GridImp::template Codim< codim >::Geometry Geometry
Definition yaspgridentity.hh:250
│ │ │ │ +
GridImp::PersistentIndexType PersistentIndexType
Definition yaspgridentity.hh:330
│ │ │ │ +
unsigned int subEntities(unsigned int cc) const
Definition yaspgridentity.hh:288
│ │ │ │ +
I & transformingsubiterator()
Definition yaspgridentity.hh:395
│ │ │ │ +
Describes the minimal information necessary to create a fully functional YaspEntity.
Definition yaspgridentityseed.hh:18
│ │ │ │ +
Iterates over entities of one grid level.
Definition yaspgridleveliterator.hh:19
│ │ │ │ +
YaspIntersectionIterator enables iteration over intersections with neighboring codim 0 entities.
Definition yaspgridintersectioniterator.hh:22
│ │ │ │ +
YaspHierarchicIterator enables iteration over son entities of codim 0.
Definition yaspgridhierarchiciterator.hh:20
│ │ │ │ +
Implementation of Level- and LeafIndexSets for YaspGrid.
Definition yaspgridindexsets.hh:25
│ │ │ │ +
persistent, globally unique Ids
Definition yaspgrididset.hh:25
│ │ │ │ +
LeafIntersectionIterator ileafbegin() const
returns intersection iterator for first intersection
Definition yaspgridentity.hh:631
│ │ │ │ +
GridImp::template Codim< 0 >::LocalGeometry LocalGeometry
Definition yaspgridentity.hh:426
│ │ │ │ +
Codim< cc >::Entity subEntity(int i) const
Definition yaspgridentity.hh:545
│ │ │ │ +
bool equals(const YaspEntity &e) const
Return true when two iterators over the same grid are equal (!).
Definition yaspgridentity.hh:464
│ │ │ │ +
YaspEntity(const YGLI &g, I &&it)
Definition yaspgridentity.hh:455
│ │ │ │ +
int count() const
Definition yaspgridentity.hh:528
│ │ │ │ +
GridImp::template Codim< 0 >::EntitySeed EntitySeed
Definition yaspgridentity.hh:435
│ │ │ │ +
GridImp::PersistentIndexType PersistentIndexType
define the type used for persistent indices
Definition yaspgridentity.hh:442
│ │ │ │ +
GridImp::YGrid::iTupel iTupel
define type used for coordinates in grid module
Definition yaspgridentity.hh:445
│ │ │ │ +
const YGLI & gridlevel() const
Definition yaspgridentity.hh:606
│ │ │ │ +
IntersectionIterator iend() const
Reference to one past the last neighbor.
Definition yaspgridentity.hh:644
│ │ │ │ +
const I & transformingsubiterator() const
Definition yaspgridentity.hh:605
│ │ │ │ +
GridImp::ctype ctype
Definition yaspgridentity.hh:420
│ │ │ │ +
YGLI & gridlevel()
Definition yaspgridentity.hh:608
│ │ │ │ +
YaspEntity(const YGLI &g, const I &it)
Definition yaspgridentity.hh:451
│ │ │ │ +
const GridImp * yaspgrid() const
Definition yaspgridentity.hh:609
│ │ │ │ +
int level() const
level of this element
Definition yaspgridentity.hh:470
│ │ │ │ +
GridImp::LeafIntersectionIterator LeafIntersectionIterator
Definition yaspgridentity.hh:438
│ │ │ │ +
YaspEntity(YGLI &&g, I &&it)
Definition yaspgridentity.hh:459
│ │ │ │ +
PartitionType partitionType() const
return partition type attribute
Definition yaspgridentity.hh:480
│ │ │ │ +
LevelIntersectionIterator ilevelend() const
Reference to one past the last neighbor.
Definition yaspgridentity.hh:656
│ │ │ │ +
bool mightVanish() const
Returns true, if entity might disappear during the next call to adapt()
Definition yaspgridentity.hh:622
│ │ │ │ +
GridImp::HierarchicIterator HierarchicIterator
Definition yaspgridentity.hh:439
│ │ │ │ +
GridImp::template Codim< 0 >::Entity Entity
Definition yaspgridentity.hh:434
│ │ │ │ +
LevelIntersectionIterator ilevelbegin() const
returns intersection iterator for first intersection
Definition yaspgridentity.hh:638
│ │ │ │ +
LeafIntersectionIterator ileafend() const
Reference to one past the last neighbor.
Definition yaspgridentity.hh:650
│ │ │ │ +
YaspEntity()
Definition yaspgridentity.hh:448
│ │ │ │ +
GridImp::YGrid::Iterator I
Definition yaspgridentity.hh:423
│ │ │ │ +
GridImp::template Codim< 0 >::Geometry Geometry
Definition yaspgridentity.hh:425
│ │ │ │ +
Geometry geometry() const
geometry of this entity
Definition yaspgridentity.hh:491
│ │ │ │ +
Entity father() const
Inter-level access to father element on coarser grid. Assumes that meshes are nested.
Definition yaspgridentity.hh:561
│ │ │ │ +
bool isNew() const
Returns true, if the entity has been created during the last call to adapt()
Definition yaspgridentity.hh:618
│ │ │ │ +
constexpr GeometryType type() const
Return the name of the reference element. The type can be used to access the Dune::ReferenceElement.
Definition yaspgridentity.hh:519
│ │ │ │ +
EntitySeed seed() const
Return the entity seed which contains sufficient information to generate the entity again and uses as...
Definition yaspgridentity.hh:475
│ │ │ │ +
HierarchicIterator hbegin(int maxlevel) const
Definition yaspgridentity.hh:665
│ │ │ │ +
IntersectionIterator ibegin() const
returns intersection iterator for first intersection
Definition yaspgridentity.hh:625
│ │ │ │ +
GridImp::YGridLevelIterator YGLI
Definition yaspgridentity.hh:422
│ │ │ │ +
unsigned int subEntities(unsigned int codim) const
Definition yaspgridentity.hh:537
│ │ │ │ +
LocalGeometry geometryInFather() const
Definition yaspgridentity.hh:588
│ │ │ │ +
GridImp::LevelIntersectionIterator IntersectionIterator
Definition yaspgridentity.hh:436
│ │ │ │ +
GridImp::LevelIntersectionIterator LevelIntersectionIterator
Definition yaspgridentity.hh:437
│ │ │ │ +
I & transformingsubiterator()
Definition yaspgridentity.hh:607
│ │ │ │ +
bool isLeaf() const
Definition yaspgridentity.hh:611
│ │ │ │ +
bool hasFather() const
returns true if father entity exists
Definition yaspgridentity.hh:581
│ │ │ │ +
GridImp::template Codim< cd >::Entity Entity
Definition yaspgridentity.hh:431
│ │ │ │ +
int level() const
level of this element
Definition yaspgridentity.hh:830
│ │ │ │ +
bool equals(const YaspEntity &e) const
Return true when two iterators over the same grid are equal (!).
Definition yaspgridentity.hh:824
│ │ │ │ +
YGLI & gridlevel()
Definition yaspgridentity.hh:941
│ │ │ │ +
Geometry geometry() const
geometry of this entity
Definition yaspgridentity.hh:849
│ │ │ │ +
YaspEntity(const YGLI &g, const I &it)
Definition yaspgridentity.hh:815
│ │ │ │ +
unsigned int subEntities(unsigned int cc) const
Definition yaspgridentity.hh:843
│ │ │ │ +
PartitionType partitionType() const
return partition type attribute
Definition yaspgridentity.hh:863
│ │ │ │ +
GridImp::ctype ctype
Definition yaspgridentity.hh:796
│ │ │ │ +
GridImp::YGridLevelIterator YGLI
Definition yaspgridentity.hh:798
│ │ │ │ +
GridImp::PersistentIndexType PersistentIndexType
define the type used for persistent indices
Definition yaspgridentity.hh:806
│ │ │ │ +
EntitySeed seed() const
Return the entity seed which contains sufficient information to generate the entity again and uses as...
Definition yaspgridentity.hh:835
│ │ │ │ +
const GridImp * yaspgrid() const
Definition yaspgridentity.hh:943
│ │ │ │ +
constexpr GeometryType type() const
Return the name of the reference element. The type can be used to access the Dune::ReferenceElement.
Definition yaspgridentity.hh:857
│ │ │ │ +
GridImp::YGrid::Iterator I
Definition yaspgridentity.hh:799
│ │ │ │ +
I & transformingsubiterator()
Definition yaspgridentity.hh:940
│ │ │ │ +
const I & transformingsubiterator() const
Definition yaspgridentity.hh:938
│ │ │ │ +
GridImp::YGrid::iTupel iTupel
define type used for coordinates in grid module
Definition yaspgridentity.hh:809
│ │ │ │ +
GridImp::template Codim< dim >::Geometry Geometry
Definition yaspgridentity.hh:801
│ │ │ │ +
const YGLI & gridlevel() const
Definition yaspgridentity.hh:939
│ │ │ │ +
YaspEntity(YGLI &&g, I &&it)
Definition yaspgridentity.hh:819
│ │ │ │ +
YaspEntity()
Definition yaspgridentity.hh:812
│ │ │ │ +
GridImp::template Codim< dim >::EntitySeed EntitySeed
Definition yaspgridentity.hh:803
│ │ │ │
1.9.8
│ │ │ │
│ │ │ │
│ │ │ │
│ │ │ │ ├── html2text {}
│ │ │ │ │ @@ -2,38 +2,22 @@
│ │ │ │ │ Loading...
│ │ │ │ │ Searching...
│ │ │ │ │ No Matches
│ │ │ │ │ * _d_u_n_e
│ │ │ │ │ * _g_r_i_d
│ │ │ │ │ * _y_a_s_p_g_r_i_d
│ │ │ │ │ _C_l_a_s_s_e_s | _N_a_m_e_s_p_a_c_e_s
│ │ │ │ │ -yaspgridentity.hh File Reference
│ │ │ │ │ -the YaspEntity class and its specializations _M_o_r_e_._._.
│ │ │ │ │ -#include