Changeset a807206 in sasmodels for sasmodels/models/pearl_necklace.c


Ignore:
Timestamp:
Oct 1, 2016 12:42:06 AM (8 years ago)
Author:
butler
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:
caddb14, 5031ca3
Parents:
2222134
Message:

updating more parameter names addressing #649

File:
1 edited

Legend:

Unmodified
Added
Removed
  • sasmodels/models/pearl_necklace.c

    r2f5c6d4 ra807206  
    1 double _pearl_necklace_kernel(double q, double radius, double edge_separation, 
     1double _pearl_necklace_kernel(double q, double radius, double edge_sep, 
    22        double thick_string, double num_pearls, double sld_pearl, 
    33        double sld_string, double sld_solv); 
    4 double form_volume(double radius, double edge_separation, 
    5         double string_thickness, double number_of_pearls); 
     4double form_volume(double radius, double edge_sep, 
     5        double thick_string, double num_pearls); 
    66 
    7 double Iq(double q, double radius, double edge_separation, 
    8         double string_thickness, double number_of_pearls, double sld,  
     7double Iq(double q, double radius, double edge_sep, 
     8        double thick_string, double num_pearls, double sld,  
    99        double string_sld, double solvent_sld); 
    1010 
    11 #define INVALID(v) (v.string_thickness >= v.radius || v.number_of_pearls <= 0) 
     11#define INVALID(v) (v.thick_string >= v.radius || v.num_pearls <= 0) 
    1212 
    1313// From Igor library 
    14 double _pearl_necklace_kernel(double q, double radius, double edge_separation, double thick_string, 
     14double _pearl_necklace_kernel(double q, double radius, double edge_sep, double thick_string, 
    1515        double num_pearls, double sld_pearl, double sld_string, double sld_solv) 
    1616{ 
     
    2727         
    2828        // center to center distance between the neighboring pearls 
    29         double A_s = edge_separation + 2.0 * radius; 
     29        double A_s = edge_sep + 2.0 * radius; 
    3030         
    3131        // Repeated Calculations 
     
    3333        double oneminussinc = 1 - sincasq; 
    3434        double q_r = q * radius; 
    35         double q_edge = q * edge_separation; 
     35        double q_edge = q * edge_sep; 
    3636         
    3737        // each volume 
    38         double string_vol = edge_separation * pi * thick_string * thick_string / 4.0; 
     38        double string_vol = edge_sep * pi * thick_string * thick_string / 4.0; 
    3939        double pearl_vol = 4.0 / 3.0 * pi * radius * radius * radius; 
    4040 
     
    105105} 
    106106 
    107 double form_volume(double radius, double edge_separation, 
    108         double string_thickness, double number_of_pearls) 
     107double form_volume(double radius, double edge_sep, 
     108        double thick_string, double num_pearls) 
    109109{ 
    110110        double total_vol; 
    111111 
    112112        double pi = 4.0*atan(1.0); 
    113         double number_of_strings = number_of_pearls - 1.0; 
     113        double number_of_strings = num_pearls - 1.0; 
    114114         
    115         double string_vol = edge_separation * pi * string_thickness * string_thickness / 4.0; 
     115        double string_vol = edge_sep * pi * thick_string * thick_string / 4.0; 
    116116        double pearl_vol = 4.0 / 3.0 * pi * radius * radius * radius; 
    117117 
    118118        total_vol = number_of_strings * string_vol; 
    119         total_vol += number_of_pearls * pearl_vol; 
     119        total_vol += num_pearls * pearl_vol; 
    120120 
    121121        return(total_vol); 
    122122} 
    123123 
    124 double Iq(double q, double radius, double edge_separation, 
    125         double string_thickness, double number_of_pearls, double sld,  
     124double Iq(double q, double radius, double edge_sep, 
     125        double thick_string, double num_pearls, double sld,  
    126126        double string_sld, double solvent_sld) 
    127127{ 
    128128        double value, tot_vol; 
    129129         
    130         value = _pearl_necklace_kernel(q, radius, edge_separation, string_thickness, 
    131                 number_of_pearls, sld, string_sld, solvent_sld); 
    132         tot_vol = form_volume(radius, edge_separation, string_thickness, number_of_pearls); 
     130        value = _pearl_necklace_kernel(q, radius, edge_sep, thick_string, 
     131                num_pearls, sld, string_sld, solvent_sld); 
     132        tot_vol = form_volume(radius, edge_sep, thick_string, num_pearls); 
    133133 
    134134        return value*tot_vol; 
Note: See TracChangeset for help on using the changeset viewer.