File: eulerc.c

package info (click to toggle)
scilab 2.2-4
  • links: PTS
  • area: non-free
  • in suites: hamm
  • size: 31,472 kB
  • ctags: 21,963
  • sloc: fortran: 110,983; ansic: 89,717; makefile: 3,016; sh: 1,892; csh: 150; cpp: 101
file content (41 lines) | stat: -rw-r--r-- 875 bytes parent folder | download
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
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
#include <string.h>
#include <malloc.h>

#include "../machine.h"

void C2F(eulerc)(alpho,la2,lp2,ls2,ma,mm,n,phiw,sigma,sigmadim)
int *alpho,*la2,*lp2,*ls2,*ma,*mm,*n,*phiw,**sigma,*sigmadim;
{
  int *alphi,*beta,*deg,*hp1,*hp2,*ihp1,*ihp2,*imin;
  int *sigmaw;
  int edg,eul,i,ma1;
  int isize = sizeof(int);
  
  ma1 = *ma + 1;
  
  if ((sigmaw = (int *)malloc(ma1 * isize)) == NULL) {
    cerro("Running out of memory");
    return;
  }

  F2C(euler)(&eul,alpho,la2,lp2,ls2,ma,&ma1,mm,n,phiw,sigmaw);

  if (eul == 0)
    *sigmadim = 0;
  else {
    if ((*sigma = (int *)malloc(*ma * isize)) == NULL) {
      cerro("Running out of memory");
      return;
    }
    *sigmadim = *ma;
    edg = sigmaw[*ma];
    (*sigma)[0] = edg;
    for (i = 2; i <= *ma; i++) {
      if (edg == -1) break;
      edg = sigmaw[edg-1];
      (*sigma)[i-1] = edg;
    }
  }

  free(sigmaw);
}