Changeset 0b082f3 in sasview
- Timestamp:
- Dec 8, 2011 4:50:20 PM (13 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:
- 01de557
- Parents:
- 0438933
- Location:
- sansmodels
- Files:
-
- 52 edited
Legend:
- Unmodified
- Added
- Removed
-
sansmodels/setup.py
rf9303b6 r0b082f3 5 5 import sys 6 6 import os 7 import platform 7 8 from numpy.distutils.misc_util import get_numpy_include_dirs 8 9 numpy_incl_path = os.path.join(get_numpy_include_dirs()[0], "numpy") … … 69 70 smearer_sources.append(os.path.join(igordir, "winFuncs.c")) 70 71 72 # Enable OpenMP 73 extra_compile_args = [] 74 extra_link_args = [] 75 if sys.platform=='linux2' or (sys.platform=='darwin' and platform.architecture()[0]=='64bit'): 76 extra_compile_args = ['-fopenmp'] 77 extra_link_args = ['-lgomp'] 78 elif os.name=='nt': 79 extra_compile_args = ['/openmp'] 80 71 81 dist = setup( 72 82 name="sansmodels", … … 90 100 91 101 ext_modules = [ Extension("sans.models.sans_extension.c_models", 92 sources=model_sources,102 sources=model_sources, 93 103 94 include_dirs=[igordir, srcdir, c_model_dir, numpy_incl_path]), 104 include_dirs=[igordir, srcdir, c_model_dir, numpy_incl_path], 105 extra_compile_args=extra_compile_args, 106 extra_link_args=extra_link_args 107 ), 108 95 109 # Smearer extension 96 110 Extension("sans.models.sans_extension.smearer", … … 101 115 "smearer2d_helper_module.cpp"), 102 116 os.path.join(smear_dir, "smearer2d_helper.cpp"),], 103 include_dirs=[smear_dir,numpy_incl_path] 117 include_dirs=[smear_dir,numpy_incl_path], 118 extra_compile_args=extra_compile_args, 119 extra_link_args=extra_link_args 104 120 ) 105 121 ] -
sansmodels/src/sans/models/ReflAdvModel.py
rd760c0a r0b082f3 2 2 3 3 ############################################################################## 4 # 5 # 6 # 4 # This software was developed by the University of Tennessee as part of the 5 # Distributed Data Analysis of Neutron Scattering Experiments (DANSE) 6 # project funded by the US National Science Foundation. 7 7 # 8 # 9 # 10 # 8 # If you use DANSE applications to do scientific research that leads to 9 # publication, we ask that you acknowledge the use of the software with the 10 # following sentence: 11 11 # 12 # 12 # "This work benefited from DANSE software developed under NSF award DMR-0520547." 13 13 # 14 # 14 # copyright 2008, University of Tennessee 15 15 ############################################################################## 16 16 … … 127 127 ## Model description 128 128 self.description ="""Calculate neutron reflectivity using the Parratt iterative formula 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 129 Parameters: 130 background:background 131 scale: scale factor 132 sld_bottom0: the SLD of the substrate 133 sld_medium: the SLD of the incident medium 134 or superstrate 135 sld_flatN: the SLD of the flat region of 136 the N'th layer 137 thick_flatN: the thickness of the flat 138 region of the N'th layer 139 func_interN: the function used to describe 140 the interface of the N'th layer 141 nu_interN: the coefficient for the func_interN 142 thick_interN: the thickness of the interface 143 of the N'th layer 144 Note: the layer number starts to increase 145 from the bottom (substrate) to the top.""" 146 146 147 147 ## Parameter details [units, min, max] … … 245 245 """ Return a identical copy of self """ 246 246 return self._clone(ReflAdvModel()) 247 247 248 248 249 249 def run(self, x=0.0): -
sansmodels/src/sans/models/ReflModel.py
r500be82 r0b082f3 2 2 3 3 ############################################################################## 4 # 5 # 6 # 4 # This software was developed by the University of Tennessee as part of the 5 # Distributed Data Analysis of Neutron Scattering Experiments (DANSE) 6 # project funded by the US National Science Foundation. 7 7 # 8 # 9 # 10 # 8 # If you use DANSE applications to do scientific research that leads to 9 # publication, we ask that you acknowledge the use of the software with the 10 # following sentence: 11 11 # 12 # 12 # "This work benefited from DANSE software developed under NSF award DMR-0520547." 13 13 # 14 # 14 # copyright 2008, University of Tennessee 15 15 ############################################################################## 16 16 … … 103 103 ## Model description 104 104 self.description ="""Calculate neutron reflectivity using the Parratt iterative formula 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 105 Parameters: 106 background:background 107 scale: scale factor 108 sld_bottom0: the SLD of the substrate 109 sld_medium: the SLD of the incident medium 110 or superstrate 111 sld_flatN: the SLD of the flat region of 112 the N'th layer 113 thick_flatN: the thickness of the flat 114 region of the N'th layer 115 func_interN: the function used to describe 116 the interface of the N'th layer 117 thick_interN: the thickness of the interface 118 of the N'th layer 119 Note: the layer number starts to increase 120 from the bottom (substrate) to the top.""" 121 121 122 122 ## Parameter details [units, min, max] … … 196 196 """ Return a identical copy of self """ 197 197 return self._clone(ReflModel()) 198 198 199 199 200 200 def run(self, x=0.0): -
sansmodels/src/sans/models/c_models/CBCCrystalModel.cpp
r2605da22 r0b082f3 179 179 return NULL; 180 180 } 181 #pragma omp parallel for 181 182 for (int i = 0; i < q->dimensions[0]; i++){ 182 183 double q_value = *(double *)(q->data + i*q->strides[0]); … … 196 197 { 197 198 PyArrayObject *result; 198 int i,x_len, y_len, dims[1];199 int x_len, y_len, dims[1]; 199 200 //check validity of input vectors 200 201 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 220 221 221 222 /* Do the calculation. */ 222 for ( i=0; i< x_len; i++) { 223 #pragma omp parallel for 224 for (int i=0; i< x_len; i++) { 223 225 double x_value = *(double *)(x->data + i*x->strides[0]); 224 226 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CBarBellModel.cpp
r2605da22 r0b082f3 181 181 return NULL; 182 182 } 183 #pragma omp parallel for 183 184 for (int i = 0; i < q->dimensions[0]; i++){ 184 185 double q_value = *(double *)(q->data + i*q->strides[0]); … … 198 199 { 199 200 PyArrayObject *result; 200 int i,x_len, y_len, dims[1];201 int x_len, y_len, dims[1]; 201 202 //check validity of input vectors 202 203 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 222 223 223 224 /* Do the calculation. */ 224 for ( i=0; i< x_len; i++) { 225 #pragma omp parallel for 226 for (int i=0; i< x_len; i++) { 225 227 double x_value = *(double *)(x->data + i*x->strides[0]); 226 228 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CBinaryHSModel.cpp
r2605da22 r0b082f3 171 171 return NULL; 172 172 } 173 #pragma omp parallel for 173 174 for (int i = 0; i < q->dimensions[0]; i++){ 174 175 double q_value = *(double *)(q->data + i*q->strides[0]); … … 188 189 { 189 190 PyArrayObject *result; 190 int i,x_len, y_len, dims[1];191 int x_len, y_len, dims[1]; 191 192 //check validity of input vectors 192 193 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 212 213 213 214 /* Do the calculation. */ 214 for ( i=0; i< x_len; i++) { 215 #pragma omp parallel for 216 for (int i=0; i< x_len; i++) { 215 217 double x_value = *(double *)(x->data + i*x->strides[0]); 216 218 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CCSParallelepipedModel.cpp
r2605da22 r0b082f3 191 191 return NULL; 192 192 } 193 #pragma omp parallel for 193 194 for (int i = 0; i < q->dimensions[0]; i++){ 194 195 double q_value = *(double *)(q->data + i*q->strides[0]); … … 208 209 { 209 210 PyArrayObject *result; 210 int i,x_len, y_len, dims[1];211 int x_len, y_len, dims[1]; 211 212 //check validity of input vectors 212 213 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 232 233 233 234 /* Do the calculation. */ 234 for ( i=0; i< x_len; i++) { 235 #pragma omp parallel for 236 for (int i=0; i< x_len; i++) { 235 237 double x_value = *(double *)(x->data + i*x->strides[0]); 236 238 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CCappedCylinderModel.cpp
r2605da22 r0b082f3 181 181 return NULL; 182 182 } 183 #pragma omp parallel for 183 184 for (int i = 0; i < q->dimensions[0]; i++){ 184 185 double q_value = *(double *)(q->data + i*q->strides[0]); … … 198 199 { 199 200 PyArrayObject *result; 200 int i,x_len, y_len, dims[1];201 int x_len, y_len, dims[1]; 201 202 //check validity of input vectors 202 203 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 222 223 223 224 /* Do the calculation. */ 224 for ( i=0; i< x_len; i++) { 225 #pragma omp parallel for 226 for (int i=0; i< x_len; i++) { 225 227 double x_value = *(double *)(x->data + i*x->strides[0]); 226 228 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CCoreFourShellModel.cpp
r2605da22 r0b082f3 185 185 return NULL; 186 186 } 187 #pragma omp parallel for 187 188 for (int i = 0; i < q->dimensions[0]; i++){ 188 189 double q_value = *(double *)(q->data + i*q->strides[0]); … … 202 203 { 203 204 PyArrayObject *result; 204 int i,x_len, y_len, dims[1];205 int x_len, y_len, dims[1]; 205 206 //check validity of input vectors 206 207 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 226 227 227 228 /* Do the calculation. */ 228 for ( i=0; i< x_len; i++) { 229 #pragma omp parallel for 230 for (int i=0; i< x_len; i++) { 229 231 double x_value = *(double *)(x->data + i*x->strides[0]); 230 232 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CCoreShellCylinderModel.cpp
r2605da22 r0b082f3 182 182 return NULL; 183 183 } 184 #pragma omp parallel for 184 185 for (int i = 0; i < q->dimensions[0]; i++){ 185 186 double q_value = *(double *)(q->data + i*q->strides[0]); … … 199 200 { 200 201 PyArrayObject *result; 201 int i,x_len, y_len, dims[1];202 int x_len, y_len, dims[1]; 202 203 //check validity of input vectors 203 204 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 223 224 224 225 /* Do the calculation. */ 225 for ( i=0; i< x_len; i++) { 226 #pragma omp parallel for 227 for (int i=0; i< x_len; i++) { 226 228 double x_value = *(double *)(x->data + i*x->strides[0]); 227 229 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CCoreShellEllipsoidModel.cpp
r2605da22 r0b082f3 186 186 return NULL; 187 187 } 188 #pragma omp parallel for 188 189 for (int i = 0; i < q->dimensions[0]; i++){ 189 190 double q_value = *(double *)(q->data + i*q->strides[0]); … … 203 204 { 204 205 PyArrayObject *result; 205 int i,x_len, y_len, dims[1];206 int x_len, y_len, dims[1]; 206 207 //check validity of input vectors 207 208 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 227 228 228 229 /* Do the calculation. */ 229 for ( i=0; i< x_len; i++) { 230 #pragma omp parallel for 231 for (int i=0; i< x_len; i++) { 230 232 double x_value = *(double *)(x->data + i*x->strides[0]); 231 233 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CCoreShellModel.cpp
r2605da22 r0b082f3 170 170 return NULL; 171 171 } 172 #pragma omp parallel for 172 173 for (int i = 0; i < q->dimensions[0]; i++){ 173 174 double q_value = *(double *)(q->data + i*q->strides[0]); … … 187 188 { 188 189 PyArrayObject *result; 189 int i,x_len, y_len, dims[1];190 int x_len, y_len, dims[1]; 190 191 //check validity of input vectors 191 192 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 211 212 212 213 /* Do the calculation. */ 213 for ( i=0; i< x_len; i++) { 214 #pragma omp parallel for 215 for (int i=0; i< x_len; i++) { 214 216 double x_value = *(double *)(x->data + i*x->strides[0]); 215 217 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CCylinderModel.cpp
r1a191e4 r0b082f3 195 195 { 196 196 PyArrayObject *result; 197 int i,x_len, y_len, dims[1];197 int x_len, y_len, dims[1]; 198 198 //check validity of input vectors 199 199 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 219 219 220 220 /* Do the calculation. */ 221 for ( i=0; i< x_len; i++) { 221 #pragma omp parallel for 222 for (int i=0; i< x_len; i++) { 222 223 double x_value = *(double *)(x->data + i*x->strides[0]); 223 224 double y_value = *(double *)(y->data + i*y->strides[0]); … … 284 285 if (PyArray_Check(pars)) { 285 286 return evaluateOneDim(self->model, (PyArrayObject*)pars); 286 287 287 } 288 }else{ 288 289 PyErr_SetString(CCylinderModelError, 289 290 "CCylinderModel.evalDistribution expect numpy array of one dimension."); 290 291 return NULL; 291 292 } 292 293 }else if( PyList_Check(pars)==1) { 293 294 // Length of list should be 2 for I(qx,qy) -
sansmodels/src/sans/models/c_models/CDiamCylFunc.cpp
r2605da22 r0b082f3 165 165 return NULL; 166 166 } 167 #pragma omp parallel for 167 168 for (int i = 0; i < q->dimensions[0]; i++){ 168 169 double q_value = *(double *)(q->data + i*q->strides[0]); … … 182 183 { 183 184 PyArrayObject *result; 184 int i,x_len, y_len, dims[1];185 int x_len, y_len, dims[1]; 185 186 //check validity of input vectors 186 187 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 206 207 207 208 /* Do the calculation. */ 208 for ( i=0; i< x_len; i++) { 209 #pragma omp parallel for 210 for (int i=0; i< x_len; i++) { 209 211 double x_value = *(double *)(x->data + i*x->strides[0]); 210 212 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CDiamEllipFunc.cpp
r2605da22 r0b082f3 165 165 return NULL; 166 166 } 167 #pragma omp parallel for 167 168 for (int i = 0; i < q->dimensions[0]; i++){ 168 169 double q_value = *(double *)(q->data + i*q->strides[0]); … … 182 183 { 183 184 PyArrayObject *result; 184 int i,x_len, y_len, dims[1];185 int x_len, y_len, dims[1]; 185 186 //check validity of input vectors 186 187 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 206 207 207 208 /* Do the calculation. */ 208 for ( i=0; i< x_len; i++) { 209 #pragma omp parallel for 210 for (int i=0; i< x_len; i++) { 209 211 double x_value = *(double *)(x->data + i*x->strides[0]); 210 212 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CEllipsoidModel.cpp
r2605da22 r0b082f3 177 177 return NULL; 178 178 } 179 #pragma omp parallel for 179 180 for (int i = 0; i < q->dimensions[0]; i++){ 180 181 double q_value = *(double *)(q->data + i*q->strides[0]); … … 194 195 { 195 196 PyArrayObject *result; 196 int i,x_len, y_len, dims[1];197 int x_len, y_len, dims[1]; 197 198 //check validity of input vectors 198 199 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 218 219 219 220 /* Do the calculation. */ 220 for ( i=0; i< x_len; i++) { 221 #pragma omp parallel for 222 for (int i=0; i< x_len; i++) { 221 223 double x_value = *(double *)(x->data + i*x->strides[0]); 222 224 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CEllipticalCylinderModel.cpp
r2605da22 r0b082f3 185 185 return NULL; 186 186 } 187 #pragma omp parallel for 187 188 for (int i = 0; i < q->dimensions[0]; i++){ 188 189 double q_value = *(double *)(q->data + i*q->strides[0]); … … 202 203 { 203 204 PyArrayObject *result; 204 int i,x_len, y_len, dims[1];205 int x_len, y_len, dims[1]; 205 206 //check validity of input vectors 206 207 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 226 227 227 228 /* Do the calculation. */ 228 for ( i=0; i< x_len; i++) { 229 #pragma omp parallel for 230 for (int i=0; i< x_len; i++) { 229 231 double x_value = *(double *)(x->data + i*x->strides[0]); 230 232 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CFCCrystalModel.cpp
r2605da22 r0b082f3 179 179 return NULL; 180 180 } 181 #pragma omp parallel for 181 182 for (int i = 0; i < q->dimensions[0]; i++){ 182 183 double q_value = *(double *)(q->data + i*q->strides[0]); … … 196 197 { 197 198 PyArrayObject *result; 198 int i,x_len, y_len, dims[1];199 int x_len, y_len, dims[1]; 199 200 //check validity of input vectors 200 201 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 220 221 221 222 /* Do the calculation. */ 222 for ( i=0; i< x_len; i++) { 223 #pragma omp parallel for 224 for (int i=0; i< x_len; i++) { 223 225 double x_value = *(double *)(x->data + i*x->strides[0]); 224 226 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CFlexCylEllipXModel.cpp
r2605da22 r0b082f3 177 177 return NULL; 178 178 } 179 #pragma omp parallel for 179 180 for (int i = 0; i < q->dimensions[0]; i++){ 180 181 double q_value = *(double *)(q->data + i*q->strides[0]); … … 194 195 { 195 196 PyArrayObject *result; 196 int i,x_len, y_len, dims[1];197 int x_len, y_len, dims[1]; 197 198 //check validity of input vectors 198 199 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 218 219 219 220 /* Do the calculation. */ 220 for ( i=0; i< x_len; i++) { 221 #pragma omp parallel for 222 for (int i=0; i< x_len; i++) { 221 223 double x_value = *(double *)(x->data + i*x->strides[0]); 222 224 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CFlexibleCylinderModel.cpp
r2605da22 r0b082f3 173 173 return NULL; 174 174 } 175 #pragma omp parallel for 175 176 for (int i = 0; i < q->dimensions[0]; i++){ 176 177 double q_value = *(double *)(q->data + i*q->strides[0]); … … 190 191 { 191 192 PyArrayObject *result; 192 int i,x_len, y_len, dims[1];193 int x_len, y_len, dims[1]; 193 194 //check validity of input vectors 194 195 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 214 215 215 216 /* Do the calculation. */ 216 for ( i=0; i< x_len; i++) { 217 #pragma omp parallel for 218 for (int i=0; i< x_len; i++) { 217 219 double x_value = *(double *)(x->data + i*x->strides[0]); 218 220 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CFractalModel.cpp
r2605da22 r0b082f3 167 167 return NULL; 168 168 } 169 #pragma omp parallel for 169 170 for (int i = 0; i < q->dimensions[0]; i++){ 170 171 double q_value = *(double *)(q->data + i*q->strides[0]); … … 184 185 { 185 186 PyArrayObject *result; 186 int i,x_len, y_len, dims[1];187 int x_len, y_len, dims[1]; 187 188 //check validity of input vectors 188 189 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 208 209 209 210 /* Do the calculation. */ 210 for ( i=0; i< x_len; i++) { 211 #pragma omp parallel for 212 for (int i=0; i< x_len; i++) { 211 213 double x_value = *(double *)(x->data + i*x->strides[0]); 212 214 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CFuzzySphereModel.cpp
r2605da22 r0b082f3 169 169 return NULL; 170 170 } 171 #pragma omp parallel for 171 172 for (int i = 0; i < q->dimensions[0]; i++){ 172 173 double q_value = *(double *)(q->data + i*q->strides[0]); … … 186 187 { 187 188 PyArrayObject *result; 188 int i,x_len, y_len, dims[1];189 int x_len, y_len, dims[1]; 189 190 //check validity of input vectors 190 191 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 210 211 211 212 /* Do the calculation. */ 212 for ( i=0; i< x_len; i++) { 213 #pragma omp parallel for 214 for (int i=0; i< x_len; i++) { 213 215 double x_value = *(double *)(x->data + i*x->strides[0]); 214 216 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CGaussian.cpp
r2605da22 r0b082f3 160 160 return NULL; 161 161 } 162 #pragma omp parallel for 162 163 for (int i = 0; i < q->dimensions[0]; i++){ 163 164 double q_value = *(double *)(q->data + i*q->strides[0]); … … 177 178 { 178 179 PyArrayObject *result; 179 int i,x_len, y_len, dims[1];180 int x_len, y_len, dims[1]; 180 181 //check validity of input vectors 181 182 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 201 202 202 203 /* Do the calculation. */ 203 for ( i=0; i< x_len; i++) { 204 #pragma omp parallel for 205 for (int i=0; i< x_len; i++) { 204 206 double x_value = *(double *)(x->data + i*x->strides[0]); 205 207 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CHardsphereStructure.cpp
r2605da22 r0b082f3 162 162 return NULL; 163 163 } 164 #pragma omp parallel for 164 165 for (int i = 0; i < q->dimensions[0]; i++){ 165 166 double q_value = *(double *)(q->data + i*q->strides[0]); … … 179 180 { 180 181 PyArrayObject *result; 181 int i,x_len, y_len, dims[1];182 int x_len, y_len, dims[1]; 182 183 //check validity of input vectors 183 184 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 203 204 204 205 /* Do the calculation. */ 205 for ( i=0; i< x_len; i++) { 206 #pragma omp parallel for 207 for (int i=0; i< x_len; i++) { 206 208 double x_value = *(double *)(x->data + i*x->strides[0]); 207 209 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CHayterMSAStructure.cpp
r2605da22 r0b082f3 166 166 return NULL; 167 167 } 168 #pragma omp parallel for 168 169 for (int i = 0; i < q->dimensions[0]; i++){ 169 170 double q_value = *(double *)(q->data + i*q->strides[0]); … … 183 184 { 184 185 PyArrayObject *result; 185 int i,x_len, y_len, dims[1];186 int x_len, y_len, dims[1]; 186 187 //check validity of input vectors 187 188 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 207 208 208 209 /* Do the calculation. */ 209 for ( i=0; i< x_len; i++) { 210 #pragma omp parallel for 211 for (int i=0; i< x_len; i++) { 210 212 double x_value = *(double *)(x->data + i*x->strides[0]); 211 213 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CHollowCylinderModel.cpp
r2605da22 r0b082f3 181 181 return NULL; 182 182 } 183 #pragma omp parallel for 183 184 for (int i = 0; i < q->dimensions[0]; i++){ 184 185 double q_value = *(double *)(q->data + i*q->strides[0]); … … 198 199 { 199 200 PyArrayObject *result; 200 int i,x_len, y_len, dims[1];201 int x_len, y_len, dims[1]; 201 202 //check validity of input vectors 202 203 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 222 223 223 224 /* Do the calculation. */ 224 for ( i=0; i< x_len; i++) { 225 #pragma omp parallel for 226 for (int i=0; i< x_len; i++) { 225 227 double x_value = *(double *)(x->data + i*x->strides[0]); 226 228 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CLamellarFFHGModel.cpp
r2605da22 r0b082f3 170 170 return NULL; 171 171 } 172 #pragma omp parallel for 172 173 for (int i = 0; i < q->dimensions[0]; i++){ 173 174 double q_value = *(double *)(q->data + i*q->strides[0]); … … 187 188 { 188 189 PyArrayObject *result; 189 int i,x_len, y_len, dims[1];190 int x_len, y_len, dims[1]; 190 191 //check validity of input vectors 191 192 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 211 212 212 213 /* Do the calculation. */ 213 for ( i=0; i< x_len; i++) { 214 #pragma omp parallel for 215 for (int i=0; i< x_len; i++) { 214 216 double x_value = *(double *)(x->data + i*x->strides[0]); 215 217 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CLamellarModel.cpp
r2605da22 r0b082f3 165 165 return NULL; 166 166 } 167 #pragma omp parallel for 167 168 for (int i = 0; i < q->dimensions[0]; i++){ 168 169 double q_value = *(double *)(q->data + i*q->strides[0]); … … 182 183 { 183 184 PyArrayObject *result; 184 int i,x_len, y_len, dims[1];185 int x_len, y_len, dims[1]; 185 186 //check validity of input vectors 186 187 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 206 207 207 208 /* Do the calculation. */ 208 for ( i=0; i< x_len; i++) { 209 #pragma omp parallel for 210 for (int i=0; i< x_len; i++) { 209 211 double x_value = *(double *)(x->data + i*x->strides[0]); 210 212 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CLamellarPCrystalModel.cpp
r2605da22 r0b082f3 168 168 return NULL; 169 169 } 170 #pragma omp parallel for 170 171 for (int i = 0; i < q->dimensions[0]; i++){ 171 172 double q_value = *(double *)(q->data + i*q->strides[0]); … … 185 186 { 186 187 PyArrayObject *result; 187 int i,x_len, y_len, dims[1];188 int x_len, y_len, dims[1]; 188 189 //check validity of input vectors 189 190 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 209 210 210 211 /* Do the calculation. */ 211 for ( i=0; i< x_len; i++) { 212 #pragma omp parallel for 213 for (int i=0; i< x_len; i++) { 212 214 double x_value = *(double *)(x->data + i*x->strides[0]); 213 215 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CLamellarPSHGModel.cpp
r2605da22 r0b082f3 176 176 return NULL; 177 177 } 178 #pragma omp parallel for 178 179 for (int i = 0; i < q->dimensions[0]; i++){ 179 180 double q_value = *(double *)(q->data + i*q->strides[0]); … … 193 194 { 194 195 PyArrayObject *result; 195 int i,x_len, y_len, dims[1];196 int x_len, y_len, dims[1]; 196 197 //check validity of input vectors 197 198 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 217 218 218 219 /* Do the calculation. */ 219 for ( i=0; i< x_len; i++) { 220 #pragma omp parallel for 221 for (int i=0; i< x_len; i++) { 220 222 double x_value = *(double *)(x->data + i*x->strides[0]); 221 223 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CLamellarPSModel.cpp
r2605da22 r0b082f3 171 171 return NULL; 172 172 } 173 #pragma omp parallel for 173 174 for (int i = 0; i < q->dimensions[0]; i++){ 174 175 double q_value = *(double *)(q->data + i*q->strides[0]); … … 188 189 { 189 190 PyArrayObject *result; 190 int i,x_len, y_len, dims[1];191 int x_len, y_len, dims[1]; 191 192 //check validity of input vectors 192 193 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 212 213 213 214 /* Do the calculation. */ 214 for ( i=0; i< x_len; i++) { 215 #pragma omp parallel for 216 for (int i=0; i< x_len; i++) { 215 217 double x_value = *(double *)(x->data + i*x->strides[0]); 216 218 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CLogNormal.cpp
r2605da22 r0b082f3 160 160 return NULL; 161 161 } 162 #pragma omp parallel for 162 163 for (int i = 0; i < q->dimensions[0]; i++){ 163 164 double q_value = *(double *)(q->data + i*q->strides[0]); … … 177 178 { 178 179 PyArrayObject *result; 179 int i,x_len, y_len, dims[1];180 int x_len, y_len, dims[1]; 180 181 //check validity of input vectors 181 182 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 201 202 202 203 /* Do the calculation. */ 203 for ( i=0; i< x_len; i++) { 204 #pragma omp parallel for 205 for (int i=0; i< x_len; i++) { 204 206 double x_value = *(double *)(x->data + i*x->strides[0]); 205 207 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CLorentzian.cpp
r2605da22 r0b082f3 160 160 return NULL; 161 161 } 162 #pragma omp parallel for 162 163 for (int i = 0; i < q->dimensions[0]; i++){ 163 164 double q_value = *(double *)(q->data + i*q->strides[0]); … … 177 178 { 178 179 PyArrayObject *result; 179 int i,x_len, y_len, dims[1];180 int x_len, y_len, dims[1]; 180 181 //check validity of input vectors 181 182 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 201 202 202 203 /* Do the calculation. */ 203 for ( i=0; i< x_len; i++) { 204 #pragma omp parallel for 205 for (int i=0; i< x_len; i++) { 204 206 double x_value = *(double *)(x->data + i*x->strides[0]); 205 207 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CMultiShellModel.cpp
r2605da22 r0b082f3 174 174 return NULL; 175 175 } 176 #pragma omp parallel for 176 177 for (int i = 0; i < q->dimensions[0]; i++){ 177 178 double q_value = *(double *)(q->data + i*q->strides[0]); … … 191 192 { 192 193 PyArrayObject *result; 193 int i,x_len, y_len, dims[1];194 int x_len, y_len, dims[1]; 194 195 //check validity of input vectors 195 196 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 215 216 216 217 /* Do the calculation. */ 217 for ( i=0; i< x_len; i++) { 218 #pragma omp parallel for 219 for (int i=0; i< x_len; i++) { 218 220 double x_value = *(double *)(x->data + i*x->strides[0]); 219 221 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/COnionModel.cpp
r2605da22 r0b082f3 246 246 return NULL; 247 247 } 248 #pragma omp parallel for 248 249 for (int i = 0; i < q->dimensions[0]; i++){ 249 250 double q_value = *(double *)(q->data + i*q->strides[0]); … … 263 264 { 264 265 PyArrayObject *result; 265 int i,x_len, y_len, dims[1];266 int x_len, y_len, dims[1]; 266 267 //check validity of input vectors 267 268 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 287 288 288 289 /* Do the calculation. */ 289 for ( i=0; i< x_len; i++) { 290 #pragma omp parallel for 291 for (int i=0; i< x_len; i++) { 290 292 double x_value = *(double *)(x->data + i*x->strides[0]); 291 293 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CParallelepipedModel.cpp
r2605da22 r0b082f3 185 185 return NULL; 186 186 } 187 #pragma omp parallel for 187 188 for (int i = 0; i < q->dimensions[0]; i++){ 188 189 double q_value = *(double *)(q->data + i*q->strides[0]); … … 202 203 { 203 204 PyArrayObject *result; 204 int i,x_len, y_len, dims[1];205 int x_len, y_len, dims[1]; 205 206 //check validity of input vectors 206 207 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 226 227 227 228 /* Do the calculation. */ 228 for ( i=0; i< x_len; i++) { 229 #pragma omp parallel for 230 for (int i=0; i< x_len; i++) { 229 231 double x_value = *(double *)(x->data + i*x->strides[0]); 230 232 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CPearlNecklaceModel.cpp
r2605da22 r0b082f3 172 172 return NULL; 173 173 } 174 #pragma omp parallel for 174 175 for (int i = 0; i < q->dimensions[0]; i++){ 175 176 double q_value = *(double *)(q->data + i*q->strides[0]); … … 189 190 { 190 191 PyArrayObject *result; 191 int i,x_len, y_len, dims[1];192 int x_len, y_len, dims[1]; 192 193 //check validity of input vectors 193 194 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 213 214 214 215 /* Do the calculation. */ 215 for ( i=0; i< x_len; i++) { 216 #pragma omp parallel for 217 for (int i=0; i< x_len; i++) { 216 218 double x_value = *(double *)(x->data + i*x->strides[0]); 217 219 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CPoly_GaussCoil.cpp
r2605da22 r0b082f3 164 164 return NULL; 165 165 } 166 #pragma omp parallel for 166 167 for (int i = 0; i < q->dimensions[0]; i++){ 167 168 double q_value = *(double *)(q->data + i*q->strides[0]); … … 181 182 { 182 183 PyArrayObject *result; 183 int i,x_len, y_len, dims[1];184 int x_len, y_len, dims[1]; 184 185 //check validity of input vectors 185 186 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 205 206 206 207 /* Do the calculation. */ 207 for ( i=0; i< x_len; i++) { 208 #pragma omp parallel for 209 for (int i=0; i< x_len; i++) { 208 210 double x_value = *(double *)(x->data + i*x->strides[0]); 209 211 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CRPAModel.cpp
r2605da22 r0b082f3 189 189 return NULL; 190 190 } 191 #pragma omp parallel for 191 192 for (int i = 0; i < q->dimensions[0]; i++){ 192 193 double q_value = *(double *)(q->data + i*q->strides[0]); … … 206 207 { 207 208 PyArrayObject *result; 208 int i,x_len, y_len, dims[1];209 int x_len, y_len, dims[1]; 209 210 //check validity of input vectors 210 211 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 230 231 231 232 /* Do the calculation. */ 232 for ( i=0; i< x_len; i++) { 233 #pragma omp parallel for 234 for (int i=0; i< x_len; i++) { 233 235 double x_value = *(double *)(x->data + i*x->strides[0]); 234 236 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CReflAdvModel.cpp
rd760c0a r0b082f3 231 231 return NULL; 232 232 } 233 #pragma omp parallel for 233 234 for (int i = 0; i < q->dimensions[0]; i++){ 234 235 double q_value = *(double *)(q->data + i*q->strides[0]); … … 248 249 { 249 250 PyArrayObject *result; 250 int i,x_len, y_len, dims[1];251 int x_len, y_len, dims[1]; 251 252 //check validity of input vectors 252 253 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 272 273 273 274 /* Do the calculation. */ 274 for ( i=0; i< x_len; i++) { 275 #pragma omp parallel for 276 for (int i=0; i< x_len; i++) { 275 277 double x_value = *(double *)(x->data + i*x->strides[0]); 276 278 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CReflModel.cpp
rd760c0a r0b082f3 207 207 return NULL; 208 208 } 209 #pragma omp parallel for 209 210 for (int i = 0; i < q->dimensions[0]; i++){ 210 211 double q_value = *(double *)(q->data + i*q->strides[0]); … … 224 225 { 225 226 PyArrayObject *result; 226 int i,x_len, y_len, dims[1];227 int x_len, y_len, dims[1]; 227 228 //check validity of input vectors 228 229 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 248 249 249 250 /* Do the calculation. */ 250 for ( i=0; i< x_len; i++) { 251 #pragma omp parallel for 252 for (int i=0; i< x_len; i++) { 251 253 double x_value = *(double *)(x->data + i*x->strides[0]); 252 254 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CSCCrystalModel.cpp
r2605da22 r0b082f3 179 179 return NULL; 180 180 } 181 #pragma omp parallel for 181 182 for (int i = 0; i < q->dimensions[0]; i++){ 182 183 double q_value = *(double *)(q->data + i*q->strides[0]); … … 196 197 { 197 198 PyArrayObject *result; 198 int i,x_len, y_len, dims[1];199 int x_len, y_len, dims[1]; 199 200 //check validity of input vectors 200 201 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 220 221 221 222 /* Do the calculation. */ 222 for ( i=0; i< x_len; i++) { 223 #pragma omp parallel for 224 for (int i=0; i< x_len; i++) { 223 225 double x_value = *(double *)(x->data + i*x->strides[0]); 224 226 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CSLDCalFunc.cpp
r2605da22 r0b082f3 166 166 return NULL; 167 167 } 168 #pragma omp parallel for 168 169 for (int i = 0; i < q->dimensions[0]; i++){ 169 170 double q_value = *(double *)(q->data + i*q->strides[0]); … … 183 184 { 184 185 PyArrayObject *result; 185 int i,x_len, y_len, dims[1];186 int x_len, y_len, dims[1]; 186 187 //check validity of input vectors 187 188 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 207 208 208 209 /* Do the calculation. */ 209 for ( i=0; i< x_len; i++) { 210 #pragma omp parallel for 211 for (int i=0; i< x_len; i++) { 210 212 double x_value = *(double *)(x->data + i*x->strides[0]); 211 213 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CSchulz.cpp
r2605da22 r0b082f3 160 160 return NULL; 161 161 } 162 #pragma omp parallel for 162 163 for (int i = 0; i < q->dimensions[0]; i++){ 163 164 double q_value = *(double *)(q->data + i*q->strides[0]); … … 177 178 { 178 179 PyArrayObject *result; 179 int i,x_len, y_len, dims[1];180 int x_len, y_len, dims[1]; 180 181 //check validity of input vectors 181 182 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 201 202 202 203 /* Do the calculation. */ 203 for ( i=0; i< x_len; i++) { 204 #pragma omp parallel for 205 for (int i=0; i< x_len; i++) { 204 206 double x_value = *(double *)(x->data + i*x->strides[0]); 205 207 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CSphereModel.cpp
r2605da22 r0b082f3 165 165 return NULL; 166 166 } 167 #pragma omp parallel for 167 168 for (int i = 0; i < q->dimensions[0]; i++){ 168 169 double q_value = *(double *)(q->data + i*q->strides[0]); … … 182 183 { 183 184 PyArrayObject *result; 184 int i,x_len, y_len, dims[1];185 int x_len, y_len, dims[1]; 185 186 //check validity of input vectors 186 187 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 206 207 207 208 /* Do the calculation. */ 208 for ( i=0; i< x_len; i++) { 209 #pragma omp parallel for 210 for (int i=0; i< x_len; i++) { 209 211 double x_value = *(double *)(x->data + i*x->strides[0]); 210 212 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CSphereSLDModel.cpp
r2605da22 r0b082f3 223 223 return NULL; 224 224 } 225 #pragma omp parallel for 225 226 for (int i = 0; i < q->dimensions[0]; i++){ 226 227 double q_value = *(double *)(q->data + i*q->strides[0]); … … 240 241 { 241 242 PyArrayObject *result; 242 int i,x_len, y_len, dims[1];243 int x_len, y_len, dims[1]; 243 244 //check validity of input vectors 244 245 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 264 265 265 266 /* Do the calculation. */ 266 for ( i=0; i< x_len; i++) { 267 #pragma omp parallel for 268 for (int i=0; i< x_len; i++) { 267 269 double x_value = *(double *)(x->data + i*x->strides[0]); 268 270 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CSquareWellStructure.cpp
r2605da22 r0b082f3 164 164 return NULL; 165 165 } 166 #pragma omp parallel for 166 167 for (int i = 0; i < q->dimensions[0]; i++){ 167 168 double q_value = *(double *)(q->data + i*q->strides[0]); … … 181 182 { 182 183 PyArrayObject *result; 183 int i,x_len, y_len, dims[1];184 int x_len, y_len, dims[1]; 184 185 //check validity of input vectors 185 186 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 205 206 206 207 /* Do the calculation. */ 207 for ( i=0; i< x_len; i++) { 208 #pragma omp parallel for 209 for (int i=0; i< x_len; i++) { 208 210 double x_value = *(double *)(x->data + i*x->strides[0]); 209 211 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CStackedDisksModel.cpp
r2605da22 r0b082f3 184 184 return NULL; 185 185 } 186 #pragma omp parallel for 186 187 for (int i = 0; i < q->dimensions[0]; i++){ 187 188 double q_value = *(double *)(q->data + i*q->strides[0]); … … 201 202 { 202 203 PyArrayObject *result; 203 int i,x_len, y_len, dims[1];204 int x_len, y_len, dims[1]; 204 205 //check validity of input vectors 205 206 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 225 226 226 227 /* Do the calculation. */ 227 for ( i=0; i< x_len; i++) { 228 #pragma omp parallel for 229 for (int i=0; i< x_len; i++) { 228 230 double x_value = *(double *)(x->data + i*x->strides[0]); 229 231 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CStickyHSStructure.cpp
r2605da22 r0b082f3 164 164 return NULL; 165 165 } 166 #pragma omp parallel for 166 167 for (int i = 0; i < q->dimensions[0]; i++){ 167 168 double q_value = *(double *)(q->data + i*q->strides[0]); … … 181 182 { 182 183 PyArrayObject *result; 183 int i,x_len, y_len, dims[1];184 int x_len, y_len, dims[1]; 184 185 //check validity of input vectors 185 186 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 205 206 206 207 /* Do the calculation. */ 207 for ( i=0; i< x_len; i++) { 208 #pragma omp parallel for 209 for (int i=0; i< x_len; i++) { 208 210 double x_value = *(double *)(x->data + i*x->strides[0]); 209 211 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CTriaxialEllipsoidModel.cpp
r2605da22 r0b082f3 185 185 return NULL; 186 186 } 187 #pragma omp parallel for 187 188 for (int i = 0; i < q->dimensions[0]; i++){ 188 189 double q_value = *(double *)(q->data + i*q->strides[0]); … … 202 203 { 203 204 PyArrayObject *result; 204 int i,x_len, y_len, dims[1];205 int x_len, y_len, dims[1]; 205 206 //check validity of input vectors 206 207 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 226 227 227 228 /* Do the calculation. */ 228 for ( i=0; i< x_len; i++) { 229 #pragma omp parallel for 230 for (int i=0; i< x_len; i++) { 229 231 double x_value = *(double *)(x->data + i*x->strides[0]); 230 232 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/CVesicleModel.cpp
r2605da22 r0b082f3 169 169 return NULL; 170 170 } 171 #pragma omp parallel for 171 172 for (int i = 0; i < q->dimensions[0]; i++){ 172 173 double q_value = *(double *)(q->data + i*q->strides[0]); … … 186 187 { 187 188 PyArrayObject *result; 188 int i,x_len, y_len, dims[1];189 int x_len, y_len, dims[1]; 189 190 //check validity of input vectors 190 191 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 210 211 211 212 /* Do the calculation. */ 212 for ( i=0; i< x_len; i++) { 213 #pragma omp parallel for 214 for (int i=0; i< x_len; i++) { 213 215 double x_value = *(double *)(x->data + i*x->strides[0]); 214 216 double y_value = *(double *)(y->data + i*y->strides[0]); -
sansmodels/src/sans/models/c_models/classTemplate.txt
r2605da22 r0b082f3 150 150 return NULL; 151 151 } 152 #pragma omp parallel for 152 153 for (int i = 0; i < q->dimensions[0]; i++){ 153 154 double q_value = *(double *)(q->data + i*q->strides[0]); … … 167 168 { 168 169 PyArrayObject *result; 169 int i,x_len, y_len, dims[1];170 int x_len, y_len, dims[1]; 170 171 //check validity of input vectors 171 172 if (x->nd != 1 || x->descr->type_num != PyArray_DOUBLE … … 191 192 192 193 /* Do the calculation. */ 193 for ( i=0; i< x_len; i++) { 194 #pragma omp parallel for 195 for (int i=0; i< x_len; i++) { 194 196 double x_value = *(double *)(x->data + i*x->strides[0]); 195 197 double y_value = *(double *)(y->data + i*y->strides[0]);
Note: See TracChangeset
for help on using the changeset viewer.