PhotonDiffusion::OpCameraInteg Struct Reference

Inheritance diagram for PhotonDiffusion::OpCameraInteg:

Collaboration diagram for PhotonDiffusion::OpCameraInteg:

Public Member Functions
	OpCameraInteg (boost::shared_ptr< CommonData > common_data_ptr)
MoFEMErrorCode	doWork (int side, EntityType type, EntitiesFieldData::EntData &data)
	[Integral_calc]

Public Attributes
boost::shared_ptr< CommonData >	commonDataPtr

Detailed Description

Definition at line 151 of file photon_diffusion.cpp.

Constructor & Destructor Documentation

OpCameraInteg()

PhotonDiffusion::OpCameraInteg::OpCameraInteg ( boost::shared_ptr< CommonData >  common_data_ptr )

inline

Definition at line 153 of file photon_diffusion.cpp.

        : BoundaryEleOp("PHOTON_FLUENCE_RATE", OPROW),
          commonDataPtr(common_data_ptr) {
      std::fill(&doEntities[MBVERTEX], &doEntities[MBMAXTYPE], false);
      doEntities[MBTRI] = doEntities[MBQUAD] = true;
    }

Member Function Documentation

doWork()

MoFEMErrorCode PhotonDiffusion::OpCameraInteg::doWork	(	int	side,
		EntityType	type,
		EntitiesFieldData::EntData &	data
	)

[Integral_calc]

Definition at line 954 of file photon_diffusion.cpp.

                                                                       {
  MoFEMFunctionBegin;
  const int nb_integration_pts = getGaussPts().size2();
  const double area = getMeasure();
  auto t_w = getFTensor0IntegrationWeight();
  auto t_val = getFTensor0FromVec(*(commonDataPtr->approxVals));
 
  double values_integ = 0;
 
#ifdef ENABLE_PYTHON_BINDING
  VectorDouble sens_vals_vec;
  if (enable_python) {
    MatrixDouble ref_coords = getCoordsAtGaussPts();
    sens_vals_vec =
        interp_function(interp_image_name, ref_coords, nb_integration_pts,
                        cam_len_x, cam_len_y, "(.*)INTERPOLATION(.*)");
  }
#endif
 
  for (int gg = 0; gg != nb_integration_pts; ++gg) {
 
    double sens = 1;
#ifdef ENABLE_PYTHON_BINDING
    if (enable_python) {
      sens = sens_vals_vec(gg);
    }
#endif
 
    const double alpha = t_w * area * sens;
 
    values_integ += alpha * t_val;
 
    ++t_w;
    ++t_val;
  }
 
  constexpr std::array<int, 1> indices = {CommonData::VALUES_INTEG};
  std::array<double, 1> values;
  values[0] = values_integ;
  CHKERR VecSetValues(commonDataPtr->petscVec, 1, indices.data(), values.data(),
                      ADD_VALUES);
  MoFEMFunctionReturn(0);
}

Member Data Documentation

commonDataPtr

boost::shared_ptr<CommonData> PhotonDiffusion::OpCameraInteg::commonDataPtr

Definition at line 152 of file photon_diffusion.cpp.

---

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

• tutorials/scl-10_photon_diffusion/photon_diffusion.cpp