anonymous_namespace{main.cpp}::ContactWitness Class Reference

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

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

Public Member Functions
void	checkStop (data::ModelData &data) override
void	assertOk () 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 Attributes
bool	d_sawPp = false
bool	d_sawWall = false
double	d_minPpGap = 1.e300
double	d_minGrainY = 1.e300
double	d_Rc = 0.

Detailed Description

Records whether PP and wall contacts engaged at any step (floor top y=0).

Definition at line 119 of file main.cpp.

Member Function Documentation

assertOk()

void anonymous_namespace{main.cpp}::ContactWitness::assertOk ( ) const

inline

Definition at line 155 of file main.cpp.

                        {
    int saw_pp = d_sawPp ? 1 : 0;
    int saw_wall = d_sawWall ? 1 : 0;
    double min_pp = d_minPpGap;
    double min_y = d_minGrainY;
    double rc = d_Rc;
    if (util::parallel::mpiSize() > 1) {
      int saw_pp_g = 0, saw_wall_g = 0;
      double min_pp_g = 0., min_y_g = 0., rc_g = 0.;
      MPI_Allreduce(&saw_pp, &saw_pp_g, 1, MPI_INT, MPI_MAX,
                    util::parallel::mpiComm());
      MPI_Allreduce(&saw_wall, &saw_wall_g, 1, MPI_INT, MPI_MAX,
                    util::parallel::mpiComm());
      MPI_Allreduce(&min_pp, &min_pp_g, 1, MPI_DOUBLE, MPI_MIN,
                    util::parallel::mpiComm());
      MPI_Allreduce(&min_y, &min_y_g, 1, MPI_DOUBLE, MPI_MIN,
                    util::parallel::mpiComm());
      MPI_Allreduce(&rc, &rc_g, 1, MPI_DOUBLE, MPI_MAX,
                    util::parallel::mpiComm());
      saw_pp = saw_pp_g;
      saw_wall = saw_wall_g;
      min_pp = min_pp_g;
      min_y = min_y_g;
      rc = rc_g;
    }
    if (!saw_pp)
      throw std::runtime_error(std::format(
          "particle-particle contact missing: min_gap={:.6e} Rc={:.6e}", min_pp,
          rc));
    if (!saw_wall)
      throw std::runtime_error(std::format(
          "particle-wall contact missing: min_grain_y={:.6e} Rc={:.6e}", min_y,
          rc));
    if (util::parallel::mpiRank() == 0)
      util::io::print(std::format(
          "contacts OK: min_pp_gap={:.6e} min_grain_y={:.6e} Rc={:.6e}\n",
          min_pp, min_y, rc));
  }

References util::parallel::mpiComm(), util::parallel::mpiRank(), util::parallel::mpiSize(), and util::io::print().

Here is the call graph for this function:

checkStop()

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

inlineoverridevirtual

Reimplemented from postprocess::Postprocess.

Definition at line 121 of file main.cpp.

                                               {
    postprocess::Postprocess::checkStop(data);
    if (data.d_particlesListTypeParticle.size() < 2)
      return;
    const auto *g0 = data.d_particlesListTypeParticle[0];
    const auto *g1 = data.d_particlesListTypeParticle[1];
    double Rc = 0.;
    try {
      Rc = data.d_particleDeck_p->d_contactDeck.getContact(0, 0).d_contactR;
    } catch (...) {
      return;
    }
    d_Rc = Rc;
    const double gap = g0->getXCenter().dist(g1->getXCenter()) -
                       g0->d_geom_p->boundingRadius() -
                       g1->d_geom_p->boundingRadius();
    d_minPpGap = std::min(d_minPpGap, gap);
    if (gap < Rc)
      d_sawPp = true;
 
    double min_y = 1.e300;
    for (size_t i = 0; i < data.d_x.size(); ++i) {
      if (data.getParticleFromAllList(data.d_ptId[i])->isWall())
        continue;
      min_y = std::min(min_y, data.d_x[i].d_y);
    }
    if (util::parallel::mpiSize() > 1)
      MPI_Allreduce(MPI_IN_PLACE, &min_y, 1, MPI_DOUBLE, MPI_MIN,
                    util::parallel::mpiComm());
    d_minGrainY = std::min(d_minGrainY, min_y);
    if (min_y < Rc)
      d_sawWall = true;
  }

References postprocess::Postprocess::checkStop(), util::parallel::mpiComm(), and util::parallel::mpiSize().

Here is the call graph for this function:

Field Documentation

d_minGrainY

double anonymous_namespace{main.cpp}::ContactWitness::d_minGrainY = 1.e300

private

Definition at line 198 of file main.cpp.

d_minPpGap

double anonymous_namespace{main.cpp}::ContactWitness::d_minPpGap = 1.e300

private

Definition at line 197 of file main.cpp.

d_Rc

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

private

Definition at line 199 of file main.cpp.

d_sawPp

bool anonymous_namespace{main.cpp}::ContactWitness::d_sawPp = false

private

Definition at line 195 of file main.cpp.

d_sawWall

bool anonymous_namespace{main.cpp}::ContactWitness::d_sawWall = false

private

Definition at line 196 of file main.cpp.

---

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

• /home/user/project/test/test_data/peridem/mpi_identity_twop_wall/main.cpp