Changeset e6fa43e in sasview for sansmodels/src/sans/models/c_models


Ignore:
Timestamp:
Aug 3, 2009 6:14:16 PM (15 years ago)
Author:
Jae Cho <jhjcho@…>
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:
0f3fefe
Parents:
ea07075
Message:

more corrections on this model

File:
1 edited

Legend:

Unmodified
Added
Removed
  • sansmodels/src/sans/models/c_models/flexiblecylinder.cpp

    rea07075 re6fa43e  
    6767        length.get_weights(weights_length); 
    6868 
     69        // Get the dispersion points for the kuhn_length 
     70        vector<WeightPoint> weights_kuhn_length; 
     71        kuhn_length.get_weights(weights_kuhn_length); 
     72 
    6973        // Get the dispersion points for the radius 
    7074        vector<WeightPoint> weights_radius; 
     
    8084 
    8185                // Loop over semi axis B weight points 
    82                 for(int j=0; j< (int)weights_radius.size(); j++) { 
    83                         dp[2] = weights_radius[j].value; 
     86                for(int j=0; j< (int)weights_kuhn_length.size(); j++) { 
     87                        dp[2] = weights_kuhn_length[j].value; 
    8488 
    85                         sum += weights_length[i].weight 
    86                                 * weights_radius[j].weight * FlexExclVolCyl(dp, q); 
    87                         norm += weights_length[i].weight 
    88                                         * weights_radius[j].weight; 
     89                        // Loop over semi axis C weight points 
     90                        for(int k=0; k< (int)weights_radius.size(); k++) { 
     91                                dp[3] = weights_radius[k].value; 
     92 
     93                                sum += weights_length[i].weight 
     94                                        * weights_kuhn_length[j].weight * weights_radius[k].weight * FlexExclVolCyl(dp, q); 
     95                                norm += weights_length[i].weight 
     96                                        * weights_kuhn_length[j]*weights_radius[k].weight; 
    8997                } 
    9098        } 
Note: See TracChangeset for help on using the changeset viewer.