15 Voronoi (amrex::Vector<amrex::Geometry> &_geom) :
IC(_geom) {}
20 Define(_number_of_grains,a_alpha);
22 Voronoi (amrex::Vector<amrex::Geometry> &a_geom,
int a_number_of_grains) :
IC(a_geom)
24 Define(a_number_of_grains,1.0);
29 Define(a_number_of_grains, std::vector<Set::Scalar>(a_number_of_grains,a_alpha),Type::Partition);
33 std::vector<Set::Scalar> a_alpha,
34 Type a_type = Type::Values)
53 AMREX_D_TERM(size(0) =
geom[0].ProbHi()[0] -
geom[0].ProbLo()[0];,
54 size(1) =
geom[0].ProbHi()[1] -
geom[0].ProbLo()[1];,
55 size(2) =
geom[0].ProbHi()[2] -
geom[0].ProbLo()[2];)
57 for (amrex::MFIter mfi(*a_field[lev],amrex::TilingIfNotGPU()); mfi.isValid(); ++mfi)
59 amrex::Box bx = mfi.tilebox();
60 bx.grow(a_field[lev]->nGrow());
61 int ncomp = a_field[lev]->nComp();
62 amrex::Array4<Set::Scalar>
const& field = a_field[lev]->array(mfi);
63 amrex::ParallelFor (bx,[=] AMREX_GPU_DEVICE(
int i,
int j,
int k) {
66 AMREX_D_TERM(x(0) =
geom[lev].ProbLo()[0] + ((amrex::Real)(i) + 0.5) *
geom[lev].CellSize()[0];,
67 x(1) =
geom[lev].ProbLo()[1] + ((amrex::Real)(j) + 0.5) *
geom[lev].CellSize()[1];,
68 x(2) =
geom[lev].ProbLo()[2] + ((amrex::Real)(k) + 0.5) *
geom[lev].CellSize()[2];);
70 amrex::Real min_distance = std::numeric_limits<amrex::Real>::infinity();
71 int min_grain_id = -1;
77 if (
geom[0].isPeriodic(0))
80 std::min( (x-
voronoi[n] + size(0)*Set::Vector::Unit(0)).lpNorm<2>(),
81 (x-
voronoi[n] - size(0)*Set::Vector::Unit(0)).lpNorm<2>()));
84 if (
geom[0].isPeriodic(1))
87 std::min( (x-
voronoi[n] + size(0)*Set::Vector::Unit(1)).lpNorm<2>(),
88 (x-
voronoi[n] - size(0)*Set::Vector::Unit(1)).lpNorm<2>()));
92 if (
geom[0].isPeriodic(2))
95 std::min( (x-
voronoi[n] + size(0)*Set::Vector::Unit(2)).lpNorm<2>(),
96 (x-
voronoi[n] - size(0)*Set::Vector::Unit(2)).lpNorm<2>()));
106 if (
type == Type::Values) field(i,j,k) =
alpha[min_grain_id];
107 else if (
type == Type::Partition) field(i,j,k,min_grain_id % ncomp) =
alpha[min_grain_id];