time_int Namespace Reference

Data Structures
class	Integrator

Functions
void	applyRigidBodyConstraint (data::ModelData &data)
void	updateCentralDifference (data::ModelData &data)
void	updateVerletHalfKickAndDrift (data::ModelData &data)
void	updateVerletSecondKick (data::ModelData &data)

Function Documentation

applyRigidBodyConstraint()

void time_int::applyRigidBodyConstraint

(

data::ModelData &

data

)

Definition at line 19 of file integrator.cpp.

                                                           {
  const auto &rigid = data.d_modelDeck_p->d_rigidParticles;
  if (rigid.empty())
    return;
 
  // Nodal force in the integrators is a force density (acceleration = f/rho).
  // A rigid translating body of mass M has acceleration a = sum(f_i V_i) / M,
  // so replacing each nodal force density by rho * a reproduces that motion
  // without touching the integrators themselves.
  for (const auto &[pid, mass] : rigid) {
    util::Point net(0., 0., 0.);
    for (std::size_t II = 0; II < data.numPdForceNodes(); II++) {
      const auto i = data.pdForceNode(II);
      if (data.getPtId(i) != pid)
        continue;
      net += data.getVol(i) * data.getF(i);
    }
 
    const auto acc = net / mass;
    for (std::size_t II = 0; II < data.numPdForceNodes(); II++) {
      const auto i = data.pdForceNode(II);
      if (data.getPtId(i) != pid)
        continue;
      data.getF(i) = data.getDensity(i) * acc;
    }
  }
}

Referenced by time_int::Integrator::integrate(), main(), and time_int::Integrator::step().

Here is the caller graph for this function:

updateCentralDifference()

void time_int::updateCentralDifference

(

data::ModelData &

data

)

Definition at line 47 of file integrator.cpp.

                                                          {
  const auto dim = data.dimension();
  const auto dt = data.currentDt();
 
  tf::Executor executor(util::parallel::getNThreads());
  tf::Taskflow taskflow;
 
  taskflow.for_each_index(
    (std::size_t) 0, data.numPdForceNodes(), (std::size_t) 1,
      [&data, dim, dt](std::size_t II) {
        const auto i = data.pdForceNode(II);
        const auto rho = data.getDensity(i);
        auto &v = data.getV(i);
        auto &u = data.getU(i);
        auto &x = data.getX(i);
        const auto &f = data.getF(i);
 
        for (int dof = 0; dof < dim; dof++) {
          if (data.isDofFree(i, dof)) {
            v[dof] += (dt / rho) * f[dof];
            u[dof] += dt * v[dof];
            x[dof] += dt * v[dof];
          }
        }
 
        data.setVMag(i, v.length());
      }
  );
 
  executor.run(taskflow).get();
}

References util::parallel::getNThreads().

Referenced by time_int::Integrator::step().

Here is the call graph for this function:

Here is the caller graph for this function:

updateVerletHalfKickAndDrift()

void time_int::updateVerletHalfKickAndDrift

(

data::ModelData &

data

)

Definition at line 79 of file integrator.cpp.

                                                               {
  const auto dim = data.dimension();
  const auto dt = data.currentDt();
 
  tf::Executor executor(util::parallel::getNThreads());
  tf::Taskflow taskflow;
 
  taskflow.for_each_index(
    (std::size_t) 0, data.numPdForceNodes(), (std::size_t) 1,
      [&data, dim, dt](std::size_t II) {
        const auto i = data.pdForceNode(II);
        const auto rho = data.getDensity(i);
        auto &v = data.getV(i);
        auto &u = data.getU(i);
        auto &x = data.getX(i);
        const auto &f = data.getF(i);
 
        for (int dof = 0; dof < dim; dof++) {
          if (data.isDofFree(i, dof)) {
            v[dof] += 0.5 * (dt / rho) * f[dof];
            u[dof] += dt * v[dof];
            x[dof] += dt * v[dof];
          }
        }
 
        data.setVMag(i, v.length());
      }
  );
 
  executor.run(taskflow).get();
}

References util::parallel::getNThreads().

Referenced by time_int::Integrator::step().

Here is the call graph for this function:

Here is the caller graph for this function:

updateVerletSecondKick()

void time_int::updateVerletSecondKick

(

data::ModelData &

data

)

Definition at line 111 of file integrator.cpp.

                                                         {
  const auto dim = data.dimension();
  const auto dt = data.currentDt();
 
  tf::Executor executor(util::parallel::getNThreads());
  tf::Taskflow taskflow;
 
  taskflow.for_each_index(
    (std::size_t) 0, data.numPdForceNodes(), (std::size_t) 1,
    [&data, dim, dt](std::size_t II) {
      const auto i = data.pdForceNode(II);
      const auto rho = data.getDensity(i);
      auto &v = data.getV(i);
      const auto &f = data.getF(i);
 
      for (int dof = 0; dof < dim; dof++) {
        if (data.isDofFree(i, dof)) {
          v[dof] += 0.5 * (dt / rho) * f[dof];
        }
      }
 
      data.setVMag(i, v.length());
    }
  );
 
  executor.run(taskflow).get();
}

References util::parallel::getNThreads().

Referenced by time_int::Integrator::step().

Here is the call graph for this function:

Here is the caller graph for this function: