Method sigma_row

scheme/structure.cpp:1138–1182

}

void geom_epsilon::sigma_row(meep::component c, double sigrow[3], const meep::vec &r) {
  vector3 p = vec_to_vector3(r);

  boolean inobject;
  material_type material =
      material_of_unshifted_point_in_tree_inobject(p, restricted_tree, &inobject);

  int destroy_material = 0;
  if (material.which_subclass == MTS::MATERIAL_TYPE_SELF) { material = default_material; }
  if (material.which_subclass == MTS::MATERIAL_FUNCTION) {
    material = eval_material_func(material.subclass.material_function_data->material_func, p);
    destroy_material = 1;
  }

  sigrow[0] = sigrow[1] = sigrow[2] = 0.0;
  if (material.which_subclass == MTS::MEDIUM) {
    susceptibility_list slist = type(c) == meep::E_stuff
                                    ? material.subclass.medium_data->E_susceptibilities
                                    : material.subclass.medium_data->H_susceptibilities;
    for (int j = 0; j < slist.num_items; ++j)
      if (susceptibility_equiv(&slist.items[j], &current_pol->user_s)) {
        int ic = meep::component_index(c);
        switch (ic) { // which row of the sigma tensor to return
          case 0:
            sigrow[0] = slist.items[j].sigma_diag.x;
            sigrow[1] = slist.items[j].sigma_offdiag.x;
            sigrow[2] = slist.items[j].sigma_offdiag.y;
            break;
          case 1:
            sigrow[0] = slist.items[j].sigma_offdiag.x;
            sigrow[1] = slist.items[j].sigma_diag.y;
            sigrow[2] = slist.items[j].sigma_offdiag.z;
            break;
          default: // case 2:
            sigrow[0] = slist.items[j].sigma_offdiag.y;
            sigrow[1] = slist.items[j].sigma_offdiag.z;
            sigrow[2] = slist.items[j].sigma_diag.z;
            break;
        }
        break;
      }
  }

  if (destroy_material) material_type_destroy(material);
}

/* make multilevel_susceptibility from scheme input data */
static meep::susceptibility *make_multilevel_sus(const multilevel_atom *d) {