Changeset 6cf1cb3 in sasmodels for sasmodels/models/hollow_cylinder.c


Ignore:
Timestamp:
Nov 17, 2015 2:57:11 PM (8 years ago)
Author:
Paul Kienzle <pkienzle@…>
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:
48be770
Parents:
143d596
Message:

hollow cylinder: scale SLD by 1e-6, mark volume and orientation parameters, scale by volume

File:
1 edited

Legend:

Unmodified
Added
Removed
  • sasmodels/models/hollow_cylinder.c

    r84e6942 r6cf1cb3  
    44double form_volume(double radius, double core_radius, double length); 
    55double Iq(double q, double radius, double core_radius, double length, double sld, 
    6         double solvent_sld, double axis_theta, double axis_phi); 
     6        double solvent_sld); 
    77double Iqxy(double qx, double qy, double radius, double core_radius, double length, double sld, 
    8         double solvent_sld, double axis_theta, double axis_phi); 
     8        double solvent_sld, double theta, double phi); 
    99 
    1010// From Igor library 
     
    4848} 
    4949 
    50 //FIXME: Returning values 10^12 times too high 
    5150double Iq(double q, double radius, double core_radius, double length, double sld, 
    52         double solvent_sld, double axis_theta, double axis_phi) 
     51        double solvent_sld) 
    5352{ 
    5453    int i; 
     
    7473        //normalize by volume 
    7574        volume = form_volume(radius, core_radius, length); 
    76         //convert to [cm-1] 
    77         convert = 1.0e8; 
    78         answer = norm*scale*volume*convert; 
     75        //convert to [cm-1] given sld*1e6 
     76        convert = 1.0e-4; 
     77        answer = norm*scale*convert*volume*volume; 
    7978         
    8079        return(answer); 
     
    8382//TODO: Add this in 
    8483double Iqxy(double qx, double qy, double radius, double core_radius, double length, double sld, 
    85         double solvent_sld, double axis_theta, double axis_phi) 
     84        double solvent_sld, double theta, double phi) 
    8685{ 
    8786    return(0.0); 
Note: See TracChangeset for help on using the changeset viewer.