9ef2dff3b26ef4651538a0bc9d9484fdeb2357be
[unres.git] / source / wham / src / proc_proc.c
1 #include <stdlib.h>
2 #include <math.h>
3
4 #ifdef LINUX
5 #ifdef PGI
6 void proc_proc_(long int *f, int *i)
7 #else
8 void proc_proc__(long int *f, int *i)
9 #endif
10 #endif
11 #ifdef SGI
12 void proc_proc_(long int *f, int *i)
13 #endif
14 #ifdef WIN
15 void _stdcall PROC_PROC(long int *f, int *i)
16 #endif
17 #if defined(AIX) || defined(WINPGI) 
18 void proc_proc(long int *f, int *i)
19 #endif
20   
21 {
22 static long int NaNQ;
23 static long int NaNQm;
24
25 if(*i==-1)
26  {
27  NaNQ=*f;
28  NaNQm=0xffffffff;
29  return;
30  }
31 *i=0;
32 if(*f==NaNQ)
33  *i=1;
34 if(*f==NaNQm)
35  *i=1;
36 }
37
38
39 #ifdef LINUX
40 void proc_conv__(char *buf, int *i, int n)
41 #endif
42 #ifdef SGI
43 void proc_conv_(char *buf, int *i, int n)
44 #endif
45 #if defined(AIX) || defined(WINPGI)
46 void proc_conv(char *buf, int *i, int n)
47 #endif
48 #ifdef WIN
49 void _stdcall PROC_CONV(char *buf, int *i, int n)
50 #endif
51 {
52 int j;
53
54 sscanf(buf,"%d",&j);
55 *i=j;
56 return;
57 }
58
59 #ifdef LINUX
60 void proc_conv_r__(char *buf, int *i, int n)
61 #endif
62 #ifdef SGI
63 void proc_conv_r_(char *buf, int *i, int n)
64 #endif
65 #if defined(AIX) || defined(WINPGI)
66 void proc_conv_r(char *buf, int *i, int n)
67 #endif
68 #ifdef WIN
69 void _stdcall PROC_CONV_R(char *buf, int *i, int n)
70 #endif
71
72 {
73
74 /* sprintf(buf,"%d",*i); */
75
76 return;
77 }
78
79 #ifndef IMSL
80 #ifdef LINUX
81 void dsvrgp__(int *n, double *tab1, double *tab2, int *itab)
82 #endif
83 #ifdef SGI
84 void dsvrgp_(int *n, double *tab1, double *tab2, int *itab)
85 #endif
86 #if defined(AIX) || defined(WINPGI)
87 void dsvrgp(int *n, double *tab1, double *tab2, int *itab)
88 #endif
89 #ifdef WIN
90 void _stdcall DSVRGP(int *n, double *tab1, double *tab2, int *itab)
91 #endif
92 {
93 double t;
94 int i,j,k;
95
96 if(tab1 != tab2)
97  {
98  for(i=0; i<*n; i++)
99   tab2[i]=tab1[i];
100  }
101 k=0;
102 while(k<*n-1)
103  {
104  j=k;
105  t=tab2[k];
106  for(i=k+1; i<*n; i++)
107   if(t>tab2[i])
108    {
109    j=i;
110    t=tab2[i];
111    }
112  if(j!=k)
113   {
114   tab2[j]=tab2[k];
115   tab2[k]=t;
116   i=itab[j];
117   itab[j]=itab[k];
118   itab[k]=i;
119   }
120  k++;
121  }
122 }
123 #endif