Changeset 4628e31 in sasview for sansmodels/src/sans/models/c_extensions/bcc.c
- Timestamp:
- Oct 22, 2010 5:50:05 PM (14 years ago)
- 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:
- d7b7156
- Parents:
- 67e258c
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
sansmodels/src/sans/models/c_extensions/bcc.c
r85bb870 r4628e31 74 74 dp[4] = 0.0; 75 75 76 //convert angle degree to radian 77 double pi = 4.0*atan(1.0); 78 double theta = pars->theta * pi/180.0; 79 double phi = pars->phi * pi/180.0; 80 double psi = pars->psi * pi/180.0; 81 76 82 aa = pars->dnn; 77 83 Da = pars->d_factor*aa; … … 85 91 /// instead of against q coordinate(PRB 36(46), 3(6), 1754(3854)) 86 92 // b3 axis orientation 87 b3_x = sin( pars->theta) * cos(pars->phi);//negative sign here???88 b3_y = sin( pars->theta) * sin(pars->phi);89 b3_z = cos( pars->theta);93 b3_x = sin(theta) * cos(phi);//negative sign here??? 94 b3_y = sin(theta) * sin(phi); 95 b3_z = cos(theta); 90 96 cos_val_b3 = b3_x*q_x + b3_y*q_y + b3_z*q_z; 91 97 92 98 alpha = acos(cos_val_b3); 93 99 // b1 axis orientation 94 b1_x = sin(p ars->psi);95 b1_y = cos(p ars->psi);100 b1_x = sin(psi); 101 b1_y = cos(psi); 96 102 cos_val_b1 = (b1_x*q_x + b1_y*q_y); 97 103 // b2 axis orientation
Note: See TracChangeset
for help on using the changeset viewer.