pdf - Madagascar

6
seabeam/interpolate.c
1 /∗ Multi−dimensional missing data i n t e r p o l a t i o n . ∗/
2
3 #include
4
5 int main ( int argc , char∗ argv [ ] )
6 {
7 int na , ia , n i t e r , n , i ;
8 float a0 , ∗mm, ∗ zero ;
9 bool prec , ∗known ;
10 char ∗ l a g f i l e ;
11 s f f i l t e r aa ;
12 s f f i l e in , out , f i l t , lag , mask ;
13
14 s f i n i t ( argc , argv ) ;
15 in = s f i n p u t ( ” in ” ) ;
16 f i l t = s f i n p u t ( ” f i l t ” ) ;
17 /∗ f i l t e r f o r i n v e r s e model covariance ∗/
18 out = s f o u t p u t ( ” out ” ) ;
19
20 n = s f f i l e s i z e ( in ) ;
21
22 i f ( ! s f g e t b o o l ( ” prec ”,& prec ) ) prec=true ;
23 /∗ I f y , use p r e c o n d i t i o n i n g ∗/
24 i f ( ! s f g e t i n t ( ” n i t e r ”,& n i t e r ) ) n i t e r =100;
25 /∗ Number o f i t e r a t i o n s ∗/
26
27 i f ( ! s f h i s t i n t ( f i l t , ”n1”,&na ) ) s f e r r o r ( ”No n1=” ) ;
28 aa = s f a l l o c a t e h e l i x ( na ) ;
29
30 i f ( ! s f h i s t f l o a t ( f i l t , ”a0”,&a0 ) ) a0 =1.;
31
32 /∗ Get f i l t e r l a g s ∗/
33 i f (NULL == ( l a g f i l e = s f h i s t s t r i n g ( f i l t , ” l a g ” ) ) ) {
34 for ( i a =0; i a < na ; i a++) {
35 aa−>l a g [ i a ]= i a +1;
36 }
37 } else {
38 l a g = s f i n p u t ( l a g f i l e ) ;
39 s f i n t r e a d ( aa−>lag , na , l a g ) ;
40 s f f i l e c l o s e ( l a g ) ;
41 }
42
43 /∗ Get f i l t e r v a l u e s ∗/