Inheritance diagram for anonymous_namespace{main.cpp}::GrainContactProbe:

Collaboration diagram for anonymous_namespace{main.cpp}::GrainContactProbe:

Public Member Functions
	GrainContactProbe (std::filesystem::path out_dir={}, size_t interval=0)

void	checkStop (data::ModelData &data) override

double	minGap () const

double	gap0 () const

double	Rc () const

double	tMin () const

int	maxPairs () const

Public Member Functions inherited from postprocess::Postprocess
virtual	~Postprocess ()=default

virtual void	close (data::ModelData &data)

virtual std::string	twoParticle (data::ModelData &data)

virtual std::string	compressive (data::ModelData &data)

Private Member Functions
void	sampleMetric (data::ModelData &data)

Static Private Member Functions
static bool	ownsNode (const data::ModelData &data, size_t i)

Private Attributes
double	d_minGap = 1.0e9

double	d_gap0 = 0.

double	d_Rc = 0.

double	d_tMin = 0.

int	d_maxPairs = 0

std::filesystem::path	d_outDir

size_t	d_interval = 0

std::ofstream	d_metricOs

bool	d_wroteXref = false

Detailed Description

Definition at line 302 of file main.cpp.

Constructor & Destructor Documentation

◆ GrainContactProbe()

anonymous_namespace{main.cpp}::GrainContactProbe::GrainContactProbe	(	std::filesystem::path	out_dir = `{}`,
		size_t	interval = `0`
	)

inline

Definition at line 304 of file main.cpp.

                                                {}, size_t interval = 0)
      : d_outDir(std::move(out_dir)), d_interval(interval) {
    if (!d_outDir.empty() && util::parallel::mpiRank() == 0) {
      std::filesystem::create_directories(d_outDir / "nodal");
      d_metricOs.open(d_outDir / "mpi_metric_ts.csv");
      d_metricOs << "step,t,max_u,com0x,com0y\n";
    }
  }

Member Function Documentation

◆ checkStop()

void anonymous_namespace{main.cpp}::GrainContactProbe::checkStop ( data::ModelData & data )

inlineoverridevirtual

Reimplemented from postprocess::Postprocess.

Definition at line 313 of file main.cpp.

                                               {
    postprocess::Postprocess::checkStop(data);
    sampleMetric(data);
    if (data.d_particlesListTypeAll.size() < 3)
      return;
    double Rc = 0.;
    try {
      Rc = data.d_particleDeck_p->d_contactDeck.getContact(0, 0).d_contactR;
    } catch (...) {
      return;
    }
    if (Rc <= 0.)
      return;
    d_Rc = Rc;
    double min_gap = 1.0e9;
    int n_in = 0;
    const auto &all = data.d_particlesListTypeAll;
    for (size_t i = 0; i < all.size(); ++i) {
      if (all[i]->getGroupId("contact_id") != 0)
        continue;
      for (size_t j = i + 1; j < all.size(); ++j) {
        if (all[j]->getGroupId("contact_id") != 0)
          continue;
        const double gap = all[i]->getXCenter().dist(all[j]->getXCenter()) -
                           all[i]->d_geom_p->boundingRadius() -
                           all[j]->d_geom_p->boundingRadius();
        min_gap = std::min(min_gap, gap);
        if (gap < Rc)
          ++n_in;
      }
    }
    if (d_gap0 == 0. && min_gap < 1.0e8)
      d_gap0 = min_gap;
    if (min_gap < d_minGap) {
      d_minGap = min_gap;
      d_tMin = data.d_time;
    }
    if (n_in > d_maxPairs) {
      d_maxPairs = n_in;
      util::io::print(std::format(
          "grain contact: {} pair(s) in Rc at t={:.6e}, min_gap={:.6e}, Rc={:.6e}\n",
          n_in, data.d_time, min_gap, Rc));
    }
  }

References postprocess::Postprocess::checkStop(), and util::io::print().

Here is the call graph for this function:

◆ gap0()

double anonymous_namespace{main.cpp}::GrainContactProbe::gap0 ( ) const

inline

Definition at line 359 of file main.cpp.

359{ return d_gap0; }

Referenced by main().

Here is the caller graph for this function:

◆ maxPairs()

int anonymous_namespace{main.cpp}::GrainContactProbe::maxPairs ( ) const

inline

Definition at line 362 of file main.cpp.

362{ return d_maxPairs; }

◆ minGap()

double anonymous_namespace{main.cpp}::GrainContactProbe::minGap ( ) const

inline

Definition at line 358 of file main.cpp.

358{ return d_minGap; }

◆ ownsNode()

static bool anonymous_namespace{main.cpp}::GrainContactProbe::ownsNode	(	const data::ModelData &	data,
		size_t	i
	)

inlinestaticprivate

Definition at line 365 of file main.cpp.

                                                            {
    const int rank = util::parallel::mpiRank();
    if (data.d_pdDofMpi) {
      if (i >= data.d_pdNodePartition.size())
        return false;
      return static_cast<int>(data.d_pdNodePartition[i]) == rank;
    }
    for (const auto *p : data.d_particlesListTypeAll) {
      const size_t i0 = p->d_globStart;
      const size_t i1 = i0 + p->getNumNodes();
      if (i < i0 || i >= i1)
        continue;
      if (p->isWall())
        return rank == 0;
      return particle::isLocallyOwned(*p);
    }
    return rank == 0;
  }

References particle::isLocallyOwned(), and util::parallel::mpiRank().

Here is the call graph for this function:

◆ Rc()

double anonymous_namespace{main.cpp}::GrainContactProbe::Rc ( ) const

inline

Definition at line 360 of file main.cpp.

360{ return d_Rc; }

◆ sampleMetric()

void anonymous_namespace{main.cpp}::GrainContactProbe::sampleMetric ( data::ModelData & data )

inlineprivate

Definition at line 384 of file main.cpp.

                                         {
    if (d_outDir.empty() || d_interval == 0)
      return;
    const size_t nstep = data.currentStep();
    if (nstep % d_interval != 0 && nstep < data.numTimeSteps())
      return;
    double max_u = 0.;
    for (const auto &u : data.d_u)
      max_u = std::max(max_u, u.length());
    if (util::parallel::mpiSize() > 1)
      MPI_Allreduce(MPI_IN_PLACE, &max_u, 1, MPI_DOUBLE, MPI_MAX,
                    util::parallel::mpiComm());
    double com0x = 0., com0y = 0.;
    if (!data.d_particlesListTypeParticle.empty()) {
      const auto c = data.d_particlesListTypeParticle[0]->getXCenter();
      com0x = c.d_x;
      com0y = c.d_y;
    }
    if (util::parallel::mpiRank() == 0) {
      d_metricOs << std::format("{},{:.12e},{:.12e},{:.12e},{:.12e}\n", nstep,
                                data.d_time, max_u, com0x, com0y);
      d_metricOs.flush();
    }
    const size_t n = data.d_u.size();
    std::vector<double> buf(6 * n, 0.);
    for (size_t i = 0; i < n; ++i) {
      if (!ownsNode(data, i))
        continue;
      buf[6 * i + 0] = data.d_u[i].d_x;
      buf[6 * i + 1] = data.d_u[i].d_y;
      buf[6 * i + 2] = data.d_u[i].d_z;
      buf[6 * i + 3] = data.d_v[i].d_x;
      buf[6 * i + 4] = data.d_v[i].d_y;
      buf[6 * i + 5] = data.d_v[i].d_z;
    }
    if (util::parallel::mpiSize() > 1)
      MPI_Allreduce(MPI_IN_PLACE, buf.data(), static_cast<int>(buf.size()),
                    MPI_DOUBLE, MPI_SUM, util::parallel::mpiComm());
    if (util::parallel::mpiRank() != 0)
      return;
    if (!d_wroteXref) {
      std::ofstream xr(d_outDir / "nodal" / "x_ref.bin", std::ios::binary);
      const uint32_t nn = static_cast<uint32_t>(n);
      xr.write(reinterpret_cast<const char *>(&nn), sizeof(nn));
      for (size_t i = 0; i < n; ++i) {
        const double p[3] = {data.d_xRef[i].d_x, data.d_xRef[i].d_y,
                             data.d_xRef[i].d_z};
        xr.write(reinterpret_cast<const char *>(p), sizeof(p));
      }
      d_wroteXref = true;
    }
    std::ofstream os(d_outDir / "nodal" / std::format("uv_{:06d}.bin", nstep),
                     std::ios::binary);
    const char magic[4] = {'P', 'D', 'U', 'V'};
    const uint32_t step32 = static_cast<uint32_t>(nstep);
    const uint32_t nn = static_cast<uint32_t>(n);
    os.write(magic, 4);
    os.write(reinterpret_cast<const char *>(&step32), sizeof(step32));
    os.write(reinterpret_cast<const char *>(&nn), sizeof(nn));
    os.write(reinterpret_cast<const char *>(buf.data()),
             static_cast<std::streamsize>(buf.size() * sizeof(double)));
  }

References util::parallel::mpiComm(), util::parallel::mpiRank(), util::parallel::mpiSize(), and anonymous_namespace{main.cpp}::ownsNode().

Here is the call graph for this function:

◆ tMin()

double anonymous_namespace{main.cpp}::GrainContactProbe::tMin ( ) const

inline

Definition at line 361 of file main.cpp.

361{ return d_tMin; }

Field Documentation

◆ d_gap0

double anonymous_namespace{main.cpp}::GrainContactProbe::d_gap0 = 0.

private

Definition at line 448 of file main.cpp.

◆ d_interval

size_t anonymous_namespace{main.cpp}::GrainContactProbe::d_interval = 0

private

Definition at line 453 of file main.cpp.

◆ d_maxPairs

int anonymous_namespace{main.cpp}::GrainContactProbe::d_maxPairs = 0

private

Definition at line 451 of file main.cpp.

◆ d_metricOs

std::ofstream anonymous_namespace{main.cpp}::GrainContactProbe::d_metricOs

private

Definition at line 454 of file main.cpp.

◆ d_minGap

double anonymous_namespace{main.cpp}::GrainContactProbe::d_minGap = 1.0e9

private

Definition at line 447 of file main.cpp.

◆ d_outDir

std::filesystem::path anonymous_namespace{main.cpp}::GrainContactProbe::d_outDir

private

Definition at line 452 of file main.cpp.

◆ d_Rc

double anonymous_namespace{main.cpp}::GrainContactProbe::d_Rc = 0.

private

Definition at line 449 of file main.cpp.

◆ d_tMin

double anonymous_namespace{main.cpp}::GrainContactProbe::d_tMin = 0.

private

Definition at line 450 of file main.cpp.

◆ d_wroteXref

bool anonymous_namespace{main.cpp}::GrainContactProbe::d_wroteXref = false

private

Definition at line 455 of file main.cpp.

The documentation for this class was generated from the following file:

/home/user/project/test/test_data/peridem/jha2021_comp_n50/main.cpp

Public Member Functions

Private Member Functions

Static Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ GrainContactProbe()

Member Function Documentation

◆ checkStop()

◆ gap0()

◆ maxPairs()

◆ minGap()

◆ ownsNode()

◆ Rc()

◆ sampleMetric()

◆ tMin()

Field Documentation

◆ d_gap0

◆ d_interval

◆ d_maxPairs

◆ d_metricOs

◆ d_minGap

◆ d_outDir

◆ d_Rc

◆ d_tMin

◆ d_wroteXref