Changeset 144e032a in sasview for src/sas/sascalc/calculator/c_extensions/sld2i.c

Timestamp:

Nov 14, 2017 8:57:07 AM (6 years ago)

Author:

Paul Kienzle <pkienzle@…>

Branches:

master, magnetic_scatt, release-4.2.2, ticket-1009, ticket-1094-headless, ticket-1242-2d-resolution, ticket-1243, ticket-1249, ticket885, unittest-saveload

Children:

f0ce6e2

Parents:

a57ae07

Message:

tinycc doesn't return structures, so must pass return structure as pointer

File:

• src/sas/sascalc/calculator/c_extensions/sld2i.c (modified) (3 diffs)

src/sas/sascalc/calculator/c_extensions/sld2i.c

@@ -54 +54 @@
         polar_sld b_sld;
         double qr = 0.0;
-        complex iqr = cassign(0.0, 0.0);
-        complex ephase = cassign(0.0, 0.0);
-        complex comp_sld = cassign(0.0, 0.0);
-
-        complex sumj_uu;
-        complex sumj_ud;
-        complex sumj_du;
-        complex sumj_dd;
-        complex temp_fi;
+        Cplx iqr;
+        Cplx ephase;
+        Cplx comp_sld;
+
+        Cplx sumj_uu;
+        Cplx sumj_ud;
+        Cplx sumj_du;
+        Cplx sumj_dd;
+        Cplx temp_fi;
 
         double count = 0.0;
         int i, j;
+
+        cassign(&iqr, 0.0, 0.0);
+        cassign(&ephase, 0.0, 0.0);
+        cassign(&comp_sld, 0.0, 0.0);
 
         //Assume that pixel volumes are given in vol_pix in A^3 unit
@@ -77 +81 @@
         for(i=0; i<npoints; i++){
                 //I_out[i] = 0.0;
-                sumj_uu = cassign(0.0, 0.0);
-                sumj_ud = cassign(0.0, 0.0);
-                sumj_du = cassign(0.0, 0.0);
-                sumj_dd = cassign(0.0, 0.0);
+                cassign(&sumj_uu, 0.0, 0.0);
+                cassign(&sumj_ud, 0.0, 0.0);
+                cassign(&sumj_du, 0.0, 0.0);
+                cassign(&sumj_dd, 0.0, 0.0);
                 //printf("i: %d\n", i);
                 //q = sqrt(qx[i]*qx[i] + qy[i]*qy[i]); // + qz[i]*qz[i]);
 
                 for(j=0; j<this->n_pix; j++){
-                        //printf("j: %d\n", j);
                         if (this->sldn_val[j]!=0.0
                                 ||this->mx_val[j]!=0.0
@@ -91 +94 @@
                                 ||this->mz_val[j]!=0.0)
                         {
+                            // printf("i,j: %d,%d\n", i,j);
                                 //anisotropic
-                                temp_fi = cassign(0.0, 0.0);
-                                b_sld = cal_msld(0, qx[i], qy[i], this->sldn_val[j],
+                                cassign(&temp_fi, 0.0, 0.0);
+                                cal_msld(&b_sld, 0, qx[i], qy[i], this->sldn_val[j],
                                                          this->mx_val[j], this->my_val[j], this->mz_val[j],
                                                          this->inspin, this->outspin, this->stheta);
                                 qr = (qx[i]*this->x_val[j] + qy[i]*this->y_val[j]);
-                                iqr = cassign(0.0, qr);
-                                ephase = cplx_exp(iqr);
+                                cassign(&iqr, 0.0, qr);
+                                cplx_exp(&ephase, iqr);
 
                                 //Let's multiply pixel(atomic) volume here
-                                ephase = rcmult(this->vol_pix[j], ephase);
+                                rcmult(&ephase, this->vol_pix[j], ephase);
                                 //up_up
                                 if (this->inspin > 0.0 && this->outspin > 0.0){
-                                        comp_sld = cassign(b_sld.uu, 0.0);
-                                        temp_fi = cplx_mult(comp_sld, ephase);
-                                        sumj_uu = cplx_add(sumj_uu, temp_fi);
+                                        cassign(&comp_sld, b_sld.uu, 0.0);
+                                        cplx_mult(&temp_fi, comp_sld, ephase);
+                                        cplx_add(&sumj_uu, sumj_uu, temp_fi);
                                 }
                                 //down_down
                                 if (this->inspin < 1.0 && this->outspin < 1.0){
-                                        comp_sld = cassign(b_sld.dd, 0.0);
-                                        temp_fi = cplx_mult(comp_sld, ephase);
-                                        sumj_dd = cplx_add(sumj_dd, temp_fi);
+                                        cassign(&comp_sld, b_sld.dd, 0.0);
+                                        cplx_mult(&temp_fi, comp_sld, ephase);
+                                        cplx_add(&sumj_dd, sumj_dd, temp_fi);
                                 }
                                 //up_down
                                 if (this->inspin > 0.0 && this->outspin < 1.0){
-                                        comp_sld = cassign(b_sld.re_ud, b_sld.im_ud);
-                                        temp_fi = cplx_mult(comp_sld, ephase);
-                                        sumj_ud = cplx_add(sumj_ud, temp_fi);
+                                        cassign(&comp_sld, b_sld.re_ud, b_sld.im_ud);
+                                        cplx_mult(&temp_fi, comp_sld, ephase);
+                                        cplx_add(&sumj_ud, sumj_ud, temp_fi);
                                 }
                                 //down_up
                                 if (this->inspin < 1.0 && this->outspin > 0.0){
-                                        comp_sld = cassign(b_sld.re_du, b_sld.im_du);
-                                        temp_fi = cplx_mult(comp_sld, ephase);
-                                        sumj_du = cplx_add(sumj_du, temp_fi);
-                                }
-
+                                        cassign(&comp_sld, b_sld.re_du, b_sld.im_du);
+                                        cplx_mult(&temp_fi, comp_sld, ephase);
+                                        cplx_add(&sumj_du, sumj_du, temp_fi);
+                                }
 
                                 if (i == 0){