wave_distort.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_DISTORTED_WAVE
#define HEX_DISTORTED_WAVE

#include <gsl/gsl_interp.h>

#include "hex-special.h"

#include "potential.h"

class DistortedWave : public special::RadialFunction<double>
{
public:
    
    // constructors
    // @{
    DistortedWave(double _kn, int _ln, DistortingPotential const & _U);
    DistortedWave(DistortedWave const& W) { *this = W; }
    // @}
    
    // destructor
    ~DistortedWave();
    
    // assignment
    DistortedWave operator= (DistortedWave const& W);
    
    double operator() (double x) const;
    
    Complex getPhasef() const { return Complex(cos(phase),sin(phase)); }
    
    double getPhase() const { return phase; }
    
    double k() const { return kn; }
    
    int l() const { return ln; }
    
    double getTurningPoint () const { return r0; }
    
    std::pair<double,int> getZeroAsymptotic (double x) const;
    
    void toFile(const char * filename) const { write_array(grid, array, filename); }
    
    double farRadius() const { return rf; }
    
    std::size_t sampleCount() const { return grid.size(); }
    
    mutable unsigned Evaluations;
    
private:
    
    // distorting potential
    DistortingPotential U;
    
    // interpolator
    gsl_interp *interpolator, *interpolator0;
    
    // distorted wave input parameters
    double kn;      // wavenumber of the distorted wave
    int ln;         // angular momentum of the distorted wave
    
    // classical turning point and far radius
    double r0, rf;
    
    // distorted wave computed attributes
    int samples, samples0;  // sample count
    double h, h0;           // discretization step
    
    // distorted wave data
    double phase;           // phase shift
    rArray grid, grid0;     // grid
    rArray array, array0;   // samples
};

#endif