hamiltonians.h++

Go to the documentation of this file.

/*
 *Central switching point for Hamiltonians.
 *

    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/>.

*/

#ifndef __HAMILTONIANS_HXX__
#define __HAMILTONIANS_HXX__

#include <libmodelxx/generic_structure.h++>

namespace kdotp {
  namespace libkdotp {
    namespace hamiltonian {

      enum gerror_domains {
  HAMILTONIAN_CARTESIAN_EFFECTIVE_MASS_GERROR_DOMAIN,
  HAMILTONIAN_CARTESIAN_8BAND_GERROR_DOMAIN
      };

      /*This is the list of all models that are supported.
       *The different models incorporate strain etc. differently.
       */
      enum model {
  KDOTP_EFFECTIVE_MASS,
  KDOTP_8BAND
      };

      enum boundary_condition {
  HARDWALL,
  PERIODIC
      };

      class hamiltonian {
      public:
  void *specific_hamiltonian;
  struct libmodelxx::structure::generic_structure *s;
  enum model model;

  /**Sets up the Hamiltonian
   *\param struct the structure to be turned into a Hamiltonian
   *\param model 
   */
  hamiltonian(struct libmodelxx::structure::generic_structure *struc, enum model model_, boundary_condition bc, bool second_neighbor);

  double* get_eigenvectors_double(double eVmin, double eVmax, double** evals, int *n_eigs, GError **err) const {
    return get_eigenvectors('V', eVmin, eVmax, 0, 0, evals, n_eigs, err);
  }
  double* get_eigenvectors_int(int imin, int imax, double** evals, int *n_eigs, GError **err) const {
    return get_eigenvectors('I', 0, 0, imin, imax, evals, n_eigs, err);
  }

  double* get_eigenvectors(char VoI, double eVmin, double eVmax, int imin, int imax, double**evals, int *n_eigs, GError **err) const;


      private:
  hamiltonian() {}
      };

    }
  }
}



#endif