File: bspline_util.c

package info (click to toggle)
python-scipy 0.5.2-0.1
  • links: PTS
  • area: main
  • in suites: etch, etch-m68k
  • size: 33,888 kB
  • ctags: 44,231
  • sloc: ansic: 156,256; cpp: 90,347; python: 89,604; fortran: 73,083; sh: 1,318; objc: 424; makefile: 342
file content (26 lines) | stat: -rw-r--r-- 549 bytes parent folder | download | duplicates (10) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
 
#include <math.h>
#include <stddef.h>
#include <stdlib.h>
#include <string.h>
#include <stdio.h>

void compute_root_from_lambda(double, double *, double *);


void
compute_root_from_lambda(lambda, r, omega)
     double lambda;
     double *r;
     double *omega;
{
    double xi;
    double tmp, tmp2;

    tmp = sqrt(3 + 144*lambda);
    xi = 1 - 96*lambda + 24*lambda * tmp;
    *omega = atan(sqrt((144*lambda - 1.0)/xi));
    tmp2 = sqrt(xi);
    *r = (24*lambda - 1 - tmp2)/(24*lambda) \
	* sqrt((48*lambda + 24*lambda*tmp))/tmp2;
    return;
}