Changeset dd60b45 in sasview for sansmodels/src/c_models


Ignore:
Timestamp:
Jan 5, 2012 11:01:26 AM (13 years ago)
Author:
Mathieu Doucet <doucetm@…>
Branches:
master, ESS_GUI, ESS_GUI_Docs, ESS_GUI_batch_fitting, ESS_GUI_bumps_abstraction, ESS_GUI_iss1116, ESS_GUI_iss879, ESS_GUI_iss959, ESS_GUI_opencl, ESS_GUI_ordering, ESS_GUI_sync_sascalc, costrafo411, magnetic_scatt, release-4.1.1, release-4.1.2, release-4.2.2, release_4.0.1, ticket-1009, ticket-1094-headless, ticket-1242-2d-resolution, ticket-1243, ticket-1249, ticket885, unittest-saveload
Children:
a8eab1c
Parents:
2d1b700
Message:

refl_adv model refactor

Location:
sansmodels/src/c_models
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • sansmodels/src/c_models/models.hh

    r2d1b700 rdd60b45  
    2727using namespace std; 
    2828 
    29 class ReflAdvModel{ 
    30 public: 
    31         // Model parameters 
    32         Parameter n_layers; 
    33         Parameter scale; 
    34         Parameter thick_inter0; 
    35         Parameter func_inter0; 
    36         Parameter sld_bottom0; 
    37         Parameter sld_medium; 
    38         Parameter background; 
    39  
    40         Parameter sld_flat1; 
    41         Parameter sld_flat2; 
    42         Parameter sld_flat3; 
    43         Parameter sld_flat4; 
    44         Parameter sld_flat5; 
    45         Parameter sld_flat6; 
    46         Parameter sld_flat7; 
    47         Parameter sld_flat8; 
    48         Parameter sld_flat9; 
    49         Parameter sld_flat10; 
    50  
    51         Parameter thick_inter1; 
    52         Parameter thick_inter2; 
    53         Parameter thick_inter3; 
    54         Parameter thick_inter4; 
    55         Parameter thick_inter5; 
    56         Parameter thick_inter6; 
    57         Parameter thick_inter7; 
    58         Parameter thick_inter8; 
    59         Parameter thick_inter9; 
    60         Parameter thick_inter10; 
    61  
    62         Parameter thick_flat1; 
    63         Parameter thick_flat2; 
    64         Parameter thick_flat3; 
    65         Parameter thick_flat4; 
    66         Parameter thick_flat5; 
    67         Parameter thick_flat6; 
    68         Parameter thick_flat7; 
    69         Parameter thick_flat8; 
    70         Parameter thick_flat9; 
    71         Parameter thick_flat10; 
    72  
    73         Parameter func_inter1; 
    74         Parameter func_inter2; 
    75         Parameter func_inter3; 
    76         Parameter func_inter4; 
    77         Parameter func_inter5; 
    78         Parameter func_inter6; 
    79         Parameter func_inter7; 
    80         Parameter func_inter8; 
    81         Parameter func_inter9; 
    82         Parameter func_inter10; 
    83  
    84         Parameter sldIM_flat1; 
    85         Parameter sldIM_flat2; 
    86         Parameter sldIM_flat3; 
    87         Parameter sldIM_flat4; 
    88         Parameter sldIM_flat5; 
    89         Parameter sldIM_flat6; 
    90         Parameter sldIM_flat7; 
    91         Parameter sldIM_flat8; 
    92         Parameter sldIM_flat9; 
    93         Parameter sldIM_flat10; 
    94  
    95         Parameter nu_inter1; 
    96         Parameter nu_inter2; 
    97         Parameter nu_inter3; 
    98         Parameter nu_inter4; 
    99         Parameter nu_inter5; 
    100         Parameter nu_inter6; 
    101         Parameter nu_inter7; 
    102         Parameter nu_inter8; 
    103         Parameter nu_inter9; 
    104         Parameter nu_inter10; 
    105  
    106         Parameter sldIM_sub0; 
    107         Parameter sldIM_medium; 
    108         Parameter npts_inter; 
    109         Parameter nu_inter0; 
    110  
    111         // Constructor 
    112         ReflAdvModel(); 
    113  
    114         // Operators to get I(Q) 
    115         double operator()(double q); 
    116         double operator()(double qx, double qy); 
    117         double calculate_ER(); 
    118         double evaluate_rphi(double q, double phi); 
    119 }; 
     29 
    12030 
    12131 

  • sansmodels/src/c_models/refl_adv.cpp

    r67424cd rdd60b45  
    11 
    22#include <math.h> 
    3 #include "models.hh" 
    43#include "parameters.hh" 
    54#include <stdio.h> 
     5#include <stdlib.h> 
     6#include "refl_adv.h" 
     7 
     8extern "C" { 
     9#include "libmultifunc/librefl.h" 
     10} 
     11 
    612using namespace std; 
    713 
    8 extern "C" { 
    9         #include "refl_adv.h" 
    10 } 
    11  
     14#define lamda 4.62 
     15 
     16 
     17double re_adv_kernel(double dp[], double q) { 
     18  int n = dp[0]; 
     19  int i,j; 
     20  double nsl; 
     21 
     22  double scale = dp[1]; 
     23  double thick_inter_sub = dp[2]; 
     24  double sld_sub = dp[4]; 
     25  double sld_super = dp[5]; 
     26  double background = dp[6]; 
     27  double npts = dp[69]; //number of sub_layers in each interface 
     28 
     29  double total_thick=0.0; 
     30 
     31  int n_s; 
     32  double sld_i,sldim_i,dz,phi,R,ko2; 
     33  double pi; 
     34 
     35  int* fun_type; 
     36  double* sld; 
     37  double* sld_im; 
     38  double* thick_inter; 
     39  double* thick; 
     40  double* fun_coef; 
     41  complex  phi1,alpha,alpha2,kn,fnm,fnp,rn,Xn,nn,nn2,an,nnp1,one,two,n_sub,n_sup,knp1,Xnp1; 
     42 
     43  fun_type = (int*)malloc((n+2)*sizeof(int)); 
     44  sld = (double*)malloc((n+2)*sizeof(double)); 
     45  sld_im = (double*)malloc((n+2)*sizeof(double)); 
     46  thick_inter = (double*)malloc((n+2)*sizeof(double)); 
     47  thick = (double*)malloc((n+2)*sizeof(double)); 
     48  fun_coef = (double*)malloc((n+2)*sizeof(double)); 
     49 
     50  fun_type[0] = dp[3]; 
     51  fun_coef[0] = fabs(dp[70]); 
     52  for (i =1; i<=n; i++){ 
     53    sld[i] = dp[i+6]; 
     54    thick_inter[i]= dp[i+16]; 
     55    thick[i] = dp[i+26]; 
     56    fun_type[i] = dp[i+36]; 
     57    sld_im[i] = dp[i+46]; 
     58    fun_coef[i] = fabs(dp[i+56]); 
     59    //printf("type_func2 =%g\n",fun_coef[i]); 
     60 
     61    total_thick += thick[i] + thick_inter[i]; 
     62  } 
     63  sld[0] = sld_sub; 
     64  sld[n+1] = sld_super; 
     65  sld_im[0] = fabs(dp[1+66]); 
     66  sld_im[n+1] = fabs(dp[2+66]); 
     67  thick[0] = total_thick/5.0; 
     68  thick[n+1] = total_thick/5.0; 
     69  thick_inter[0] = thick_inter_sub; 
     70  thick_inter[n+1] = 0.0; 
     71  fun_coef[n+1] = 0.0; 
     72 
     73  nsl=npts;//21.0; //nsl = Num_sub_layer:  MUST ODD number in double //no other number works now 
     74 
     75  pi = 4.0*atan(1.0); 
     76  one = cassign(1.0,0.0); 
     77  Xn = cassign(0.0,0.0); 
     78  two = cassign(0.0,-2.0); 
     79 
     80  //Checking if floor is available. 
     81  //no imaginary sld inputs in this function yet 
     82  n_sub=cassign(1.0-sld_sub*pow(lamda,2.0)/(2.0*pi),pow(lamda,2.0)/(2.0*pi)*sld_im[0]); 
     83  n_sup=cassign(1.0-sld_super*pow(lamda,2.0)/(2.0*pi),pow(lamda,2.0)/(2.0*pi)*sld_im[n+1]); 
     84  ko2 = pow(2.0*pi/lamda,2.0); 
     85 
     86  phi = asin(lamda*q/(4.0*pi)); 
     87  phi1 = cplx_div(rcmult(phi,one),n_sup); 
     88  alpha = cplx_mult(n_sup,cplx_cos(phi1)); 
     89  alpha2 = cplx_mult(alpha,alpha); 
     90 
     91  nnp1=n_sub; 
     92  knp1=cplx_sqrt(rcmult(ko2,cplx_sub(cplx_mult(nnp1,nnp1),alpha2)));  //nnp1*ko*sin(phinp1) 
     93  Xnp1=cassign(0.0,0.0); 
     94  dz = 0.0; 
     95  // iteration for # of layers +sub from the top 
     96  for (i=1;i<=n+1; i++){ 
     97    //if (fun_coef[i]==0.0) 
     98    //  // this condition protects an error in numerical multiplication 
     99    //  fun_coef[i] = 1e-14; 
     100    //iteration for 9 sub-layers 
     101    for (j=0;j<2;j++){ 
     102      for (n_s=0;n_s<nsl; n_s++){ 
     103        // for flat layer 
     104        if (j==1){ 
     105          if (i==n+1) 
     106            break; 
     107          dz = thick[i]; 
     108          sld_i = sld[i]; 
     109          sldim_i = sld_im[i]; 
     110        } 
     111        // for interface 
     112        else{ 
     113          dz = thick_inter[i-1]/nsl; 
     114          if (sld[i-1] == sld[i]){ 
     115            sld_i = sld[i]; 
     116          } 
     117          else{ 
     118            sld_i = intersldfunc(fun_type[i-1],nsl, n_s+0.5, fun_coef[i-1], sld[i-1], sld[i]); 
     119          } 
     120          if (sld_im[i-1] == sld_im[i]){ 
     121            sldim_i = sld_im[i]; 
     122          } 
     123          else{ 
     124            sldim_i = intersldfunc(fun_type[i-1],nsl, n_s+0.5, fun_coef[i-1], sld_im[i-1], sld_im[i]); 
     125          } 
     126        } 
     127        nn = cassign(1.0-sld_i*pow(lamda,2.0)/(2.0*pi),pow(lamda,2.0)/(2.0*pi)*sldim_i); 
     128        nn2=cplx_mult(nn,nn); 
     129 
     130        kn=cplx_sqrt(rcmult(ko2,cplx_sub(nn2,alpha2)));        //nn*ko*sin(phin) 
     131        an=cplx_exp(rcmult(dz,cplx_mult(two,kn))); 
     132 
     133        fnm=cplx_sub(kn,knp1); 
     134        fnp=cplx_add(kn,knp1); 
     135        rn=cplx_div(fnm,fnp); 
     136        Xn=cplx_mult(an,cplx_div(cplx_add(rn,Xnp1),cplx_add(one,cplx_mult(rn,Xnp1))));    //Xn=an*((rn+Xnp1*anp1)/(1+rn*Xnp1*anp1)) 
     137 
     138        Xnp1=Xn; 
     139        knp1=kn; 
     140        // no for-loop for flat layer 
     141        if (j==1) 
     142          break; 
     143      } 
     144    } 
     145  } 
     146  R=pow(Xn.re,2.0)+pow(Xn.im,2.0); 
     147  // This temperarily fixes the total reflection for Rfunction and linear. 
     148  // ToDo: Show why it happens that Xn.re=0 and Xn.im >1! 
     149  if (Xn.im == 0.0 || R > 1){ 
     150    R=1.0; 
     151  } 
     152  R *= scale; 
     153  R += background; 
     154 
     155  free(fun_type); 
     156  free(sld); 
     157  free(sld_im); 
     158  free(thick_inter); 
     159  free(thick); 
     160  free(fun_coef); 
     161 
     162  return R; 
     163 
     164} 
    12165ReflAdvModel :: ReflAdvModel() { 
    13         n_layers = Parameter(1.0); 
    14         scale = Parameter(1.0); 
    15         thick_inter0 = Parameter(1.0); 
    16         func_inter0 = Parameter(0); 
    17         sld_bottom0 = Parameter(2.07e-06); 
    18         sld_medium = Parameter(1.0e-06); 
    19     background = Parameter(0.0); 
    20  
    21  
    22     sld_flat1 = Parameter(2.7e-06); 
    23     sld_flat2 = Parameter(3.5e-06); 
    24     sld_flat3 = Parameter(4.0e-06); 
    25     sld_flat4 = Parameter(3.5e-06); 
    26     sld_flat5 = Parameter(4.0e-06); 
    27     sld_flat6 = Parameter(3.5e-06); 
    28     sld_flat7 = Parameter(4.0e-06); 
    29     sld_flat8 = Parameter(3.5e-06); 
    30     sld_flat9 = Parameter(4.0e-06); 
    31     sld_flat10 = Parameter(3.5e-06); 
    32  
    33  
    34     thick_inter1 = Parameter(1.0); 
    35     thick_inter2 = Parameter(1.0); 
    36     thick_inter3 = Parameter(1.0); 
    37     thick_inter4 = Parameter(1.0); 
    38     thick_inter5 = Parameter(1.0); 
    39     thick_inter6 = Parameter(1.0); 
    40     thick_inter7 = Parameter(1.0); 
    41     thick_inter8 = Parameter(1.0); 
    42     thick_inter9 = Parameter(1.0); 
    43     thick_inter10 = Parameter(1.0); 
    44  
    45  
    46     thick_flat1 = Parameter(15.0); 
    47     thick_flat2 = Parameter(100.0); 
    48     thick_flat3 = Parameter(100.0); 
    49     thick_flat4 = Parameter(100.0); 
    50     thick_flat5 = Parameter(100.0); 
    51     thick_flat6 = Parameter(100.0); 
    52     thick_flat7 = Parameter(100.0); 
    53     thick_flat8 = Parameter(100.0); 
    54     thick_flat9 = Parameter(100.0); 
    55     thick_flat10 = Parameter(100.0); 
    56  
    57  
    58     func_inter1 = Parameter(0); 
    59     func_inter2 = Parameter(0); 
    60     func_inter3 = Parameter(0); 
    61     func_inter4 = Parameter(0); 
    62     func_inter5 = Parameter(0); 
    63     func_inter6 = Parameter(0); 
    64     func_inter7 = Parameter(0); 
    65     func_inter8 = Parameter(0); 
    66     func_inter9 = Parameter(0); 
    67     func_inter10 = Parameter(0); 
    68  
    69     sldIM_flat1 = Parameter(0); 
    70     sldIM_flat2 = Parameter(0); 
    71     sldIM_flat3 = Parameter(0); 
    72     sldIM_flat4 = Parameter(0); 
    73     sldIM_flat5 = Parameter(0); 
    74     sldIM_flat6 = Parameter(0); 
    75     sldIM_flat7 = Parameter(0); 
    76     sldIM_flat8 = Parameter(0); 
    77     sldIM_flat9 = Parameter(0); 
    78     sldIM_flat10 = Parameter(0); 
    79  
    80     nu_inter1 = Parameter(2.5); 
    81     nu_inter2 = Parameter(2.5); 
    82     nu_inter3 = Parameter(2.5); 
    83     nu_inter4 = Parameter(2.5); 
    84     nu_inter5 = Parameter(2.5); 
    85     nu_inter6 = Parameter(2.5); 
    86     nu_inter7 = Parameter(2.5); 
    87     nu_inter8 = Parameter(2.5); 
    88     nu_inter9 = Parameter(2.5); 
    89     nu_inter10 = Parameter(2.5); 
    90  
    91     sldIM_sub0 = Parameter(0.0); 
    92     sldIM_medium = Parameter(0.0); 
    93     npts_inter = Parameter(21.0); 
    94     nu_inter0 = Parameter(2.5); 
     166  n_layers = Parameter(1.0); 
     167  scale = Parameter(1.0); 
     168  thick_inter0 = Parameter(1.0); 
     169  func_inter0 = Parameter(0); 
     170  sld_bottom0 = Parameter(2.07e-06); 
     171  sld_medium = Parameter(1.0e-06); 
     172  background = Parameter(0.0); 
     173 
     174 
     175  sld_flat1 = Parameter(2.7e-06); 
     176  sld_flat2 = Parameter(3.5e-06); 
     177  sld_flat3 = Parameter(4.0e-06); 
     178  sld_flat4 = Parameter(3.5e-06); 
     179  sld_flat5 = Parameter(4.0e-06); 
     180  sld_flat6 = Parameter(3.5e-06); 
     181  sld_flat7 = Parameter(4.0e-06); 
     182  sld_flat8 = Parameter(3.5e-06); 
     183  sld_flat9 = Parameter(4.0e-06); 
     184  sld_flat10 = Parameter(3.5e-06); 
     185 
     186 
     187  thick_inter1 = Parameter(1.0); 
     188  thick_inter2 = Parameter(1.0); 
     189  thick_inter3 = Parameter(1.0); 
     190  thick_inter4 = Parameter(1.0); 
     191  thick_inter5 = Parameter(1.0); 
     192  thick_inter6 = Parameter(1.0); 
     193  thick_inter7 = Parameter(1.0); 
     194  thick_inter8 = Parameter(1.0); 
     195  thick_inter9 = Parameter(1.0); 
     196  thick_inter10 = Parameter(1.0); 
     197 
     198 
     199  thick_flat1 = Parameter(15.0); 
     200  thick_flat2 = Parameter(100.0); 
     201  thick_flat3 = Parameter(100.0); 
     202  thick_flat4 = Parameter(100.0); 
     203  thick_flat5 = Parameter(100.0); 
     204  thick_flat6 = Parameter(100.0); 
     205  thick_flat7 = Parameter(100.0); 
     206  thick_flat8 = Parameter(100.0); 
     207  thick_flat9 = Parameter(100.0); 
     208  thick_flat10 = Parameter(100.0); 
     209 
     210 
     211  func_inter1 = Parameter(0); 
     212  func_inter2 = Parameter(0); 
     213  func_inter3 = Parameter(0); 
     214  func_inter4 = Parameter(0); 
     215  func_inter5 = Parameter(0); 
     216  func_inter6 = Parameter(0); 
     217  func_inter7 = Parameter(0); 
     218  func_inter8 = Parameter(0); 
     219  func_inter9 = Parameter(0); 
     220  func_inter10 = Parameter(0); 
     221 
     222  sldIM_flat1 = Parameter(0); 
     223  sldIM_flat2 = Parameter(0); 
     224  sldIM_flat3 = Parameter(0); 
     225  sldIM_flat4 = Parameter(0); 
     226  sldIM_flat5 = Parameter(0); 
     227  sldIM_flat6 = Parameter(0); 
     228  sldIM_flat7 = Parameter(0); 
     229  sldIM_flat8 = Parameter(0); 
     230  sldIM_flat9 = Parameter(0); 
     231  sldIM_flat10 = Parameter(0); 
     232 
     233  nu_inter1 = Parameter(2.5); 
     234  nu_inter2 = Parameter(2.5); 
     235  nu_inter3 = Parameter(2.5); 
     236  nu_inter4 = Parameter(2.5); 
     237  nu_inter5 = Parameter(2.5); 
     238  nu_inter6 = Parameter(2.5); 
     239  nu_inter7 = Parameter(2.5); 
     240  nu_inter8 = Parameter(2.5); 
     241  nu_inter9 = Parameter(2.5); 
     242  nu_inter10 = Parameter(2.5); 
     243 
     244  sldIM_sub0 = Parameter(0.0); 
     245  sldIM_medium = Parameter(0.0); 
     246  npts_inter = Parameter(21.0); 
     247  nu_inter0 = Parameter(2.5); 
    95248} 
    96249 
    … …  
    101254 */ 
    102255double ReflAdvModel :: operator()(double q) { 
    103         double dp[71]; 
    104         // Fill parameter array for IGOR library 
    105         // Add the background after averaging 
    106         dp[0] = n_layers(); 
    107         dp[1] = scale(); 
    108         dp[2] = thick_inter0(); 
    109         dp[3] = func_inter0(); 
    110         dp[4] = sld_bottom0(); 
    111         dp[5] = sld_medium(); 
    112         dp[6] = background(); 
    113  
    114         dp[7] = sld_flat1(); 
    115         dp[8] = sld_flat2(); 
    116         dp[9] = sld_flat3(); 
    117         dp[10] = sld_flat4(); 
    118         dp[11] = sld_flat5(); 
    119         dp[12] = sld_flat6(); 
    120         dp[13] = sld_flat7(); 
    121         dp[14] = sld_flat8(); 
    122         dp[15] = sld_flat9(); 
    123         dp[16] = sld_flat10(); 
    124  
    125         dp[17] = thick_inter1(); 
    126         dp[18] = thick_inter2(); 
    127         dp[19] = thick_inter3(); 
    128         dp[20] = thick_inter4(); 
    129         dp[21] = thick_inter5(); 
    130         dp[22] = thick_inter6(); 
    131         dp[23] = thick_inter7(); 
    132         dp[24] = thick_inter8(); 
    133         dp[25] = thick_inter9(); 
    134         dp[26] = thick_inter10(); 
    135  
    136         dp[27] = thick_flat1(); 
    137         dp[28] = thick_flat2(); 
    138         dp[29] = thick_flat3(); 
    139         dp[30] = thick_flat4(); 
    140         dp[31] = thick_flat5(); 
    141         dp[32] = thick_flat6(); 
    142         dp[33] = thick_flat7(); 
    143         dp[34] = thick_flat8(); 
    144         dp[35] = thick_flat9(); 
    145         dp[36] = thick_flat10(); 
    146  
    147         dp[37] = func_inter1(); 
    148         dp[38] = func_inter2(); 
    149         dp[39] = func_inter3(); 
    150         dp[40] = func_inter4(); 
    151         dp[41] = func_inter5(); 
    152         dp[42] = func_inter6(); 
    153         dp[43] = func_inter7(); 
    154         dp[44] = func_inter8(); 
    155         dp[45] = func_inter9(); 
    156         dp[46] = func_inter10(); 
    157  
    158         dp[47] = sldIM_flat1(); 
    159         dp[48] = sldIM_flat2(); 
    160         dp[49] = sldIM_flat3(); 
    161         dp[50] = sldIM_flat4(); 
    162         dp[51] = sldIM_flat5(); 
    163         dp[52] = sldIM_flat6(); 
    164         dp[53] = sldIM_flat7(); 
    165         dp[54] = sldIM_flat8(); 
    166         dp[55] = sldIM_flat9(); 
    167         dp[56] = sldIM_flat10(); 
    168  
    169         dp[57] = nu_inter1(); 
    170         dp[58] = nu_inter2(); 
    171         dp[59] = nu_inter3(); 
    172         dp[60] = nu_inter4(); 
    173         dp[61] = nu_inter5(); 
    174         dp[62] = nu_inter6(); 
    175         dp[63] = nu_inter7(); 
    176         dp[64] = nu_inter8(); 
    177         dp[65] = nu_inter9(); 
    178         dp[66] = nu_inter10(); 
    179  
    180         dp[67] = sldIM_sub0(); 
    181         dp[68] = sldIM_medium(); 
    182         dp[69] = npts_inter(); 
    183         dp[70] = nu_inter0(); 
    184         // Get the dispersion points for the radius 
    185         //vector<WeightPoint> weights_thick; 
    186         //thick_inter0.get_weights(weights_thick_inter0); 
    187  
    188  
    189         return re_adv_kernel(dp,q); 
     256  double dp[71]; 
     257  // Fill parameter array for IGOR library 
     258  // Add the background after averaging 
     259  dp[0] = n_layers(); 
     260  dp[1] = scale(); 
     261  dp[2] = thick_inter0(); 
     262  dp[3] = func_inter0(); 
     263  dp[4] = sld_bottom0(); 
     264  dp[5] = sld_medium(); 
     265  dp[6] = background(); 
     266 
     267  dp[7] = sld_flat1(); 
     268  dp[8] = sld_flat2(); 
     269  dp[9] = sld_flat3(); 
     270  dp[10] = sld_flat4(); 
     271  dp[11] = sld_flat5(); 
     272  dp[12] = sld_flat6(); 
     273  dp[13] = sld_flat7(); 
     274  dp[14] = sld_flat8(); 
     275  dp[15] = sld_flat9(); 
     276  dp[16] = sld_flat10(); 
     277 
     278  dp[17] = thick_inter1(); 
     279  dp[18] = thick_inter2(); 
     280  dp[19] = thick_inter3(); 
     281  dp[20] = thick_inter4(); 
     282  dp[21] = thick_inter5(); 
     283  dp[22] = thick_inter6(); 
     284  dp[23] = thick_inter7(); 
     285  dp[24] = thick_inter8(); 
     286  dp[25] = thick_inter9(); 
     287  dp[26] = thick_inter10(); 
     288 
     289  dp[27] = thick_flat1(); 
     290  dp[28] = thick_flat2(); 
     291  dp[29] = thick_flat3(); 
     292  dp[30] = thick_flat4(); 
     293  dp[31] = thick_flat5(); 
     294  dp[32] = thick_flat6(); 
     295  dp[33] = thick_flat7(); 
     296  dp[34] = thick_flat8(); 
     297  dp[35] = thick_flat9(); 
     298  dp[36] = thick_flat10(); 
     299 
     300  dp[37] = func_inter1(); 
     301  dp[38] = func_inter2(); 
     302  dp[39] = func_inter3(); 
     303  dp[40] = func_inter4(); 
     304  dp[41] = func_inter5(); 
     305  dp[42] = func_inter6(); 
     306  dp[43] = func_inter7(); 
     307  dp[44] = func_inter8(); 
     308  dp[45] = func_inter9(); 
     309  dp[46] = func_inter10(); 
     310 
     311  dp[47] = sldIM_flat1(); 
     312  dp[48] = sldIM_flat2(); 
     313  dp[49] = sldIM_flat3(); 
     314  dp[50] = sldIM_flat4(); 
     315  dp[51] = sldIM_flat5(); 
     316  dp[52] = sldIM_flat6(); 
     317  dp[53] = sldIM_flat7(); 
     318  dp[54] = sldIM_flat8(); 
     319  dp[55] = sldIM_flat9(); 
     320  dp[56] = sldIM_flat10(); 
     321 
     322  dp[57] = nu_inter1(); 
     323  dp[58] = nu_inter2(); 
     324  dp[59] = nu_inter3(); 
     325  dp[60] = nu_inter4(); 
     326  dp[61] = nu_inter5(); 
     327  dp[62] = nu_inter6(); 
     328  dp[63] = nu_inter7(); 
     329  dp[64] = nu_inter8(); 
     330  dp[65] = nu_inter9(); 
     331  dp[66] = nu_inter10(); 
     332 
     333  dp[67] = sldIM_sub0(); 
     334  dp[68] = sldIM_medium(); 
     335  dp[69] = npts_inter(); 
     336  dp[70] = nu_inter0(); 
     337  // Get the dispersion points for the radius 
     338  //vector<WeightPoint> weights_thick; 
     339  //thick_inter0.get_weights(weights_thick_inter0); 
     340 
     341 
     342  return re_adv_kernel(dp,q); 
    190343} 
    191344 
    … …  
    197350 */ 
    198351double ReflAdvModel :: operator()(double qx, double qy) { 
    199         // For 2D set qy as q, ignoring qx. 
    200         double q = qy;//sqrt(qx*qx + qy*qy); 
    201         if (q < 0.0){ 
    202                 return 0.0; 
    203         } 
    204         return (*this).operator()(q); 
     352  // For 2D set qy as q, ignoring qx. 
     353  double q = qy;//sqrt(qx*qx + qy*qy); 
     354  if (q < 0.0){ 
     355    return 0.0; 
     356  } 
     357  return (*this).operator()(q); 
    205358} 
    206359 
    … …  
    213366 */ 
    214367double ReflAdvModel :: evaluate_rphi(double q, double phi) { 
    215         return (*this).operator()(q); 
     368  return (*this).operator()(q); 
    216369} 
    217370 
    … …  
    221374 */ 
    222375double ReflAdvModel :: calculate_ER() { 
    223         //NOT implemented yet!!! 
    224         return 0.0; 
    225 } 
     376  //NOT implemented yet!!! 
     377  return 0.0; 
     378} 
Note: See TracChangeset for help on using the changeset viewer.