Changeset 8ff5cb3 in sasview for sansmodels/src/sans


Ignore:
Timestamp:
Dec 14, 2011 3:26:29 PM (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:
96814e1
Parents:
5d75d65a
Message:

Re #4 Fix a few model warnings

Location:
sansmodels/src/sans/models/c_extensions
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • sansmodels/src/sans/models/c_extensions/libmultifunc/librefl.c

    rb9765ad r8ff5cb3  
    2121 
    2222 
    23 complex cadd(x,y) 
     23complex cplx_add(x,y) 
    2424        complex x,y; 
    2525{ 
     
    4040} 
    4141 
    42 complex csub(x,y) 
     42complex cplx_sub(x,y) 
    4343        complex x,y; 
    4444{ 
     
    5050 
    5151 
    52 complex cmult(x,y) 
     52complex cplx_mult(x,y) 
    5353        complex x,y; 
    5454{ 
     
    5959} 
    6060 
    61 complex cdiv(x,y) 
     61complex cplx_div(x,y) 
    6262        complex x,y; 
    6363{ 
     
    6868} 
    6969 
    70 complex cexp(b) 
     70complex cplx_exp(b) 
    7171        complex b; 
    7272{ 
     
    8181 
    8282 
    83 complex csqrt(z)    //see Schaum`s Math Handbook p. 22, 6.6 and 6.10 
     83complex cplx_sqrt(z)    //see Schaum`s Math Handbook p. 22, 6.6 and 6.10 
    8484        complex z; 
    8585{ 
     
    124124} 
    125125 
    126 complex ccos(b) 
     126complex cplx_cos(b) 
    127127        complex b; 
    128128{ 
     
    131131        two = cassign(2.0,0.0); 
    132132        i = cassign(0.0,1.0); 
    133         bi = cmult(b,i); 
    134         negbi = csub(zero,bi); 
    135         z = cdiv(cadd(cexp(bi),cexp(negbi)),two); 
     133        bi = cplx_mult(b,i); 
     134        negbi = cplx_sub(zero,bi); 
     135        z = cplx_div(cplx_add(cplx_exp(bi),cplx_exp(negbi)),two); 
    136136        return z; 
    137137} 
  • sansmodels/src/sans/models/c_extensions/libmultifunc/librefl.h

    ra24f530 r8ff5cb3  
    1616complex cassign(double real, double imag); 
    1717 
    18 complex cadd(complex x,complex y); 
     18complex cplx_add(complex x,complex y); 
    1919 
    2020complex rcmult(double x,complex y); 
    2121 
    22 complex csub(complex x,complex y); 
     22complex cplx_sub(complex x,complex y); 
    2323 
    24 complex cmult(complex x,complex y); 
     24complex cplx_mult(complex x,complex y); 
    2525 
    26 complex cdiv(complex x,complex y); 
     26complex cplx_div(complex x,complex y); 
    2727 
    28 complex cexp(complex b); 
     28complex cplx_exp(complex b); 
    2929 
    30 complex csqrt(complex z); 
     30complex cplx_sqrt(complex z); 
    3131 
    32 complex ccos(complex b); 
     32complex cplx_cos(complex b); 
    3333 
    3434double intersldfunc(int fun_type, double n_sub, double i, double nu, double sld_l, double sld_r); 
  • sansmodels/src/sans/models/c_extensions/refl.c

    r33aea7f r8ff5cb3  
    6969 
    7070    phi = asin(lamda*q/(4.0*pi)); 
    71     phi1 = cdiv(rcmult(phi,one),n_sup); 
    72     alpha = cmult(n_sup,ccos(phi1)); 
    73         alpha2 = cmult(alpha,alpha); 
     71    phi1 = cplx_div(rcmult(phi,one),n_sup); 
     72    alpha = cplx_mult(n_sup,cplx_cos(phi1)); 
     73        alpha2 = cplx_mult(alpha,alpha); 
    7474 
    7575    nnp1=n_sub; 
    76     knp1=csqrt(rcmult(ko2,csub(cmult(nnp1,nnp1),alpha2)));  //nnp1*ko*sin(phinp1) 
     76    knp1=cplx_sqrt(rcmult(ko2,cplx_sub(cplx_mult(nnp1,nnp1),alpha2)));  //nnp1*ko*sin(phinp1) 
    7777    Xnp1=cassign(0.0,0.0); 
    7878    dz = 0.0; 
     
    102102                                } 
    103103                                nn = cassign(1.0-sld_i*pow(lamda,2.0)/(2.0*pi),0.0); 
    104                                 nn2=cmult(nn,nn); 
    105  
    106                                 kn=csqrt(rcmult(ko2,csub(nn2,alpha2)));        //nn*ko*sin(phin) 
    107                                 an=cexp(rcmult(dz,cmult(two,kn))); 
    108  
    109                                 fnm=csub(kn,knp1); 
    110                                 fnp=cadd(kn,knp1); 
    111                                 rn=cdiv(fnm,fnp); 
    112                                 Xn=cmult(an,cdiv(cadd(rn,Xnp1),cadd(one,cmult(rn,Xnp1))));    //Xn=an*((rn+Xnp1*anp1)/(1+rn*Xnp1*anp1)) 
     104                                nn2=cplx_mult(nn,nn); 
     105 
     106                                kn=cplx_sqrt(rcmult(ko2,cplx_sub(nn2,alpha2)));        //nn*ko*sin(phin) 
     107                                an=cplx_exp(rcmult(dz,cplx_mult(two,kn))); 
     108 
     109                                fnm=cplx_sub(kn,knp1); 
     110                                fnp=cplx_add(kn,knp1); 
     111                                rn=cplx_div(fnm,fnp); 
     112                                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)) 
    113113 
    114114                                Xnp1=Xn; 
  • sansmodels/src/sans/models/c_extensions/refl_adv.c

    ra24f530 r8ff5cb3  
    8383 
    8484    phi = asin(lamda*q/(4.0*pi)); 
    85     phi1 = cdiv(rcmult(phi,one),n_sup); 
    86     alpha = cmult(n_sup,ccos(phi1)); 
    87         alpha2 = cmult(alpha,alpha); 
     85    phi1 = cplx_div(rcmult(phi,one),n_sup); 
     86    alpha = cplx_mult(n_sup,cplx_cos(phi1)); 
     87        alpha2 = cplx_mult(alpha,alpha); 
    8888 
    8989    nnp1=n_sub; 
    90     knp1=csqrt(rcmult(ko2,csub(cmult(nnp1,nnp1),alpha2)));  //nnp1*ko*sin(phinp1) 
     90    knp1=cplx_sqrt(rcmult(ko2,cplx_sub(cplx_mult(nnp1,nnp1),alpha2)));  //nnp1*ko*sin(phinp1) 
    9191    Xnp1=cassign(0.0,0.0); 
    9292    dz = 0.0; 
     
    124124                                } 
    125125                                nn = cassign(1.0-sld_i*pow(lamda,2.0)/(2.0*pi),pow(lamda,2.0)/(2.0*pi)*sldim_i); 
    126                                 nn2=cmult(nn,nn); 
    127  
    128                                 kn=csqrt(rcmult(ko2,csub(nn2,alpha2)));        //nn*ko*sin(phin) 
    129                                 an=cexp(rcmult(dz,cmult(two,kn))); 
    130  
    131                                 fnm=csub(kn,knp1); 
    132                                 fnp=cadd(kn,knp1); 
    133                                 rn=cdiv(fnm,fnp); 
    134                                 Xn=cmult(an,cdiv(cadd(rn,Xnp1),cadd(one,cmult(rn,Xnp1))));    //Xn=an*((rn+Xnp1*anp1)/(1+rn*Xnp1*anp1)) 
     126                                nn2=cplx_mult(nn,nn); 
     127 
     128                                kn=cplx_sqrt(rcmult(ko2,cplx_sub(nn2,alpha2)));        //nn*ko*sin(phin) 
     129                                an=cplx_exp(rcmult(dz,cplx_mult(two,kn))); 
     130 
     131                                fnm=cplx_sub(kn,knp1); 
     132                                fnp=cplx_add(kn,knp1); 
     133                                rn=cplx_div(fnm,fnp); 
     134                                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)) 
    135135 
    136136                                Xnp1=Xn; 
Note: See TracChangeset for help on using the changeset viewer.