include/qc2componentbec.h

/***************************************************************************
 *   Quantum Construct (qC++)                                                                  *
 *   The Quantum Physics Computational Library                                         *
 *   Copyright (C) 2005 by Shekhar S. Chandra                                      *
 *   Shekhar.Chandra@sci.monash,edu.au                                     *
 *                                                                         *
 *   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 2 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, write to the                         *
 *   Free Software Foundation, Inc.,                                       *
 *   59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.             *
 ***************************************************************************/
 /*
        Ver Info - QC2ComponentBEC.h
        1.0.0 - Original Implementation
 */

#ifndef QC2COMPONENTBEC_H
#define QC2COMPONENTBEC_H

#include "qcbec.h"

template <int rank>
class QC2ComponentBEC : public QCBEC<rank>
{
public:
// -------------------------------------------------------
    QC2ComponentBEC(void);
    ~QC2ComponentBEC(void);
        void setSize(int noOfRows, int noOfCols);

    // -------------------------------------------------------
    void setg11(parameterType g11)
    {   this->SystemParameters.setParameter(7,g11);     }
    void setg12(parameterType g12)
    {   this->SystemParameters.setParameter(8,g12);     }
    void setg22(parameterType g22)
    {   this->SystemParameters.setParameter(9,g22);     }
    void setg21(parameterType g21)
    {   this->SystemParameters.setParameter(10,g21);     }

        // -------------------------------------------------------
    parameterType getg11()
    {   return this->SystemParameters.getParameter(7);     }
    parameterType getg12()
    {   return this->SystemParameters.getParameter(8);     }
    parameterType getg22()
    {   return this->SystemParameters.getParameter(9);     }
    parameterType getg21()
    {   return this->SystemParameters.getParameter(10);     }
        Array<complexType,rank>* accessComp1()
        {   return this->OrderParameter.accessArray();     }
        Array<complexType,rank>* accessComp2()
        {   return OrderParameter2.accessArray();     }

    // -------------------------------------------------------
    void windPhase();
    void constantPhase();
    void randomPhase();
    void applyPhaseDisturbance();
    void evolve(bool forwards);
    void setupSteadyState();
    void writeResults();
        void writeImage(int shot, string filePrefix, string extension);

protected:
    QCComplexField<rank> OrderParameter2; 
};

#endif // QC2COMPONENTBEC_H

template <int rank>
QC2ComponentBEC<rank>::QC2ComponentBEC()
{
    this->SystemParameters.setNoOfParameters(11);
}

template <int rank>
QC2ComponentBEC<rank>::~QC2ComponentBEC()
{
    //dtor
}

template <int rank>
void QC2ComponentBEC<rank>::setSize(int noOfRows, int noOfCols)
{
    this->OrderParameter.setSize(noOfRows,noOfCols);
    OrderParameter2.setSize(noOfRows,noOfCols);
}

template <int rank>
void QC2ComponentBEC<rank>::windPhase()
{
    this->Phase.applyVortex(this->OrderParameter,this->deltaX,this->deltaY,this->xOffset,this->yOffset);
    this->Phase.applyVortex(OrderParameter2,this->deltaX,this->deltaY,this->xOffset,this->yOffset);
}

template <int rank>
void QC2ComponentBEC<rank>::constantPhase()
{
    this->Phase.applyUniform(this->OrderParameter,this->deltaX,this->deltaY,this->xOffset,this->yOffset);
    this->Phase.applyUniform(OrderParameter2,this->deltaX,this->deltaY,this->xOffset,this->yOffset);
}

template <int rank>
void QC2ComponentBEC<rank>::randomPhase()
{
    this->Phase.applyRandom(this->OrderParameter,this->deltaX,this->deltaY,this->xOffset,this->yOffset);
    this->Phase.applyRandom(OrderParameter2,this->deltaX,this->deltaY,this->xOffset,this->yOffset);
}

template <int rank>
void QC2ComponentBEC<rank>::applyPhaseDisturbance()
{
    this->Phase.applyDisturbanceTypeII(this->OrderParameter,this->deltaX,this->deltaY,this->xOffset,this->yOffset);
    this->Phase.applyDisturbanceTypeII(OrderParameter2,this->deltaX,this->deltaY,this->xOffset,this->yOffset);
}

template <int rank>
void QC2ComponentBEC<rank>::evolve(bool forwards)
{
    this->Evolve.evolveForwards(forwards);
    this->Evolve.usingTDGP2Component(this->OrderParameter,OrderParameter2,this->xOffset,this->yOffset,this->deltaX,this->deltaY,this->SystemParameters);
}

template <int rank>
void QC2ComponentBEC<rank>::setupSteadyState()
{
    this->SteadyState.usingTDGPVortexFixedPoint(this->OrderParameter,this->deltaR,this->SystemParameters);
        this->SteadyState.transferSolutionTo(this->OrderParameter,this->xOffset,this->yOffset,this->deltaX,this->deltaY,this->deltaR);
        this->SteadyState.usingTDGPVortexFixedPoint(OrderParameter2,this->deltaR,this->SystemParameters);
        this->SteadyState.transferSolutionTo(OrderParameter2,this->xOffset,this->yOffset,this->deltaX,this->deltaY,this->deltaR);
}

template <int rank>
void QC2ComponentBEC<rank>::writeResults()
{
    this->OrderParameter.writeNormToFile("FinalNormComp1.csv");
        this->Phase.outputPhaseFunction(this->OrderParameter,"FinalPhaseComp1.csv");
        OrderParameter2.writeNormToFile("FinalNormComp2.csv");
        this->Phase.outputPhaseFunction(OrderParameter2,"FinalPhaseComp2.csv");
}

template <int rank>
void QC2ComponentBEC<rank>::writeImage(int shot, string filePrefix, string extension)
{
    char digit = static_cast<char> (shot + 48);
    string filename1 = filePrefix + '1' + digit + extension;
    string filename2 = filePrefix + '2' + digit + extension;

    this->ptrField->writeToFile(filename1.c_str());
    OrderParameter2.writeToFile(filename2.c_str());
}

---