1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
|
#include <math.h>
double Iq(double q, double porod_scale, double porod_exp, double peak_scale,
double correlation_length, double peak_pos, double width_exp,
double shape_exp);
#pragma acc routine seq
double Iq(double q, double porod_scale, double porod_exp, double peak_scale,
double correlation_length, double peak_pos, double width_exp,
double shape_exp) {
double z = fabs(q - peak_pos) * correlation_length;
double inten = pow((porod_scale / pow(q, porod_exp)) + (peak_scale / (1.0 + pow(z, width_exp))), shape_exp);
return inten;
}
|
