LP2R.h

Go to the documentation of this file.

#ifndef _LinLin_H_
#define _LinLin_H_
#include <iostream>
#include <fstream>
#include <sstream>
#include <vector>
#include <string>
#include <cmath>
#include <algorithm>
 
class C_LPoly{
public:
C_LPoly()=default;
C_LPoly(const double M, const double W, const double  M_e) :
 mass(M), wt(W), Z_chain(M/M_e), t_FRouse(Z_chain*Z_chain) {}
 
double mass=0.0; 
double wt=0.0; 
double Z_chain=0.0; 
double z=0.0; 
bool alive=true; 
bool relax_free_Rouse=false; 
bool rept_set=false; 
double tau_d_0=1.0e22; 
double Z_rept=0.0; 
double rept_wt=0.0; 
int p_max=0; 
int p_next=0; 
double t_FRouse; 
             };
 
 
class InvSqSum{ // sum(1/p^2)
double psum[500];
public:
  InvSqSum(){
   psum[0]=0.0; psum[1]=1.0;
   for(int i=2; i<500; i++){psum[i]=psum[i-1]+1.0/((double) (i*i));}
            }
  double intg_psum(int n1, int n2){
   double n1d=((double) n1); double n2d=((double) n2);
   double n1dsq=n1d*n1d; double n2dsq=n2d*n2d;
   double val=1.0/n1d - 1.0/n2d + 0.50*(1.0/n1dsq + 1.0/n2dsq);
     val+=(1.0/(n1d*n1dsq) - 1.0/(n2d*n2dsq))/6.0;
     return val;
                                  }
   double operator()(int Z){
     if(Z < 500){return psum[Z];}
     else{ return psum[499]+intg_psum(500,Z);}
                        }
   double operator()(int Z1, int Z2){
      double s1, s2;
      if(Z1 <= 500){s1=psum[Z1-1];}
      else{s1=psum[499]+intg_psum(500,Z1-1);}
      if(Z2 < 500){s2=psum[Z2];}
      else{s2=psum[499]+intg_psum(500,Z2);}
      return s2-s1;
                                 }
 
              };
 
#include "./LP2R_NS.h"
#include "./routines.h"
 
 
#endif