File: mass_surface_fractal.c

package info (click to toggle)
mccode 3.5.19%2Bds5-2
  • links: PTS, VCS
  • area: main
  • in suites: forky, 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 (17 lines) | stat: -rw-r--r-- 513 bytes parent folder | download | duplicates (6) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
 
static double
Iq(double q,
          double fractal_dim_mass,
          double fractal_dim_surf,
          double rg_cluster,
          double rg_primary)
{
     //computation
    const double Dm = 0.5*fractal_dim_mass;
    const double Dt = 0.5*(6.0 - (fractal_dim_mass + fractal_dim_surf));

    const double t1 = Dm==0. ? 1.0 : pow(1.0 + square(q*rg_cluster)/(3.0*Dm), -Dm);
    const double t2 = Dt==0. ? 1.0 : pow(1.0 + square(q*rg_primary)/(3.0*Dt), -Dt);
    const double form = t1*t2;

    return form;
}