Discontinuous Galerkin Library
#include "dg/algorithm.h"
Loading...
Searching...
No Matches
gridX.h
Go to the documentation of this file.
1#pragma once
2#include <cassert>
3#include <cmath> //floor function
4#include "dlt.h"
5#include "grid.h"
6#include "../enums.h"
7
45namespace dg{
47template<class real_type>
48struct RealGridX2d;
49template<class real_type>
50struct RealGridX3d;
52
66template<class real_type>
67struct RealGridX1d
68{
69 using value_type = real_type;
71 using host_vector = thrust::host_vector<real_type>;
72 using host_grid = RealGridX1d<real_type>;
74 constexpr static unsigned ndim() {return 1;}
75 RealGridX1d() = default;
86 RealGridX1d( real_type x0, real_type x1, real_type f, unsigned n, unsigned N, bc bcx = NEU):
87 x0_(x0), x1_(x1), f_(f),
88 n_(n), Nx_(N), bcx_(bcx)
89 {
90 assert( (f >= 0) && (f < 0.5) );
91 assert( fabs(outer_N() - f*(real_type)N) < 1e-14);
92 assert( x1 > x0 );
93 assert( N > 0 );
94 assert( n != 0 );
95 assert( bcx != PER);
96 }
97 host_vector abscissas(unsigned u=0) const { return grid().abscissas(0);}
98 host_vector weights(unsigned u=0) const { return grid().weights(0);}
99 unsigned shape(unsigned u=0) const { return grid().shape(0);}
100 auto get_shape() const { return grid().get_shape();}
106 real_type x0() const {return x0_;}
112 real_type x1() const {return x1_;}
118 real_type f() const {return f_;}
124 real_type lx() const {return x1_-x0_;}
130 real_type h() const {return lx()/(real_type)Nx_;}
136 unsigned N() const {return Nx_;}
142 unsigned outer_N() const {return (unsigned)(round(f_*(real_type)Nx_));}
148 unsigned inner_N() const {return N()-2*outer_N();}
154 unsigned n() const {return n_;}
160 bc bcx() const {return bcx_;}
166 unsigned size() const { return n_*Nx_;}
172 void display( std::ostream& os = std::cout) const
173 {
174 os << "RealGrid parameters are: \n"
175 <<" n = "<<n_<<"\n"
176 <<" N = "<<Nx_<<"\n"
177 <<" inner N = "<<inner_N()<<"\n"
178 <<" outer N = "<<outer_N()<<"\n"
179 <<" h = "<<h()<<"\n"
180 <<" x0 = "<<x0_<<"\n"
181 <<" x1 = "<<x1_<<"\n"
182 <<" lx = "<<lx()<<"\n"
183 <<"Boundary conditions in x are: \n"
184 <<" "<<bc2str(bcx_)<<"\n";
185 }
186 RealGrid<real_type,1> grid() const{return RealGrid<real_type,1>( x0_, x1_, n_, Nx_, bcx_);}
187
196 void shift_topologic( real_type x0, real_type& x1) const
197 {
198 assert( contains(x0));
199 real_type deltaX;
200 real_type xleft = x0_ + f_*lx();
201 real_type xright = x1_ - f_*lx();
202 if( x0 >= xleft && x0<xright)
203 {
204 if( x1 > xleft) deltaX = (x1 -xleft);
205 else deltaX = xright - x1;
206 unsigned N = floor(deltaX/(xright-xleft));
207 if( x1 > xright ) x1 -= N*lx();
208 if( x1 < xleft ) x1 += N*lx();
209 }
210 else if( x0 < xleft && x1 >=xleft)
211 x1 += (xright-xleft);
212 else if( x0 >= xright && x1 < xright)
213 x1 -= (xright-xleft);
214
215 }
216
225 bool contains( real_type x) const
226 {
227 if( (x>=x0_ && x <= x1_)) return true;
228 return false;
229 }
230 private:
231 real_type x0_, x1_, f_;
232 unsigned n_, Nx_;
233 bc bcx_;
234};
235
236//template<class real_type>
237//struct aRealTopologyX3d; //forward declare 3d version
238
255template<class real_type>
256struct aRealTopologyX2d
257{
258 using value_type = real_type;
260 using host_vector = thrust::host_vector<real_type>;
261 using host_grid = RealGridX2d<real_type>;
263 constexpr static unsigned ndim() {return 2;}
264 host_vector abscissas(unsigned u=0) const { return grid().abscissas(u);}
265 host_vector weights(unsigned u=0) const { return grid().weights(u);}
266 unsigned shape(unsigned u=0) const { return grid().shape(u);}
267 auto get_shape() const { return grid().get_shape();}
268
274 real_type x0() const {return x0_;}
280 real_type x1() const {return x1_;}
286 real_type y0() const {return y0_;}
292 real_type y1() const {return y1_;}
298 real_type lx() const {return x1_-x0_;}
304 real_type ly() const {return y1_-y0_;}
310 real_type hx() const {return lx()/(real_type)Nx_;}
316 real_type hy() const {return ly()/(real_type)Ny_;}
322 real_type fx() const {return fx_;}
328 real_type fy() const {return fy_;}
334 unsigned n() const {return n_;}
340 unsigned Nx() const {return Nx_;}
346 unsigned inner_Nx() const {return Nx_ - outer_Nx();}
352 unsigned outer_Nx() const {return (unsigned)round(fx_*(real_type)Nx_);}
358 unsigned Ny() const {return Ny_;}
364 unsigned inner_Ny() const {return Ny_-2*outer_Ny();}
370 unsigned outer_Ny() const {return (unsigned)round(fy_*(real_type)Ny_);}
376 bc bcx() const {return bcx_;}
382 bc bcy() const {return bcy_;}
388 RealGrid<real_type,2> grid() const {return RealGrid<real_type,2>( x0_,x1_,y0_,y1_,n_,Nx_,Ny_,bcx_,bcy_);}
394 unsigned size() const { return n_*n_*Nx_*Ny_;}
400 void display( std::ostream& os = std::cout) const
401 {
402 os << "Grid parameters are: \n"
403 <<" n = "<<n_<<"\n"
404 <<" Nx = "<<Nx_<<"\n"
405 <<" inner Nx = "<<inner_Nx()<<"\n"
406 <<" outer Nx = "<<outer_Nx()<<"\n"
407 <<" Ny = "<<Ny_<<"\n"
408 <<" inner Ny = "<<inner_Ny()<<"\n"
409 <<" outer Ny = "<<outer_Ny()<<"\n"
410 <<" hx = "<<hx()<<"\n"
411 <<" hy = "<<hy()<<"\n"
412 <<" x0 = "<<x0_<<"\n"
413 <<" x1 = "<<x1_<<"\n"
414 <<" y0 = "<<y0_<<"\n"
415 <<" y1 = "<<y1_<<"\n"
416 <<" lx = "<<lx()<<"\n"
417 <<" ly = "<<ly()<<"\n"
418 <<"Boundary conditions in x are: \n"
419 <<" "<<bc2str(bcx_)<<"\n"
420 <<"Boundary conditions in y are: \n"
421 <<" "<<bc2str(bcy_)<<"\n";
422 }
423
434 void shift_topologic( real_type x0, real_type y0, real_type& x1, real_type& y1) const
435 {
436 assert( contains(x0, y0));
437 real_type deltaX;
438 if( x1 > x0_) deltaX = (x1 -x0_);
439 else deltaX = x1_ - x1;
440 unsigned N = floor(deltaX/lx());
441 if( x1 > x1_ && bcx_ == dg::PER) x1 -= N*lx();
442 if( x1 < x0_ && bcx_ == dg::PER) x1 += N*lx();
443
444 if( x0 < x1_ - fx_*(x1_-x0_) ) //if x0 is one of the inner points
445 {
446 real_type deltaY;
447 real_type yleft = y0_ + fy_*ly();
448 real_type yright = y1_ - fy_*ly();
449 if( y0 >= yleft && y0<yright)
450 {
451 if( y1 > yleft) deltaY = (y1 -yleft);
452 else deltaY = yright - y1;
453 unsigned N = floor(deltaY/(yright-yleft));
454 if( y1 > yright ) y1 -= N*ly();
455 if( y1 < yleft ) y1 += N*ly();
456 }
457 else if( y0 < yleft && y1 >=yleft)
458 y1 += (yright-yleft);
459 else if( y0 >= yright && y1 < yright)
460 y1 -= (yright-yleft);
461 }
462
463 }
464
474 bool contains( real_type x, real_type y)const
475 {
476 if( (x>=x0_ && x <= x1_) && (y>=y0_ && y <= y1_)) return true;
477 return false;
478 }
479 protected:
481 ~aRealTopologyX2d() = default;
482 aRealTopologyX2d() = default;
485 aRealTopologyX2d( real_type x0, real_type x1, real_type y0, real_type y1, real_type fx, real_type fy, unsigned n, unsigned Nx, unsigned Ny, bc bcx, bc bcy):
486 x0_(x0), x1_(x1), y0_(y0), y1_(y1), fx_(fx), fy_(fy),
487 n_(n), Nx_(Nx), Ny_(Ny), bcx_(bcx), bcy_( bcy)
488 {
489 assert( (fy_ >= 0.) && (fy_ < 0.5) );
490 assert( (fx_ >= 0.) && (fx_ < 1.) );
491 assert( fabs(outer_Nx() - fx_*(real_type)Nx) < 1e-14);
492 assert( fabs(outer_Ny() - fy_*(real_type)Ny) < 1e-14);
493 assert( n != 0);
494 assert( x1 > x0 && y1 > y0);
495 assert( Nx_ > 0 && Ny > 0 );
496 assert( bcy != PER);
497 }
499 aRealTopologyX2d(const aRealTopologyX2d& src) = default;
501 aRealTopologyX2d& operator=(const aRealTopologyX2d& src) = default;
502 private:
503 real_type x0_, x1_, y0_, y1_;
504 real_type fx_, fy_;
505 unsigned n_, Nx_, Ny_;
506 bc bcx_, bcy_;
507};
512template<class real_type>
513struct RealGridX2d : public aRealTopologyX2d<real_type>
514{
517 RealGridX2d( real_type x0, real_type x1, real_type y0, real_type y1, real_type fx, real_type fy, unsigned n, unsigned Nx, unsigned Ny, bc bcx=PER, bc bcy=NEU):
518 aRealTopologyX2d<real_type>(x0,x1,y0,y1,fx,fy,n,Nx,Ny,bcx,bcy) { }
520 explicit RealGridX2d( const aRealTopologyX2d<real_type>& src): aRealTopologyX2d<real_type>(src){}
521};
522
539template<class real_type>
540struct aRealTopologyX3d
541{
542 using value_type = real_type;
544 using host_vector = thrust::host_vector<real_type>;
545 using host_grid = RealGridX3d<real_type>;
547 constexpr static unsigned ndim() {return 3;}
548 host_vector abscissas(unsigned u=0) const { return grid().abscissas(u);}
549 host_vector weights(unsigned u=0) const { return grid().weights(u);}
550 unsigned shape(unsigned u=0) const { return grid().shape(u);}
551 auto get_shape() const { return grid().get_shape();}
557 real_type x0() const {return x0_;}
563 real_type x1() const {return x1_;}
564
570 real_type y0() const {return y0_;}
576 real_type y1() const {return y1_;}
577
583 real_type z0() const {return z0_;}
589 real_type z1() const {return z1_;}
590
596 real_type lx() const {return x1_-x0_;}
602 real_type ly() const {return y1_-y0_;}
608 real_type lz() const {return z1_-z0_;}
609
615 real_type hx() const {return lx()/(real_type)Nx_;}
621 real_type hy() const {return ly()/(real_type)Ny_;}
627 real_type hz() const {return lz()/(real_type)Nz_;}
633 real_type fx() const {return fx_;}
639 real_type fy() const {return fy_;}
645 unsigned n() const {return n_;}
651 unsigned Nx() const {return Nx_;}
657 unsigned inner_Nx() const {return Nx_ - outer_Nx();}
663 unsigned outer_Nx() const {return (unsigned)round(fx_*(real_type)Nx_);}
669 unsigned Ny() const {return Ny_;}
675 unsigned inner_Ny() const {return Ny_-2*outer_Ny();}
681 unsigned outer_Ny() const {return (unsigned)round(fy_*(real_type)Ny_);}
687 unsigned Nz() const {return Nz_;}
693 bc bcx() const {return bcx_;}
699 bc bcy() const {return bcy_;}
705 bc bcz() const {return bcz_;}
711 RealGrid<real_type,3> grid() const {
712 return RealGrid<real_type,3>( x0_,x1_,y0_,y1_,z0_,z1_,n_,Nx_,Ny_,Nz_,bcx_,bcy_,bcz_);
713 }
719 unsigned size() const { return n_*n_*Nx_*Ny_*Nz_;}
725 void display( std::ostream& os = std::cout) const
726 {
727 os << "Grid parameters are: \n"
728 <<" n = "<<n_<<"\n"
729 <<" Nx = "<<Nx_<<"\n"
730 <<" inner Nx = "<<inner_Nx()<<"\n"
731 <<" outer Nx = "<<outer_Nx()<<"\n"
732 <<" Ny = "<<Ny_<<"\n"
733 <<" inner Ny = "<<inner_Ny()<<"\n"
734 <<" outer Ny = "<<outer_Ny()<<"\n"
735 <<" Nz = "<<Nz_<<"\n"
736 <<" hx = "<<hx()<<"\n"
737 <<" hy = "<<hy()<<"\n"
738 <<" hz = "<<hz()<<"\n"
739 <<" x0 = "<<x0_<<"\n"
740 <<" x1 = "<<x1_<<"\n"
741 <<" y0 = "<<y0_<<"\n"
742 <<" y1 = "<<y1_<<"\n"
743 <<" z0 = "<<z0_<<"\n"
744 <<" z1 = "<<z1_<<"\n"
745 <<" lx = "<<lx()<<"\n"
746 <<" ly = "<<ly()<<"\n"
747 <<" lz = "<<lz()<<"\n"
748 <<"Boundary conditions in x are: \n"
749 <<" "<<bc2str(bcx_)<<"\n"
750 <<"Boundary conditions in y are: \n"
751 <<" "<<bc2str(bcy_)<<"\n"
752 <<"Boundary conditions in z are: \n"
753 <<" "<<bc2str(bcz_)<<"\n";
754 }
765 bool contains( real_type x, real_type y, real_type z)const
766 {
767 if( (x>=x0_ && x <= x1_) && (y>=y0_ && y <= y1_) && (z>=z0_ && z<=z1_))
768 return true;
769 return false;
770 }
771 protected:
773 ~aRealTopologyX3d() = default;
774 aRealTopologyX3d() = default;
777 aRealTopologyX3d( real_type x0, real_type x1, real_type y0, real_type y1, real_type z0, real_type z1, real_type fx, real_type fy, unsigned n, unsigned Nx, unsigned Ny, unsigned Nz, bc bcx, bc bcy, bc bcz):
778 x0_(x0), x1_(x1), y0_(y0), y1_(y1), z0_(z0), z1_(z1), fx_(fx), fy_(fy),
779 n_(n), Nx_(Nx), Ny_(Ny), Nz_(Nz), bcx_(bcx), bcy_( bcy), bcz_( bcz)
780 {
781 assert( (fy_ >= 0.) && (fy_ < 0.5) );
782 assert( (fx_ >= 0.) && (fx_ < 1.) );
783 assert( fabs(outer_Nx() - fx_*(real_type)Nx) < 1e-14);
784 assert( fabs(outer_Ny() - fy_*(real_type)Ny) < 1e-14);
785 assert( n != 0);
786 assert( x1 > x0 && y1 > y0 ); assert( z1 > z0 );
787 assert( Nx_ > 0 && Ny > 0); assert( Nz > 0);
788 }
790 aRealTopologyX3d(const aRealTopologyX3d& src) = default;
792 aRealTopologyX3d& operator=(const aRealTopologyX3d& src) = default;
793 private:
794 real_type x0_, x1_, y0_, y1_, z0_, z1_;
795 real_type fx_,fy_;
796 unsigned n_, Nx_, Ny_, Nz_;
797 bc bcx_, bcy_, bcz_;
798};
799
804template<class real_type>
805struct RealGridX3d : public aRealTopologyX3d<real_type>
806{
809 RealGridX3d( real_type x0, real_type x1, real_type y0, real_type y1, real_type z0, real_type z1, real_type fx, real_type fy, unsigned n, unsigned Nx, unsigned Ny, unsigned Nz, bc bcx=PER, bc bcy=NEU, bc bcz=PER):
810 aRealTopologyX3d<real_type>(x0,x1,y0,y1,z0,z1,fx,fy,n,Nx,Ny,Nz,bcx,bcy,bcz) { }
812 explicit RealGridX3d( const aRealTopologyX3d<real_type>& src): aRealTopologyX3d<real_type>(src){}
813};
814
817using GridX1d = dg::RealGridX1d<double>;
818using GridX2d = dg::RealGridX2d<double>;
819using GridX3d = dg::RealGridX3d<double>;
820using aTopologyX2d = dg::aRealTopologyX2d<double>;
821using aTopologyX3d = dg::aRealTopologyX3d<double>;
823
824}// namespace dg
dlt.h
The discrete legendre trafo class.
enums.h
enums
grid.h
base topology classes
dg::bc2str
std::string bc2str(bc bcx)
write a string describing boundary condition to an output stream
Definition enums.h:37
dg::bc
bc
Switch between boundary conditions.
Definition enums.h:15
dg::coo3d::z
@ z
z direction
dg::PER
@ PER
periodic boundaries
Definition enums.h:16
dg::NEU
@ NEU
Neumann on both boundaries.
Definition enums.h:20
dg::coo2d::y
@ y
y direction
dg::coo2d::x
@ x
x direction
dg
This is the namespace for all functions and classes defined and used by the discontinuous Galerkin li...
dg::RealGrid
The simplest implementation of aRealTopology.
Definition grid.h:710
dg::RealGridX1d
1D grid for X-point topology
Definition gridX.h:68
dg::RealGridX1d::RealGridX1d
RealGridX1d(real_type x0, real_type x1, real_type f, unsigned n, unsigned N, bc bcx=NEU)
1D X-point grid
Definition gridX.h:86
dg::RealGridX1d::inner_N
unsigned inner_N() const
number of cells in the inner region
Definition gridX.h:148
dg::RealGridX1d::contains
bool contains(real_type x) const
Check if the grid contains a point.
Definition gridX.h:225
dg::RealGridX1d::bcx
bc bcx() const
boundary conditions
Definition gridX.h:160
dg::RealGridX1d::value_type
real_type value_type
Definition gridX.h:69
dg::RealGridX1d::get_shape
auto get_shape() const
Definition gridX.h:100
dg::RealGridX1d::weights
host_vector weights(unsigned u=0) const
Definition gridX.h:98
dg::RealGridX1d::size
unsigned size() const
the total number of points
Definition gridX.h:166
dg::RealGridX1d::abscissas
host_vector abscissas(unsigned u=0) const
Definition gridX.h:97
dg::RealGridX1d::RealGridX1d
RealGridX1d()=default
dg::RealGridX1d::outer_N
unsigned outer_N() const
number of cells in one of the outer regions
Definition gridX.h:142
dg::RealGridX1d::N
unsigned N() const
number of cells
Definition gridX.h:136
dg::RealGridX1d::shift_topologic
void shift_topologic(real_type x0, real_type &x1) const
Shifts a point coordinate due to topology.
Definition gridX.h:196
dg::RealGridX1d::display
void display(std::ostream &os=std::cout) const
Display.
Definition gridX.h:172
dg::RealGridX1d::x0
real_type x0() const
left boundary
Definition gridX.h:106
dg::RealGridX1d::n
unsigned n() const
number of polynomial coefficients
Definition gridX.h:154
dg::RealGridX1d::grid
RealGrid< real_type, 1 > grid() const
Definition gridX.h:186
dg::RealGridX1d::x1
real_type x1() const
right boundary
Definition gridX.h:112
dg::RealGridX1d::f
real_type f() const
Factor.
Definition gridX.h:118
dg::RealGridX1d::ndim
static constexpr unsigned ndim()
number of dimensions : 1
Definition gridX.h:74
dg::RealGridX1d::host_vector
thrust::host_vector< real_type > host_vector
The host vector type used by host functions like evaluate.
Definition gridX.h:71
dg::RealGridX1d::lx
real_type lx() const
total length of interval
Definition gridX.h:124
dg::RealGridX1d::h
real_type h() const
cell size
Definition gridX.h:130
dg::RealGridX1d::shape
unsigned shape(unsigned u=0) const
Definition gridX.h:99
dg::RealGridX2d
The simplest implementation of aRealTopologyX2d.
Definition gridX.h:514
dg::RealGridX2d::RealGridX2d
RealGridX2d(const aRealTopologyX2d< real_type > &src)
allow explicit type conversion from any other topology
Definition gridX.h:520
dg::RealGridX2d::RealGridX2d
RealGridX2d(real_type x0, real_type x1, real_type y0, real_type y1, real_type fx, real_type fy, unsigned n, unsigned Nx, unsigned Ny, bc bcx=PER, bc bcy=NEU)
Construct a 2D X-point grid.
Definition gridX.h:517
dg::RealGridX3d
The simplest implementation of aRealTopologyX3d.
Definition gridX.h:806
dg::RealGridX3d::RealGridX3d
RealGridX3d(real_type x0, real_type x1, real_type y0, real_type y1, real_type z0, real_type z1, real_type fx, real_type fy, unsigned n, unsigned Nx, unsigned Ny, unsigned Nz, bc bcx=PER, bc bcy=NEU, bc bcz=PER)
Construct a 3D X-point grid.
Definition gridX.h:809
dg::RealGridX3d::RealGridX3d
RealGridX3d(const aRealTopologyX3d< real_type > &src)
allow explicit type conversion from any other topology
Definition gridX.h:812
dg::aRealTopologyX2d
A 2D grid class with X-point topology.
Definition gridX.h:257
dg::aRealTopologyX2d::host_vector
thrust::host_vector< real_type > host_vector
The host vector type used by host functions like evaluate.
Definition gridX.h:260
dg::aRealTopologyX2d::size
unsigned size() const
real_typehe total number of points
Definition gridX.h:394
dg::aRealTopologyX2d::contains
bool contains(real_type x, real_type y) const
Check if the grid contains a point.
Definition gridX.h:474
dg::aRealTopologyX2d::hx
real_type hx() const
cell size in x
Definition gridX.h:310
dg::aRealTopologyX2d::weights
host_vector weights(unsigned u=0) const
Definition gridX.h:265
dg::aRealTopologyX2d::get_shape
auto get_shape() const
Definition gridX.h:267
dg::aRealTopologyX2d::aRealTopologyX2d
aRealTopologyX2d()=default
dg::aRealTopologyX2d::grid
RealGrid< real_type, 2 > grid() const
Return a copy without topology.
Definition gridX.h:388
dg::aRealTopologyX2d::shape
unsigned shape(unsigned u=0) const
Definition gridX.h:266
dg::aRealTopologyX2d::y0
real_type y0() const
left boundary in y
Definition gridX.h:286
dg::aRealTopologyX2d::bcy
bc bcy() const
boundary conditions in y
Definition gridX.h:382
dg::aRealTopologyX2d::lx
real_type lx() const
length of x
Definition gridX.h:298
dg::aRealTopologyX2d::y1
real_type y1() const
Right boundary in y.
Definition gridX.h:292
dg::aRealTopologyX2d::operator=
aRealTopologyX2d & operator=(const aRealTopologyX2d &src)=default
dg::aRealTopologyX2d::Nx
unsigned Nx() const
number of cells in x
Definition gridX.h:340
dg::aRealTopologyX2d::x1
real_type x1() const
Right boundary in x.
Definition gridX.h:280
dg::aRealTopologyX2d::x0
real_type x0() const
Left boundary in x.
Definition gridX.h:274
dg::aRealTopologyX2d::display
void display(std::ostream &os=std::cout) const
Display.
Definition gridX.h:400
dg::aRealTopologyX2d::fx
real_type fx() const
partition factor in x
Definition gridX.h:322
dg::aRealTopologyX2d::aRealTopologyX2d
aRealTopologyX2d(real_type x0, real_type x1, real_type y0, real_type y1, real_type fx, real_type fy, unsigned n, unsigned Nx, unsigned Ny, bc bcx, bc bcy)
Construct a 2D X-point grid.
Definition gridX.h:485
dg::aRealTopologyX2d::bcx
bc bcx() const
boundary conditions in x
Definition gridX.h:376
dg::aRealTopologyX2d::inner_Nx
unsigned inner_Nx() const
number of topological cells in x
Definition gridX.h:346
dg::aRealTopologyX2d::Ny
unsigned Ny() const
number of cells in y
Definition gridX.h:358
dg::aRealTopologyX2d::ly
real_type ly() const
length of y
Definition gridX.h:304
dg::aRealTopologyX2d::fy
real_type fy() const
partition factor in y
Definition gridX.h:328
dg::aRealTopologyX2d::outer_Ny
unsigned outer_Ny() const
number of cells in one of the outer regions of y
Definition gridX.h:370
dg::aRealTopologyX2d::hy
real_type hy() const
cell size in y
Definition gridX.h:316
dg::aRealTopologyX2d::n
unsigned n() const
number of polynomial coefficients in x and y
Definition gridX.h:334
dg::aRealTopologyX2d::outer_Nx
unsigned outer_Nx() const
number of smooth rows in x
Definition gridX.h:352
dg::aRealTopologyX2d::shift_topologic
void shift_topologic(real_type x0, real_type y0, real_type &x1, real_type &y1) const
Shifts a point coordinate due to topology.
Definition gridX.h:434
dg::aRealTopologyX2d::inner_Ny
unsigned inner_Ny() const
number of cells in the inner region of y
Definition gridX.h:364
dg::aRealTopologyX2d::aRealTopologyX2d
aRealTopologyX2d(const aRealTopologyX2d &src)=default
dg::aRealTopologyX2d::~aRealTopologyX2d
~aRealTopologyX2d()=default
disallow destruction through base class pointer
dg::aRealTopologyX2d::ndim
static constexpr unsigned ndim()
number of dimensions : 2
Definition gridX.h:263
dg::aRealTopologyX2d::value_type
real_type value_type
Definition gridX.h:258
dg::aRealTopologyX2d::abscissas
host_vector abscissas(unsigned u=0) const
Definition gridX.h:264
dg::aRealTopologyX3d
A 3D grid class with X-point topology.
Definition gridX.h:541
dg::aRealTopologyX3d::bcx
bc bcx() const
boundary conditions in x
Definition gridX.h:693
dg::aRealTopologyX3d::ndim
static constexpr unsigned ndim()
number of dimensions : 3
Definition gridX.h:547
dg::aRealTopologyX3d::operator=
aRealTopologyX3d & operator=(const aRealTopologyX3d &src)=default
dg::aRealTopologyX3d::n
unsigned n() const
number of polynomial coefficients in x and y
Definition gridX.h:645
dg::aRealTopologyX3d::y1
real_type y1() const
right boundary in y
Definition gridX.h:576
dg::aRealTopologyX3d::grid
RealGrid< real_type, 3 > grid() const
Return a copy without topology.
Definition gridX.h:711
dg::aRealTopologyX3d::z1
real_type z1() const
right boundary in z
Definition gridX.h:589
dg::aRealTopologyX3d::lz
real_type lz() const
length in z
Definition gridX.h:608
dg::aRealTopologyX3d::hx
real_type hx() const
cell size in x
Definition gridX.h:615
dg::aRealTopologyX3d::aRealTopologyX3d
aRealTopologyX3d(real_type x0, real_type x1, real_type y0, real_type y1, real_type z0, real_type z1, real_type fx, real_type fy, unsigned n, unsigned Nx, unsigned Ny, unsigned Nz, bc bcx, bc bcy, bc bcz)
Construct a 3D X-point grid.
Definition gridX.h:777
dg::aRealTopologyX3d::display
void display(std::ostream &os=std::cout) const
Display.
Definition gridX.h:725
dg::aRealTopologyX3d::fx
real_type fx() const
partition factor in x
Definition gridX.h:633
dg::aRealTopologyX3d::Nz
unsigned Nz() const
number of points in z
Definition gridX.h:687
dg::aRealTopologyX3d::get_shape
auto get_shape() const
Definition gridX.h:551
dg::aRealTopologyX3d::weights
host_vector weights(unsigned u=0) const
Definition gridX.h:549
dg::aRealTopologyX3d::inner_Ny
unsigned inner_Ny() const
number of cells in the inner region of y
Definition gridX.h:675
dg::aRealTopologyX3d::hy
real_type hy() const
cell size in y
Definition gridX.h:621
dg::aRealTopologyX3d::Ny
unsigned Ny() const
number of cells in y
Definition gridX.h:669
dg::aRealTopologyX3d::hz
real_type hz() const
cell size in z
Definition gridX.h:627
dg::aRealTopologyX3d::shape
unsigned shape(unsigned u=0) const
Definition gridX.h:550
dg::aRealTopologyX3d::Nx
unsigned Nx() const
number of points in x
Definition gridX.h:651
dg::aRealTopologyX3d::bcy
bc bcy() const
boundary conditions in y
Definition gridX.h:699
dg::aRealTopologyX3d::size
unsigned size() const
real_typehe total number of points
Definition gridX.h:719
dg::aRealTopologyX3d::~aRealTopologyX3d
~aRealTopologyX3d()=default
disallow destruction through base class pointer
dg::aRealTopologyX3d::z0
real_type z0() const
left boundary in z
Definition gridX.h:583
dg::aRealTopologyX3d::abscissas
host_vector abscissas(unsigned u=0) const
Definition gridX.h:548
dg::aRealTopologyX3d::y0
real_type y0() const
left boundary in y
Definition gridX.h:570
dg::aRealTopologyX3d::host_vector
thrust::host_vector< real_type > host_vector
The host vector type used by host functions like evaluate.
Definition gridX.h:544
dg::aRealTopologyX3d::inner_Nx
unsigned inner_Nx() const
number of topological cells in x
Definition gridX.h:657
dg::aRealTopologyX3d::value_type
real_type value_type
Definition gridX.h:542
dg::aRealTopologyX3d::x0
real_type x0() const
left boundary in x
Definition gridX.h:557
dg::aRealTopologyX3d::aRealTopologyX3d
aRealTopologyX3d()=default
dg::aRealTopologyX3d::outer_Ny
unsigned outer_Ny() const
number of cells in one of the outer regions of y
Definition gridX.h:681
dg::aRealTopologyX3d::x1
real_type x1() const
right boundary in x
Definition gridX.h:563
dg::aRealTopologyX3d::ly
real_type ly() const
length in y
Definition gridX.h:602
dg::aRealTopologyX3d::fy
real_type fy() const
partition factor in y
Definition gridX.h:639
dg::aRealTopologyX3d::bcz
bc bcz() const
boundary conditions in z
Definition gridX.h:705
dg::aRealTopologyX3d::contains
bool contains(real_type x, real_type y, real_type z) const
Check if the grid contains a point.
Definition gridX.h:765
dg::aRealTopologyX3d::aRealTopologyX3d
aRealTopologyX3d(const aRealTopologyX3d &src)=default
dg::aRealTopologyX3d::lx
real_type lx() const
length in x
Definition gridX.h:596
dg::aRealTopologyX3d::outer_Nx
unsigned outer_Nx() const
number of smooth rows in x
Definition gridX.h:663