source: sasview/sansmodels/src/igor_wrapper/test/test.c @ df4702f

ESS_GUIESS_GUI_DocsESS_GUI_batch_fittingESS_GUI_bumps_abstractionESS_GUI_iss1116ESS_GUI_iss879ESS_GUI_iss959ESS_GUI_openclESS_GUI_orderingESS_GUI_sync_sascalccostrafo411magnetic_scattrelease-4.1.1release-4.1.2release-4.2.2release_4.0.1ticket-1009ticket-1094-headlessticket-1242-2d-resolutionticket-1243ticket-1249ticket885unittest-saveload
Last change on this file since df4702f was ae3ce4e, checked in by Mathieu Doucet <doucetm@…>, 17 years ago

Moving sansmodels to trunk

  • Property mode set to 100644
File size: 1.6 KB
Line 
1#include "../include/danse.h"
2#include <stdio.h>
3
4int main(){
5        double output;
6        double pars[11];       
7        double q, phi;
8        double theta_vals[17];
9        double theta_w[17];
10        double phi_vals[17];
11        double phi_w[17];
12        int i;
13       
14        pars[0] = 1.0;
15        pars[1] = 20.0; // radius
16        pars[2] = 400.0; // length
17        pars[3] = 1.0;  // contrast
18        pars[4] = 0.0;  // bck
19        pars[5] = 1.57; //theta
20        pars[6] = 0.0;  // phi
21        pars[7] = 0.0;  // disp theta
22        pars[8] = 0.0;  // disp phi
23        pars[9] = 0.0;  // disp radius
24        pars[10] = 0.0;
25       
26        theta_vals[0] = 0.01;
27        theta_vals[1] = 0.1;
28        theta_vals[2] = 0.2;
29        theta_vals[3] = 0.3;
30        theta_vals[4] = 0.4;
31        theta_vals[5] = 0.5;
32        theta_vals[6] = 0.6;
33        theta_vals[7] = 0.7;
34        theta_vals[8] = 0.8;
35        theta_vals[9] = 0.9;
36        theta_vals[10] = 1.0;
37        theta_vals[11] = 1.1;
38        theta_vals[12] = 1.2;
39        theta_vals[13] = 1.3;
40        theta_vals[14] = 1.4;
41        theta_vals[15] = 1.5;
42        theta_vals[16] = 1.6;
43       
44        for(i=0; i<17; i++) {
45           theta_w[i] = 1.0;
46        }
47       
48        q = 0.025;
49        phi = 0.0;
50       
51        // Testing oriented cylinder
52        output = oriented_cylinder_2D(pars, q, phi);
53        printf("oriented output = %g\n", output);
54       
55        // Testing dispersed with no dispersion
56        output = disperse_cylinder_analytical_2D(pars, q, phi);
57        printf("disp output     = %g\n", output);
58
59        // Testing dispersed with dispersion
60        pars[10] = 100.0;
61        pars[7]  = 3.0;
62        output = disperse_cylinder_analytical_2D(pars, q, phi);
63        printf("disp output     = %g\n", output);
64
65        // Testing input dist
66        pars[10] = 0.0;
67        pars[7]  = 0.0;
68        pars[1]  = 20.0;
69        output = cylinder_Weights(pars, phi_vals, phi_w, 0, 
70                                                                theta_vals, theta_w, 17, 
71                                                                q, phi); 
72        printf("weight output     = %g\n", output);
73
74}
Note: See TracBrowser for help on using the repository browser.