File: sphere.c

package info (click to toggle)
mccode 3.5.19%2Bds5-2
  • links: PTS, VCS
  • area: main
  • in suites: sid, trixie
  • size: 1,113,256 kB
  • sloc: ansic: 40,697; python: 25,137; yacc: 8,438; sh: 5,405; javascript: 4,596; lex: 1,632; cpp: 742; perl: 296; lisp: 273; makefile: 226; fortran: 132
file content (20 lines) | stat: -rw-r--r-- 488 bytes parent folder | download | duplicates (6) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
 
static double form_volume(double radius)
{
    return M_4PI_3*cube(radius);
}

static double
radius_effective(int mode, double radius)
{
    // case 1: radius
    return radius;
}

static void Fq(double q, double *f1, double *f2, double sld, double sld_solvent, double radius)
{
    const double bes = sas_3j1x_x(q*radius);
    const double contrast = (sld - sld_solvent);
    const double form = contrast * form_volume(radius) * bes;
    *f1 = 1.0e-2*form;
    *f2 = 1.0e-4*form*form;
}