Changeset c036ddb in sasmodels for sasmodels/models/ellipsoid.c


Ignore:
Timestamp:
Aug 7, 2018 8:45:45 PM (6 years ago)
Author:
Paul Kienzle <pkienzle@…>
Branches:
master, core_shell_microgels, magnetic_model, ticket-1257-vesicle-product, ticket_1156, ticket_1265_superball, ticket_822_more_unit_tests
Children:
7e923c2
Parents:
7b0abf8
Message:

refactor so Iq is not needed if Fq is defined

File:
1 edited

Legend:

Unmodified
Added
Removed
  • sasmodels/models/ellipsoid.c

    r01c8d9e rc036ddb  
    55} 
    66 
     7/* Fq overrides Iq 
    78static  double 
    89Iq(double q, 
     
    2021    //     i(h) = int_0^1 Phi^2(h a sqrt(1 + u^2(v^2-1)) du 
    2122    const double v_square_minus_one = square(radius_polar/radius_equatorial) - 1.0; 
    22    
     23 
    2324    // translate a point in [-1,1] to a point in [0, 1] 
    2425    // const double u = GAUSS_Z[i]*(upper-lower)/2 + (upper+lower)/2; 
     
    3637    const double s = (sld - sld_solvent) * form_volume(radius_polar, radius_equatorial); 
    3738    return 1.0e-4 * s * s * form; 
    38 }  
     39} 
     40*/ 
    3941 
    4042static void 
     
    7274    const double form_avg = total_F1*zm; 
    7375    const double s = (sld - sld_solvent) * form_volume(radius_polar, radius_equatorial); 
     76    *F1 = 1e-2 * s * form_avg; 
    7477    *F2 = 1e-4 * s * s * form_squared_avg; 
    75     *F1 = 1e-2 * s * form_avg; 
    7678} 
    77  
    78  
    79  
    80  
    81  
    8279 
    8380static double 
Note: See TracChangeset for help on using the changeset viewer.