anonymous_namespace{vtkParticleWriter.cpp} Namespace Reference

Functions
void	appendPointArraysForNodes (vtkUnstructuredGrid *grid, const data::ModelData *model, const std::vector< size_t > &gids, const std::vector< std::string > &tags)

Function Documentation

appendPointArraysForNodes()

void anonymous_namespace{vtkParticleWriter.cpp}::appendPointArraysForNodes	(	vtkUnstructuredGrid *	grid,
		const data::ModelData *	model,
		const std::vector< size_t > &	gids,
		const std::vector< std::string > &	tags
	)

Definition at line 389 of file vtkParticleWriter.cpp.

                                                                 {
  double value[3] = {0., 0., 0.};
  auto add_vec3 = [&](const char *name, const std::vector<util::Point> &field) {
    auto array = vtkSmartPointer<vtkDoubleArray>::New();
    array->SetNumberOfComponents(3);
    array->SetName(name);
    array->SetComponentName(0, "x");
    array->SetComponentName(1, "y");
    array->SetComponentName(2, "z");
    for (size_t g : gids) {
      const auto &ui = field[g];
      value[0] = ui.d_x;
      value[1] = ui.d_y;
      value[2] = ui.d_z;
      array->InsertNextTuple(value);
    }
    grid->GetPointData()->AddArray(array);
  };
  auto add_scalar = [&](const char *name, auto getter) {
    auto array = vtkSmartPointer<vtkDoubleArray>::New();
    array->SetNumberOfComponents(1);
    array->SetName(name);
    for (size_t g : gids) {
      value[0] = static_cast<double>(getter(g));
      array->InsertNextTuple(value);
    }
    grid->GetPointData()->AddArray(array);
  };
 
  if (util::methods::isTagInList("Displacement", tags))
    add_vec3("Displacement", model->d_u);
  if (util::methods::isTagInList("Velocity", tags))
    add_vec3("Velocity", model->d_v);
  if (util::methods::isTagInList("Force_Density", tags))
    add_vec3("Force_Density", model->d_f);
  if (util::methods::isTagInList("Force", tags)) {
    auto array = vtkSmartPointer<vtkDoubleArray>::New();
    array->SetNumberOfComponents(3);
    array->SetName("Force");
    array->SetComponentName(0, "x");
    array->SetComponentName(1, "y");
    array->SetComponentName(2, "z");
    for (size_t g : gids) {
      const auto &fi = model->d_f[g];
      const double vol = model->d_vol[g];
      value[0] = fi.d_x * vol;
      value[1] = fi.d_y * vol;
      value[2] = fi.d_z * vol;
      array->InsertNextTuple(value);
    }
    grid->GetPointData()->AddArray(array);
  }
  if (util::methods::isTagInList("Damage_Z", tags) && !model->d_Z.empty())
    add_scalar("Damage_Z", [&](size_t g) { return model->d_Z[g]; });
  if (util::methods::isTagInList("Damage", tags) && !model->d_phi.empty())
    add_scalar("Damage", [&](size_t g) { return model->d_phi[g]; });
  if (util::methods::isTagInList("Damage_Bond", tags) && !model->d_phiBond.empty())
    add_scalar("Damage_Bond", [&](size_t g) { return model->d_phiBond[g]; });
  if (util::methods::isTagInList("Particle_ID", tags))
    add_scalar("Particle_ID", [&](size_t g) {
      return static_cast<double>(
          model->getParticleFromAllList(model->d_ptId[g])->getId());
    });
}

References appendPointArraysForNodes(), data::ModelData::d_f, data::ModelData::d_phi, data::ModelData::d_phiBond, data::ModelData::d_ptId, data::ModelData::d_u, data::ModelData::d_v, data::ModelData::d_vol, data::ModelData::d_Z, particle::BaseParticle::getId(), data::ModelData::getParticleFromAllList(), and util::methods::isTagInList().

Referenced by appendPointArraysForNodes().

Here is the call graph for this function:

Here is the caller graph for this function: