QOperators.h

Go to the documentation of this file.

/****************************************************************************
*   CQEDSimulator, a dynamic library for C/C++ programs                     
*   Copyright (C) 2009 Stefano Carrazza. All rights reserved.               
*                                                                                                                                                       
*   This program is free software: you can redistribute it and/or modify    
*   it under the terms of the GNU 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 General Public License for more details.                                                        
*                                                                                                                                                       
*   You should have received a copy of the GNU General Public License       
*   along with this program.  If not, see <http://www.gnu.org/licenses/>.   
*****************************************************************************/

#ifndef QOPERATORS_H
#define QOPERATORS_H

#include "SMatrix.h"
#include "Bra.h"


class EXPORTED_FUNCTION QOperator: public SMatrix
{
public:
        QOperator(int dim_i = 1):SMatrix(dim_i) {} 
        QOperator(int dim_i,complex<double> diag):SMatrix(dim_i,diag) {} 
        QOperator(int dim_i,const complex<double> *Mtab_i):SMatrix(dim_i,Mtab_i) {} 
        QOperator(int dim_i,double *real,double *imag):SMatrix(dim_i,real,imag) {} 
        QOperator(const SMatrix & m):SMatrix(m) {} 
        QOperator(const QOperator & m):SMatrix(m) {} 

        //Operators 
        friend EXPORTED_FUNCTION Ket operator*(const QOperator & m,const Ket & ket); 
        friend EXPORTED_FUNCTION Bra operator*(const Bra & bra, const QOperator & m); 
};


class EXPORTED_FUNCTION QCreation: public QOperator
{
public:
        QCreation(int dim_i = 1); 
};


class EXPORTED_FUNCTION QAnnihilation: public QOperator
{
public:
        QAnnihilation(int dim_i = 1);
};


class EXPORTED_FUNCTION QDisplacement: public QOperator
{
private:
        QCreation *ac;
        QAnnihilation *aa;
public:
        QDisplacement(int dim_i = 1); 
        ~QDisplacement(); 

        QOperator operator() (complex<double> alpha) const {
          return exp(-pow(abs(alpha),2)/2.0)*(alpha*ac).getexpm()*(-conj(alpha)*aa).getexpm();
        } 
};


class EXPORTED_FUNCTION QParity: public QOperator
{
public:
        QParity(int dim_i = 1); 
};

#endif

---