Changeset 4a82d4d in sasmodels for sasmodels/models/spherical_sld.c


Ignore:
Timestamp:
Feb 17, 2016 11:42:35 AM (8 years ago)
Author:
wojciech
Branches:
master, core_shell_microgels, costrafo411, magnetic_model, release_v0.94, release_v0.95, ticket-1257-vesicle-product, ticket_1156, ticket_1265_superball, ticket_822_more_unit_tests
Children:
2556095
Parents:
7ef3589
Message:

Spherical SLD working expect low q region

File:
1 edited

Legend:

Unmodified
Added
Removed

  • sasmodels/models/spherical_sld.c

    r7ef3589 r4a82d4d  
    22 
    33double Iq(double q, 
    4     int n_shells, double sld_solve, int npts_inter, 
    5     double sld_core_0, double rad_core_0, double thinck_inter_0, int func_inter_0, double nu_inter_0, 
    6     double sld_flat_1, double thick_inter_1, double thick_flat_1, int func_inter_1, double nu_inter_1, 
    7     double sld_flat_2, double thick_inter_2, double thick_flat_2, int func_inter_2, double nu_inter_2, 
    8     double sld_flat_3, double thick_inter_3, double thick_flat_3, int func_inter_3, double nu_inter_3, 
    9     double sld_flat_4, double thick_inter_4, double thick_flat_4, int func_inter_4, double nu_inter_4, 
    10     double sld_flat_5, double thick_inter_5, double thick_flat_5, int func_inter_5, double nu_inter_5, 
    11     double sld_flat_6, double thick_inter_6, double thick_flat_6, int func_inter_6, double nu_inter_6, 
    12     double sld_flat_7, double thick_inter_7, double thick_flat_7, int func_inter_7, double nu_inter_7, 
    13     double sld_flat_8, double thick_inter_8, double thick_flat_8, int func_inter_8, double nu_inter_8, 
    14     double sld_flat_9, double thick_inter_9, double thick_flat_9, int func_inter_9, double nu_inter_9, 
    15     double sld_flat_10, double thick_inter_10, double thick_flat_10, int func_inter_10, double nu_inter_10); 
     4    int n_shells, double thick_inter_0, int func_inter_0, double sld_core_0, double sld_solv, 
     5    double sld_flat_1, double sld_flat_2, double sld_flat_3, double sld_flat_4, double sld_flat_5, 
     6    double sld_flat_6, double sld_flat_7, double sld_flat_8, double sld_flat_9, double sld_flat_10, 
     7    double thick_inter_1, double thick_inter_2, double thick_inter_3, double thick_inter_4, double thick_inter_5, 
     8    double thick_inter_6, double thick_inter_7, double thick_inter_8, double thick_inter_9, double thick_inter_10, 
     9    double thick_flat_1, double thick_flat_2, double thick_flat_3, double thick_flat_4, double thick_flat_5, 
     10    double thick_flat_6, double thick_flat_7, double thick_flat_8, double thick_flat_9, double thick_flat_10, 
     11    int func_inter_1, int func_inter_2, int func_inter_3, int func_inter_4, int func_inter_5, 
     12    int func_inter_6, int func_inter_7, int func_inter_8, int func_inter_9, int func_inter_10, 
     13    double nu_inter_1, double nu_inter_2,double nu_inter_3, double nu_inter_4, double nu_inter_5, 
     14    double nu_inter_6, double nu_inter_7, double nu_inter_8, double nu_inter_9, double nu_inter_10, 
     15    int npts_inter, double nu_inter_0, double rad_core_0); 
    1616 
    1717double Iqxy(double qx, double qy, 
    18     int n_shells, double sld_solve, int npts_inter, 
    19     double sld_core_0, double rad_core_0, double thinck_inter_0, int func_inter_0, double nu_inter_0, 
    20     double sld_flat_1, double thick_inter_1, double thick_flat_1, int func_inter_1, double nu_inter_1, 
    21     double sld_flat_2, double thick_inter_2, double thick_flat_2, int func_inter_2, double nu_inter_2, 
    22     double sld_flat_3, double thick_inter_3, double thick_flat_3, int func_inter_3, double nu_inter_3, 
    23     double sld_flat_4, double thick_inter_4, double thick_flat_4, int func_inter_4, double nu_inter_4, 
    24     double sld_flat_5, double thick_inter_5, double thick_flat_5, int func_inter_5, double nu_inter_5, 
    25     double sld_flat_6, double thick_inter_6, double thick_flat_6, int func_inter_6, double nu_inter_6, 
    26     double sld_flat_7, double thick_inter_7, double thick_flat_7, int func_inter_7, double nu_inter_7, 
    27     double sld_flat_8, double thick_inter_8, double thick_flat_8, int func_inter_8, double nu_inter_8, 
    28     double sld_flat_9, double thick_inter_9, double thick_flat_9, int func_inter_9, double nu_inter_9, 
    29     double sld_flat_10, double thick_inter_10, double thick_flat_10, int func_inter_10, double nu_inter_10); 
     18    int n_shells, double thick_inter_0, int func_inter_0, double sld_core_0, double sld_solv, 
     19    double sld_flat_1, double sld_flat_2, double sld_flat_3, double sld_flat_4, double sld_flat_5, 
     20    double sld_flat_6, double sld_flat_7, double sld_flat_8, double sld_flat_9, double sld_flat_10, 
     21    double thick_inter_1, double thick_inter_2, double thick_inter_3, double thick_inter_4, double thick_inter_5, 
     22    double thick_inter_6, double thick_inter_7, double thick_inter_8, double thick_inter_9, double thick_inter_10, 
     23    double thick_flat_1, double thick_flat_2, double thick_flat_3, double thick_flat_4, double thick_flat_5, 
     24    double thick_flat_6, double thick_flat_7, double thick_flat_8, double thick_flat_9, double thick_flat_10, 
     25    int func_inter_1, int func_inter_2, int func_inter_3, int func_inter_4, int func_inter_5, 
     26    int func_inter_6, int func_inter_7, int func_inter_8, int func_inter_9, int func_inter_10, 
     27    double nu_inter_1, double nu_inter_2,double nu_inter_3, double nu_inter_4, double nu_inter_5, 
     28    double nu_inter_6, double nu_inter_7, double nu_inter_8, double nu_inter_9, double nu_inter_10, 
     29    int npts_inter, double nu_inter_0, double rad_core_0); 
    3030 
    3131//TODO: Check what is for volume for this model 
    … …  
    189189  f2 *= scale; 
    190190  f2 += background; 
    191  
    192191  //free(fun_type); 
    193192  //free(sld); 
    … …  
    206205 */ 
    207206double Iq(double q, 
    208     int n_shells, double sld_solv, int npts_inter, 
    209     double sld_core_0, double rad_core_0, double thick_inter_0, int func_inter_0, double nu_inter_0, 
    210     double sld_flat_1, double thick_inter_1, double thick_flat_1, int func_inter_1, double nu_inter_1, 
    211     double sld_flat_2, double thick_inter_2, double thick_flat_2, int func_inter_2, double nu_inter_2, 
    212     double sld_flat_3, double thick_inter_3, double thick_flat_3, int func_inter_3, double nu_inter_3, 
    213     double sld_flat_4, double thick_inter_4, double thick_flat_4, int func_inter_4, double nu_inter_4, 
    214     double sld_flat_5, double thick_inter_5, double thick_flat_5, int func_inter_5, double nu_inter_5, 
    215     double sld_flat_6, double thick_inter_6, double thick_flat_6, int func_inter_6, double nu_inter_6, 
    216     double sld_flat_7, double thick_inter_7, double thick_flat_7, int func_inter_7, double nu_inter_7, 
    217     double sld_flat_8, double thick_inter_8, double thick_flat_8, int func_inter_8, double nu_inter_8, 
    218     double sld_flat_9, double thick_inter_9, double thick_flat_9, int func_inter_9, double nu_inter_9, 
    219     double sld_flat_10, double thick_inter_10, double thick_flat_10, int func_inter_10, double nu_inter_10) { 
    220  
     207    int n_shells, double thick_inter_0, int func_inter_0, double sld_core_0, double sld_solv, 
     208    double sld_flat_1, double sld_flat_2, double sld_flat_3, double sld_flat_4, double sld_flat_5, 
     209    double sld_flat_6, double sld_flat_7, double sld_flat_8, double sld_flat_9, double sld_flat_10, 
     210    double thick_inter_1, double thick_inter_2, double thick_inter_3, double thick_inter_4, double thick_inter_5, 
     211    double thick_inter_6, double thick_inter_7, double thick_inter_8, double thick_inter_9, double thick_inter_10, 
     212    double thick_flat_1, double thick_flat_2, double thick_flat_3, double thick_flat_4, double thick_flat_5, 
     213    double thick_flat_6, double thick_flat_7, double thick_flat_8, double thick_flat_9, double thick_flat_10, 
     214    int func_inter_1, int func_inter_2, int func_inter_3, int func_inter_4, int func_inter_5, 
     215    int func_inter_6, int func_inter_7, int func_inter_8, int func_inter_9, int func_inter_10, 
     216    double nu_inter_1, double nu_inter_2,double nu_inter_3, double nu_inter_4, double nu_inter_5, 
     217    double nu_inter_6, double nu_inter_7, double nu_inter_8, double nu_inter_9, double nu_inter_10, 
     218    int npts_inter, double nu_inter_0, double rad_core_0) { 
     219 
     220    //printf("Number of points %d\n",npts_inter); 
    221221    double intensity; 
    222222    //TODO: Remove this container at later stage. It is only kept to minimize stupid errors now 
    … …  
    290290    dp[59] = rad_core_0; 
    291291 
    292  
    293292    intensity = sphere_sld_kernel(dp,q); 
    294  
     293    //printf("%10d\n",intensity); 
    295294    return intensity; 
    296295} 
    … …  
    303302 */ 
    304303double Iqxy(double qx, double qy, 
    305     int n_shells, double sld_solve, int npts_inter, 
    306     double sld_core_0, double rad_core_0, double thinck_inter_0, int func_inter_0, double nu_inter_0, 
    307     double sld_flat_1, double thick_inter_1, double thick_flat_1, int func_inter_1, double nu_inter_1, 
    308     double sld_flat_2, double thick_inter_2, double thick_flat_2, int func_inter_2, double nu_inter_2, 
    309     double sld_flat_3, double thick_inter_3, double thick_flat_3, int func_inter_3, double nu_inter_3, 
    310     double sld_flat_4, double thick_inter_4, double thick_flat_4, int func_inter_4, double nu_inter_4, 
    311     double sld_flat_5, double thick_inter_5, double thick_flat_5, int func_inter_5, double nu_inter_5, 
    312     double sld_flat_6, double thick_inter_6, double thick_flat_6, int func_inter_6, double nu_inter_6, 
    313     double sld_flat_7, double thick_inter_7, double thick_flat_7, int func_inter_7, double nu_inter_7, 
    314     double sld_flat_8, double thick_inter_8, double thick_flat_8, int func_inter_8, double nu_inter_8, 
    315     double sld_flat_9, double thick_inter_9, double thick_flat_9, int func_inter_9, double nu_inter_9, 
    316     double sld_flat_10, double thick_inter_10, double thick_flat_10, int func_inter_10, double nu_inter_10) { 
     304    int n_shells, double thick_inter_0, int func_inter_0, double sld_core_0, double sld_solv, 
     305    double sld_flat_1, double sld_flat_2, double sld_flat_3, double sld_flat_4, double sld_flat_5, 
     306    double sld_flat_6, double sld_flat_7, double sld_flat_8, double sld_flat_9, double sld_flat_10, 
     307    double thick_inter_1, double thick_inter_2, double thick_inter_3, double thick_inter_4, double thick_inter_5, 
     308    double thick_inter_6, double thick_inter_7, double thick_inter_8, double thick_inter_9, double thick_inter_10, 
     309    double thick_flat_1, double thick_flat_2, double thick_flat_3, double thick_flat_4, double thick_flat_5, 
     310    double thick_flat_6, double thick_flat_7, double thick_flat_8, double thick_flat_9, double thick_flat_10, 
     311    int func_inter_1, int func_inter_2, int func_inter_3, int func_inter_4, int func_inter_5, 
     312    int func_inter_6, int func_inter_7, int func_inter_8, int func_inter_9, int func_inter_10, 
     313    double nu_inter_1, double nu_inter_2,double nu_inter_3, double nu_inter_4, double nu_inter_5, 
     314    double nu_inter_6, double nu_inter_7, double nu_inter_8, double nu_inter_9, double nu_inter_10, 
     315    int npts_inter, double nu_inter_0, double rad_core_0) { 
    317316 
    318317    double q = sqrt(qx*qx + qy*qy); 
    319     return Iq(q,n_shells, sld_solve, npts_inter, 
    320     sld_core_0, rad_core_0, thinck_inter_0, func_inter_0, nu_inter_0, 
    321     sld_flat_1, thick_inter_1, thick_flat_1, func_inter_1, nu_inter_1, 
    322     sld_flat_2, thick_inter_2, thick_flat_2, func_inter_2, nu_inter_2, 
    323     sld_flat_3, thick_inter_3, thick_flat_3, func_inter_3, nu_inter_3, 
    324     sld_flat_4, thick_inter_4, thick_flat_4, func_inter_4, nu_inter_4, 
    325     sld_flat_5, thick_inter_5, thick_flat_5, func_inter_5, nu_inter_5, 
    326     sld_flat_6, thick_inter_6, thick_flat_6, func_inter_6, nu_inter_6, 
    327     sld_flat_7, thick_inter_7, thick_flat_7, func_inter_7, nu_inter_7, 
    328     sld_flat_8, thick_inter_8, thick_flat_8, func_inter_8, nu_inter_8, 
    329     sld_flat_9, thick_inter_9, thick_flat_9, func_inter_9, nu_inter_9, 
    330     sld_flat_10, thick_inter_10, thick_flat_10, func_inter_10, nu_inter_10); 
     318    return Iq(q, n_shells, thick_inter_0, func_inter_0, sld_core_0, sld_solv, 
     319    sld_flat_1, sld_flat_2, sld_flat_3, sld_flat_4, sld_flat_5, 
     320    sld_flat_6, sld_flat_7, sld_flat_8, sld_flat_9, sld_flat_10, 
     321    thick_inter_1, thick_inter_2, thick_inter_3, thick_inter_4, thick_inter_5, 
     322    thick_inter_6, thick_inter_7, thick_inter_8, thick_inter_9, thick_inter_10, 
     323    thick_flat_1, thick_flat_2, thick_flat_3, thick_flat_4, thick_flat_5, 
     324    thick_flat_6, thick_flat_7, thick_flat_8, thick_flat_9, thick_flat_10, 
     325    func_inter_1, func_inter_2, func_inter_3, func_inter_4, func_inter_5, 
     326    func_inter_6, func_inter_7, func_inter_8, func_inter_9, func_inter_10, 
     327    nu_inter_1, nu_inter_2, nu_inter_3, nu_inter_4, nu_inter_5, 
     328    nu_inter_6, nu_inter_7, nu_inter_8, nu_inter_9, nu_inter_10, 
     329    npts_inter, nu_inter_0, rad_core_0); 
    331330 
    332331} 
Note: See TracChangeset for help on using the changeset viewer.