potential.h

Go to the documentation of this file.

//  * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *  //
//                                                                                   //
//                       / /   / /    __    \ \  / /                                 //
//                      / /__ / /   / _ \    \ \/ /                                  //
//                     /  ___  /   | |/_/    / /\ \                                  //
//                    / /   / /    \_\      / /  \ \                                 //
//                                                                                   //
//                                                                                   //
//  Copyright (c) 2016, Jakub Benda, Charles University in Prague                    //
//                                                                                   //
// MIT License:                                                                      //
//                                                                                   //
//  Permission is hereby granted, free of charge, to any person obtaining a          //
// copy of this software and associated documentation files (the "Software"),        //
// to deal in the Software without restriction, including without limitation         //
// the rights to use, copy, modify, merge, publish, distribute, sublicense,          //
// and/or sell copies of the Software, and to permit persons to whom the             //
// Software is furnished to do so, subject to the following conditions:              //
//                                                                                   //
//  The above copyright notice and this permission notice shall be included          //
// in all copies or substantial portions of the Software.                            //
//                                                                                   //
//  THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS          //
// OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,       //
// FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE       //
// AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, //
// WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF         //
// OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.  //
//                                                                                   //
//  * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *  //

#ifndef HEX_DISTORTING_POTENTIAL
#define HEX_DISTORTING_POTENTIAL

#include "hex-arrays.h"
#include "hex-hydrogen.h"
#include "hex-special.h"

extern const rArrays Ucoeffs;

class DistortingPotential : public special::RadialFunction<double>
{
public:
    // constructors
    // @{
    DistortingPotential() : n_(0), k_(0.), rmax_(0) {}
    DistortingPotential(int n, double rmax = 0.) : n_(n), k_(0.), rmax_(rmax) {}
    DistortingPotential(double k, double rmax = 0.) : n_(0), k_(k), rmax_(rmax) {}
    DistortingPotential(DistortingPotential const & U) : n_(U.n_), k_(U.k_), rmax_(U.rmax_) {}
    // @}
    
    DistortingPotential operator= (DistortingPotential const& V);
    
    bool operator== (DistortingPotential const & V) const;
    
    double operator() (double x) const;
    
    double getTurningPoint () const { return 0.; }
    
    std::pair<double,int> getZeroAsymptotic (double x) const { return std::make_pair(-1,-1); }
    
    double plusMonopole(double x) const;
    
    double getConstant() const;
    
    double getFarRadius() const;
    
    void toFile(const char * filename) const;
    int n() const { return n_; }
    double k() const { return k_; }
    
private:
    int n_;        // principal quantum number of distorting state
    double k_;     // wavenumber of distorting state
    double rmax_;  // far radius
};

#endif