hamiltonians.c++

Go to the documentation of this file.

/*
 *

 Copyright (C) 2010 Joseph Pingenot

 This program is free software: you can redistribute it and/or modify
 it under the terms of the GNU Affero General Public License as published by
 the Free Software Foundation, either version 3 of the License, or
 (at your option) any later version.

 This program is distributed in the hope that it will be useful,
 but WITHOUT ANY WARRANTY; without even the implied warranty of
 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 GNU Affero General Public License for more details.

 You should have received a copy of the GNU Affero General Public License
 along with this program.  If not, see <http://www.gnu.org/licenses/>.

*/
#include <libkdotp/hamiltonians.h++>
#include <libmodelxx/grid.h++>
#include <libkdotp/hamiltonian_cartesian_effective_mass.h++>
#include <libkdotp/hamiltonian_cartesian_8band.h++>

namespace kdotp {
  namespace libkdotp {
    namespace hamiltonian {

      hamiltonian::hamiltonian(struct libmodelxx::structure::generic_structure *struc, enum model model_, boundary_condition bc, bool second_neighbor) : s(struc), model(model_) {
  switch(model) {
  case KDOTP_EFFECTIVE_MASS:
    switch(struc->g->t) {
    case libmodelxx::grid::CARTESIAN:
      specific_hamiltonian = new cartesian_effective_mass(struc, bc, second_neighbor);
      break;
    }
    break;
  case KDOTP_8BAND:
    switch(struc->g->t) {
    case libmodelxx::grid::CARTESIAN:
      specific_hamiltonian = new cartesian_8band(struc, bc, second_neighbor);
      break;
    }
    break;
  }
      }

      double* hamiltonian::get_eigenvectors(char VoI, double eVmin, double eVmax, int imin, int imax, double**evals, int *n_eigs, GError **err) const {
  switch(model) {
  case KDOTP_EFFECTIVE_MASS:
    switch(s->g->t) {
    case libmodelxx::grid::CARTESIAN:
      return ((cartesian_effective_mass*)specific_hamiltonian)->get_eigenvectors(VoI, eVmin, eVmax, imin, imax, evals, n_eigs, err);
      break;
    }
    break;
  case KDOTP_8BAND:
    switch(s->g->t) {
    case libmodelxx::grid::CARTESIAN:
      return ((cartesian_8band*)specific_hamiltonian)->get_eigenvectors(VoI, eVmin, eVmax, imin, imax, evals, n_eigs, err);
      break;
    }
    break;
  }
  return NULL;
      }
    }
  }
}