DiagonalizationSolver.h

Go to the documentation of this file.

/* Copyright 2016 Kristofer Björnson
 *
 * Licensed under the Apache License, Version 2.0 (the "License");
 * you may not use this file except in compliance with the License.
 * You may obtain a copy of the License at
 *
 *     http://www.apache.org/licenses/LICENSE-2.0
 *
 * Unless required by applicable law or agreed to in writing, software
 * distributed under the License is distributed on an "AS IS" BASIS,
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 * See the License for the specific language governing permissions and
 * limitations under the License.
 */

#ifndef COM_DAFER45_TBTK_DIAGONALIZATION_SOLVER
#define COM_DAFER45_TBTK_DIAGONALIZATION_SOLVER

#include "Model.h"
#include <complex>

namespace TBTK{

class DiagonalizationSolver{
public:
    DiagonalizationSolver();

    ~DiagonalizationSolver();

    void setModel(Model *model);

    void setSCCallback(
        bool (*scCallback)(
            DiagonalizationSolver *diagonalizationSolver
        )
    );

    void setMaxIterations(int maxIterations);

    void run();

    const double* getEigenValues();

    double* getEigenValuesRW();

    const std::complex<double>* getEigenVectors();

    std::complex<double>* getEigenVectorsRW();

    const double getEigenValue(int state);

    const std::complex<double> getAmplitude(int state, const Index &index);

    Model *getModel();
private:
    Model *model;

    std::complex<double> *hamiltonian;

    double *eigenValues;

    std::complex<double> *eigenVectors;

    int maxIterations;

    bool (*scCallback)(DiagonalizationSolver *diagonalizationSolver);

    void init();

    void update();

    void solve();
};

inline void DiagonalizationSolver::setModel(Model *model){
    this->model = model;
}

inline void DiagonalizationSolver::setSCCallback(
    bool (*scCallback)(
        DiagonalizationSolver *diagonalizationSolver
    )
){
    this->scCallback = scCallback;
}

inline void DiagonalizationSolver::setMaxIterations(int maxIterations){
    this->maxIterations = maxIterations;
}

inline const double* DiagonalizationSolver::getEigenValues(){
    return eigenValues;
}

inline double* DiagonalizationSolver::getEigenValuesRW(){
    return eigenValues;
}

inline const std::complex<double>* DiagonalizationSolver::getEigenVectors(){
    return eigenVectors;
}

inline std::complex<double>* DiagonalizationSolver::getEigenVectorsRW(){
    return eigenVectors;
}

inline const std::complex<double> DiagonalizationSolver::getAmplitude(
    int state,
    const Index &index
){
    return eigenVectors[model->getBasisSize()*state + model->getBasisIndex(index)];
}

inline const double DiagonalizationSolver::getEigenValue(int state){
    return eigenValues[state];
}

inline Model* DiagonalizationSolver::getModel(){
    return model;
}

};  //End of namespace TBTK

#endif