message: Core: added atom type checking in Hamiltonian_Heisenberg_Pairs::Hessian.
diff: @@ -572,6 +572,8 @@ namespace Engine
  {
  int ispin = 3 * (exchange_pairs[i_pair].i + Vectormath::idx_from_translations(geometry->n_cells, geometry->n_spins_basic_domain, translations)) + alpha;
  int jspin = 3 * (exchange_pairs[i_pair].j + Vectormath::idx_from_translations(geometry->n_cells, geometry->n_spins_basic_domain, translations, exchange_pairs[i_pair].translations)) + alpha;
  + if ( check_atom_type(this->geometry->atom_types[ispin]) && check_atom_type(this->geometry->atom_types[jspin]) )
  + {
  hessian(ispin, jspin) += -exchange_magnitudes[i_pair];
  hessian(jspin, ispin) += -exchange_magnitudes[i_pair];
  }
  @@ -580,6 +582,7 @@ namespace Engine
  }
  }
  }
  + }
  // DMI
  for (unsigned int i_pair = 0; i_pair < this->dmi_pairs.size(); ++i_pair)
  @@ -599,6 +602,8 @@ namespace Engine
  {
  int ispin = 3 * (dmi_pairs[i_pair].i + Vectormath::idx_from_translations(geometry->n_cells, geometry->n_spins_basic_domain, translations)) + alpha;
  int jspin = 3 * (dmi_pairs[i_pair].j + Vectormath::idx_from_translations(geometry->n_cells, geometry->n_spins_basic_domain, translations, dmi_pairs[i_pair].translations)) + beta;
  + if ( check_atom_type(this->geometry->atom_types[ispin]) && check_atom_type(this->geometry->atom_types[jspin]) )
  + {
  if ((alpha == 0 && beta == 1))
  {
  @@ -649,6 +654,7 @@ namespace Engine
  }
  }
  }
  + }
  //// Dipole-Dipole
  //for (unsigned int i_pair = 0; i_pair < this->DD_indices.size(); ++i_pair)
  