home *** CD-ROM | disk | FTP | other *** search
/ Source Code 1994 March / Source_Code_CD-ROM_Walnut_Creek_March_1994.iso / compsrcs / misc / volume38 / lic / part06 < prev    next >
Encoding:
Text File  |  1993-08-11  |  64.9 KB  |  1,581 lines
  1. Newsgroups: comp.sources.misc
  2. From: casey@gauss.llnl.gov (Casey Leedom)
  3. Subject: v38i109:  lic - LLNL Line Integral Convolution, v1.2, Part06/10
  4. Message-ID: <1993Aug12.013908.14233@sparky.sterling.com>
  5. X-Md4-Signature: 64cf321283b575bf139fd82cad4b7ae6
  6. Sender: kent@sparky.sterling.com (Kent Landfield)
  7. Organization: Sterling Software
  8. Date: Thu, 12 Aug 1993 01:39:08 GMT
  9. Approved: kent@sparky.sterling.com
  10.  
  11. Submitted-by: casey@gauss.llnl.gov (Casey Leedom)
  12. Posting-number: Volume 38, Issue 109
  13. Archive-name: lic/part06
  14. Environment: UNIX
  15.  
  16. #! /bin/sh
  17. # This is a shell archive.  Remove anything before this line, then feed it
  18. # into a shell via "sh file" or similar.  To overwrite existing files,
  19. # type "sh file -c".
  20. # Contents:  lic.1.2/doc/siggraph93/p.frame.UU.C
  21. #   lic.1.2/liblic/Filters.c lic.1.2/lic/lic.c
  22. # Wrapped by kent@sparky on Wed Aug 11 19:38:05 1993
  23. PATH=/bin:/usr/bin:/usr/ucb:/usr/local/bin:/usr/lbin ; export PATH
  24. echo If this archive is complete, you will see the following message:
  25. echo '          "shar: End of archive 6 (of 10)."'
  26. if test -f 'lic.1.2/doc/siggraph93/p.frame.UU.C' -a "${1}" != "-c" ; then 
  27.   echo shar: Will not clobber existing file \"'lic.1.2/doc/siggraph93/p.frame.UU.C'\"
  28. else
  29.   echo shar: Extracting \"'lic.1.2/doc/siggraph93/p.frame.UU.C'\" \(32376 characters\)
  30.   sed "s/^X//" >'lic.1.2/doc/siggraph93/p.frame.UU.C' <<'END_OF_FILE'
  31. XM   $+', !554  555  "JJH * ')#Z@    W' CZ'%1H92!U<'!E<B!L969T
  32. XM(&AA;F0@<75A<G1E<B!O9B!T:&4@8VER8W5L87(@=F5C=&]R(!P+<@ .55( 
  33. XM!554  *JJ@   <D/J [9 $!F:65L9"!I<R!C;VYV;VQV960@=7-I;F<@3$E#
  34. XM(&]V97(@0G5T=&5R=V]R=&@@;&]W+7!A<W,@9FEL=&5R960@    %U50  55
  35. XM5  "JJH   ')#ZA=+ !!=VAI=&4@;F]I<V4@=VET:"!C=71O9F8@9G)E<75E
  36. XM;F-I97,@;V8@,3(X+" X-BP@-C0L(&%N9" S,B H;&5F="!H "!53@ %550 
  37. XM JJJ $ !R0^H#N( &G1O(')I9VAT+"!T;W @=&\@8F]T=&]M*2X+ #8    V
  38. XM    [_G; (N)-P[$#ML.F0[+#I<        . '!O           .U0[R    
  39. XM         '); !YO4 !EFCP 9_.N  K+RP[%#N(.P [1#K8        0 ');
  40. XM !YO4 !KO_@       %A<@"597%U86Q;8VAA<EML72QL97%;:6YD97AE<ULP
  41. XM+#$L=&EM97-;8VAA<EMI72QS=')I;F=;(A%S=6-H$71H8701(ETL8VAA<EMS
  42. XM75TL8VAA<EMI75TL;&5S<W1H86Y;8VAA<EM,72QI;F1E>&5S6S L,2QC:&%R
  43. XM6W-=+'!L=7-;8VAA<EMI72QN=6U;,2PB,2)=75U=75TL #8    V    [_G;
  44. XM 0_S@0[&#N$.E@  #I0        . &%R           .U@[N            
  45. XM &UA #8#$@%9"H( [_G; #&M3@['#MH   [/#IH        . &%R     /_^
  46. XM'_\.R0[3   .FPBS U0  &%R !W5P !-0XL J?'Q !4X' [(#MD.T@  #K8 
  47. XM       / "QC#K\.P  ? %8 !+"Z )YXS@ 5.$P.R0[B   .V@['        
  48. XM$ !N*@!N/+T $<67      (!9ET PF5Q=6%L6W1I;65S6V-H87);:UTL:61;
  49. XM8VAA<EMW75U=+&-R;W-S6V]V97);<&QU<UMN=6U;,2PB,2)=+&-O<UMI9%MT
  50. XM:6UE<ULH*FXJ*6-H87);8UTL8VAA<EMW75U=75TL;G5M6S(L(C(B75TL;W9E
  51. XM<EMP;'5S6VYU;5LQ+"(Q(ETL8V]S6VED6W!L=7-;=&EM97-;8VAA<EMD72QC
  52. XM:&%R6W==72QC:&%R6V)E=&%=75U=72QN=6U;,BPB,B)=75U= !Z",  E2%  
  53. XMQI/" &16P@[*#N(.S0[4#L\        0     (',$0!99*D      @%55 ,]
  54. XM97%U86Q;<')O;7!T6UTL=&EM97-;;W9E<EMN=6U;,2PB,2)=+&YU;5LT+"(T
  55. XM(EU=+&ED6V%T;W!;<&QU<UMC:&%R6V)=+&UI;G5S6W!L=7-;*"IN*BEC:&%R
  56. XM6V%=+&]V97);<VEN6W!L=7-;*"IN*BEI9%MT:6UE<UMC:&%R6V)=+&-H87);
  57. XM8UU=72QM:6YU<UMS:6Y;:61;=&EM97-;8VAA<EMA72QC:&%R6V-=75U=75U=
  58. XM+&-H87);8UU=75U=+'!L=7-;<W1R:6YG6R(1(ETL;W9E<EMS:6Y;<&QU<ULH
  59. XM*FXJ*6ED6W!L=7-;=&EM97-;8VAA<EMB72QC:&%R6V1=72QC:&%R6V)E=&%=
  60. XM75TL;6EN=7-;<VEN6VED6W!L=7-;=&EM97-;8VAA<EMA72QC:&%R6V1=72QC
  61. XM:&%R6V)E=&%=75U=75U=+&-H87);9%U=72QP;'5S6W-T<FEN9ULB$2)=+&]V
  62. XM97);<VEN6W!L=7-;*"IN*BEI9%MP;'5S6W1I;65S6V-H87);8ETL:61;<&QU
  63. XM<ULH*FXJ*6-H87);8UTL;6EN=7-;8VAA<EMD75U=75TL;6EN=7-;8VAA<EMB
  64. XM971A75U=72QM:6YU<UMS:6Y;:61;<&QU<UMT:6UE<UMC:&%R6V%=+&ED6W!L
  65. XM=7-;*"IN*BEC:&%R6V-=+&UI;G5S6V-H87);9%U=75U=+&UI;G5S6V-H87);
  66. XM8F5T85U=75U=75U=+'1I;65S6VYU;5LR+"(R(ETL:61;<&QU<UMC:&%R6V-=
  67. XM+&UI;G5S6V-H87);9%U=75U=75TL<&QU<UMS=')I;F=;(A$B72QO=F5R6W-I
  68. XM;EMP;'5S6R@J;BHI:61;<&QU<UMT:6UE<UMC:&%R6V)=+&ED6W!L=7-;8VAA
  69. XM<EMC72QC:&%R6V1=75U=+&-H87);8F5T85U=72QM:6YU<UMS:6Y;:61;<&QU
  70. XM<UMT:6UE<UMC:&%R6V%=+&ED6W!L=7-;8VAA<EMC72QC:&%R6V1=75U=+&-H
  71. XM87);8F5T85U=75U=75TL=&EM97-;;G5M6S(L(C(B72QI9%MP;'5S6V-H87);
  72. XM8UTL8VAA<EMD75U=75U=75U=73$ -@   ;=\[0#O^=L!!H,2#LL.VP[$#K4.
  73. XMEP        X 87(           [9#O$             97(  P=D &5G+P!M
  74. XM!W, 1019#LP.X@[!#M<.J@  #M<.ZQ  6V, .8L= (G:4P     " 5TL B]E
  75. XM<75A;%MI;F1E>&5S6S L,2QC:&%R6W-=+&-H87);95U=+&QP87)E;ELH*FDR
  76. XM:2HI871O<%MT:6UE<UMC:&%R6VEN9G1Y72QS=')I;F=;(A%I9A$B72QC:&%R
  77. XM6U9=+'-T<FEN9ULB$7P2?!$B72QC:&%R6R@J;BHI95TL<W1R:6YG6R(1$1$1
  78. XM$1$1$1$1$1$1$1$1$2)=72QT:6UE<UMN=6U;,"PB,")=+'-T<FEN9ULB$6EF
  79. XM$2)=+&QE<W-T:&%N6R@J;BHI;W9E<ELH*FXJ*7!L=7-;*"IN*BEF;&]O<ELH
  80. XM*FXJ*6EN9&5X97-;*"IN*BDP+#$L8VAA<ELH*FXJ*5!=+&-H87);*"IN*BEC
  81. XM75U=+&UI;G5S6R@J;BHI:6YD97AE<ULH*FXJ*3 L,2QC:&%R6R@J;BHI4%TL
  82. XM8VAA<ELH*FXJ*6-=75U=+&EN9&5X97-;*"IN*BDP+#$L8VAA<ELH*FXJ*59=
  83. XM+&-H87);*"IN*BEC75U=+'1I;65S6VYU;5LH*FXJ*3 L(C B72QS=')I;F=;
  84. XM(A$1$2)=75U=+'1I;65S6V]V97);<&QU<UMF;&]O<EMI;F1E>&5S6S L,2QC
  85. XM:&%R6U!=+&-H87);8UU=72QM:6YU<UMI;F1E>&5S6S L,2QC:&%R6U!=+&-H
  86. XM87);8UU=75TL:6YD97AE<ULP+#$L8VAA<EM672QC:&%R6V-=75TL<W1R:6YG
  87. XM6R(1;W1H97)W:7-E(EU=75U=97-; !Z?Z0 %EJP )2JS "! \P[-#N(   [*
  88. XM#L\        0 %MC #$U0@ 7:!T      @%U<P _:6YT6W1I;65S6V-H87);
  89. XM:UTL:61;8VAA<EMW75TL8VAA<EMD72QC:&%R6W==72QC:&%R6V%=+&-H87);
  90. XM8EU=75U= #8   (@=L@ [_G; *^)-P[.#ML.W@[;#IT        . %LH    
  91. XM       .WP[B             %MB #0CV &S"H( \[A1 (_K!P[/#MH.QP  
  92. XM#IH        . &)E     /_LA L.S0[4   .FPBT U0  '-; -M@#  0$;8 
  93. XM%B28  JVW@[0#ML.J0  #JL        , %MD -M@#  0$;8 %B28  JVW@  
  94. XM                   !_P          <&P       %R6P &    !@    , 
  95. XM  !H :0/HETQ  4<",H<"^T VVGD !\MH0 6))@ "'WZ#M$.VP[%#M(.M@  
  96. XM      P     VVGD !\MH0 6))@ "'WZ                      (!    
  97. XM                 5MC  8    &     P   &@!IP^B97$ !1P(UQP+90#;
  98. XM48$ 9K_Y !8DF  '[\ .T@[;#M$.R ZV        #     #;48$ 9K_Y !8D
  99. XMF  '[\                       @(          &AA       !='( !@  
  100. XM  8    #    : &[#Z(H*@ %' C8' MT -N .  )6VP %B28  JVW@[3#ML.
  101. XMV@  #L<        , '5M -N .  )6VP %B28  JVW@                  
  102. XM   " P          *"H       %S6P &    !@    ,   !H :L/HF1E  4<
  103. XM"+X<"RH W  X !!9K0 6))@ "K;>#M0.VP[*   .SP        P :6X W  X
  104. XM !!9K0 6))@ "K;>                      ($          !;*       
  105. XM 5U=  8    &     P   &@!K@^B+&, !1P(W1P+*@      >8DW ._YVP 2
  106. XM   .U0[;   .\@[$        #  P+       >8DW ._YVP 2            
  107. XM             @@          ');       "<UL !554  555  "JJH * &W
  108. XM#ZA=70 X' B\'$-I<F-U;&%R(&%N9"!T=7)B=6QE;G0@9FQU:60@9'EN86UI
  109. XM8W,@=F5C=&]R(&9I96QD<R L8P .55( !554  *JJ@!  ;</J&]T "-I;6%G
  110. XM960@=7-I;F<@3$E#(&]V97(@=VAI=&4@;F]I<V4N"P .R@      ](,2 ._"
  111. XMCP ;   .U@[;   .[@[&        # !T:0      ](,2 ._"CP ;        
  112. XM                 @D          &AA       #:&$ !554  555  "JJH 
  113. XM* &Y#ZCYVP U' C '$$@='=O+61I;65N<VEO;F%L('9E8W1O<B!F:65L9"!S
  114. XM:&]W:6YG('1H92!L;V-A;"""  Y54@ %550  JJJ   !N0^H    4W-T<F5A
  115. XM;2!L:6YE('-T87)T:6YG(&EN(&-E;&P@*!L/JAMX&P^H&RP@&P^J&WD;#Z@;
  116. XM*2X@5&AE('9E8W1O<B!F:65L9"!I<R!T:&4@=7!P97(@     !=54  %550 
  117. XM JJJ $@!N0^H    3&QE9G0@8V]R;F5R(&]F('1H92!F;'5I9"!D>6YA;6EC
  118. XM<R!F:65L9"!I;B!F:6=U<F5S(!P)$APR' D6'"!A;F0@' DO'#0<"3 <+@MI
  119. XMY !M,RT ;[3U &CSOP R#Y0.UP[B#LP.Z ZJ#LP   [K$     "AK0P BJVW
  120. XM      (!   !C&EN6R@J;BHI=&EM97-;<W1R:6YG6R)F;W(B72QI9%MC;VUM
  121. XM85LH*FXJ*6-H87);*"IN*BEE72QC:&%R6R@J;BHI8UU=75TL:61;*"II,FDJ
  122. XM*6%T;W!;:61;=&EM97-;8VAA<EMT72QC:&%R6V]=+&-O;6UA6R@J;BHI8VAA
  123. XM<ELH*FXJ*7!=+&-H87);*"IN*BEY75U=72QI9%MT:6UE<UMC:&%R6V)=+&-H
  124. XM87);;UTL8VAA<EMT72QC:&%R6W1=+&-H87);;UTL8V]M;6%;*"IN*BEC:&%R
  125. XM6R@J;BHI;5TL8VAA<ELH*FXJ*7E=75U=+&ED6W1I;65S6V-H87);;%TL8VAA
  126. XM<EME72QC:&%R6V9=+&-O;6UA6R@J;BHI8VAA<ELH*FXJ*71=+&-H87);*"IN
  127. XM*BEX75U=72QI9%MT:6UE<UMC:&%R6W)=+&-H87);:5TL8VAA<EMG72QC:&%R
  128. XM6VA=+&-O;6UA6R@J;BHI8VAA<ELH*FXJ*71=+&-H87);*"IN*BEX75U=75U=
  129. XM70    "X40  N%$ [HDW .Z)-P[8#N<.PP  #K4        &            
  130. XM    ](,2 ._YVP 2   .V0[;   .\0[+        #  !MP      ](,2 ._Y
  131. XMVP 2                         @H          &%M       "8W0 !554
  132. XM  555  "JJH * '!#ZA55  ]' CS'%!H;W1O9W)A<&@@;V8@9FQO=V5R<R!P
  133. XM<F]C97-S960@=7-I;F<@3$E#('=I=&@@&P^J&TP;#Z@;((\ #E52  555  "
  134. XMJJH 0 '!#Z@    Y97%U86P@=&\@,"P@-2P@,3 @86YD(#(P("AL969T('1O
  135. XM(')I9VAT+"!T;W @=&\@8F]T=&]M*2X+5  +51X '#JB .'4@P 4*,X.V@[B
  136. XM#LD.TP['        $ !I;P!\/U\ *$ !      (!;&0!671I;65S6W-T<FEN
  137. XM9ULB$3TB72QO=F5R6VYU;5LQ+"(Q(ETL;G5M6S0L(C0B75TL:61;<&QU<ULH
  138. XM*FXJ*6YU;5LQ+"(Q(ETL8V]S6W!L=7-;*"IN*BEI9%MT:6UE<ULH*FXJ*6-H
  139. XM87);8UTL8VAA<EMW75U=+&-O<UMT:6UE<UMI9%MP;'5S6W1I;65S6V-H87);
  140. XM9%TL8VAA<EMW75TL8VAA<EMB971A75U=+'1I;65S6R@J;BHI<W1R:6YG6R@J
  141. XM;BHI(A$K(ETL8V]S6R@J;BHI=&EM97-;*"IN*BEI9%MT:6UE<UMC:&%R6R@J
  142. XM;BHI8UTL8VAA<ELH*FXJ*7==75TL8V]S6R@J;BHI:61;<&QU<ULH*FXJ*71I
  143. XM;65S6R@J;BHI8VAA<ELH*FXJ*61=+&-H87);*"IN*BEW75TL8VAA<ELH*FXJ
  144. XM*6)E=&%=75U=75U=75U=75U=76,!/4W2 B!VR #O^=L KXDW#ML.VP[.   .
  145. XMG0        X =%T           [D#N4             87( -@   #8   #O
  146. XM^=L!1!)N#MP.VPZ?#MX.G0        X 8V@           [@#O,         
  147. XM    8V@      -SWS@#O^=L &P  #MT.VP  #OL.P@        P *6T     
  148. XM -SWS@#O^=L &P                        (,           L8P      
  149. XM RQC  555  %550  JJJ "@!P@^H*BD -QP(^!Q0:&%S92!S:&EF=&5D($AA
  150. XM;FYI;F<@<FEP<&QE(&9U;F-T:6]N<RAT;W I+"!A($AA;@9H72P #E52  55
  151. XM5  "JJH   '"#ZAN*@ Y;FEN9R!W:6YD;W=I;F<@9G5N8W1I;VXH;6ED9&QE
  152. XM*2P@86YD($AA;FYI;F<@<FEP<&QE(&9U;F,&   755  !554  *JJ@!  <(/
  153. XMJ    #ET:6]N<R!M=6QT:7!L:65D(&)Y('1H92!(86YN:6YG('=I;F1O=R!F
  154. XM=6YC=&EO;BAB;W1T;VTI+@L  3U-T@ V    [_G; )V)-P[>#ML.W [.#IT 
  155. XM       . !P(           .X0\              ',@      "4B3< [_G;
  156. XM !L   [?#ML   [B#LX        , ""/      "4B3< [_G; !L         
  157. XM               "#P          9"        -F=  %550 !554  *JJ@ H
  158. XM = /J&)O #D<")D<02!W:6YD('9E;&]C:71Y('9I<W5A;&EZ871I;VX@:7,@
  159. XM8W)E871E9"!B>2!C;VUP;W-I= 9D  Y54@ %550  JJJ   !T ^H(ET .FEN
  160. XM9R!A;B!I;6%G92!O9B!.;W)T:"!!;65R:6-A('5N9&5R(&%N(&EM86=E(&]F
  161. XM('1H92!V96QO8P8 %U50  555  "JJH 0 '0#ZAN*@!%:71Y(&9I96QD(')E
  162. XM;F1E<F5D('5S:6YG('9A<FEA8FQE(&QE;F=T:"!,24,@;W9E<B Q+QL/JAMF
  163. XM&P^H&R!N;VES92X+70     !*1)N ._YVP ;   .X [;   .\P[<        
  164. XM#  H*@     !*1)N ._YVP ;                         @T         
  165. XM "IN       #+&, !554  555  "JJH * ',#ZAN*@ U' D '%=H:71E(&YO
  166. XM:7-E(&-O;G9O;'9E9"!W:71H(&9L=6ED(&1Y;F%M:6-S('9E8W1O<B!=  Y5
  167. XM4@ %550  JJJ   !S ^H=L@ 0&9I96QD('5S:6YG('9A<FEA8FQE(&YO<FUA
  168. XM;&EZ871I;VX@*'1O<"D@=F5R<W5S(&9I>&5D(&YO<FUA;&EZ8082;@ 755  
  169. XM!554  *JJ@!  <P/J&-H  ]T:6]N("AB;W1T;VTI+@L          'F)-P#O
  170. XM^=L )   #N$.VP  #P .W@        P          'F)-P#O^=L )       
  171. XM                  (.                    !     555  %550  JJJ
  172. XM "@!S@^H <( ,QP(EAQ7:&ET92!N;VES92!C;VYV;VQV960@=VET:"!C:&5C
  173. XM:V5R8F]A<F0@=F5C=&]R("P@80 .55( !554  *JJ@   <X/J*JJ $%F:65L
  174. XM9"!U<VEN9R!F:7AE9"!N;W)M86QI>F%T:6]N("AL969T*2P@86YD('1H96X@
  175. XM9W)A9&EE;G0@<VAA9&5D(   %U50  555  "JJH   '.#Z@    \*')I9VAT
  176. XM*2!T;R!G:79E('1H92!A<'!E87)A;F-E(&]F(&$@<F]U9V@@=V]V96X@<W5R
  177. XM9F%C92!T97@&3=( (%5.  555  "JJH 0 '.#Z@.S@ &='5R92X+  "X40  
  178. XMN%$ [HDW (Z/7 [B#N<.WP  #LX        & (DW ._YVP V    -@   ._Y
  179. XMVP%$$FX.XP[;#J(.Y@Z@        #@#YVP          #NP/ 0          
  180. XM            E(DW ._YVP ;   .Y [;   .Y0[;        # !B;P      
  181. XME(DW ._YVP ;                         A           &%T       #
  182. XM8V\ !554  555  "JJH * '1#ZBJJ@ X' D)'%1H92!F:7AE9"!N;W)M86QI
  183. XM>F%T:6]N(&9L=6ED(&1Y;F%M:6-S(&9I96QD(&EM86=E9"!T:  .55( !554
  184. XM  *JJ@ ( =$/J*JJ $=I;B!F:6=U<F4@' C@'#@<".(<(&ES(&UU;'1I<&QI
  185. XM960@8GD@82!C;VQO<B!I;6%G92!O9B!T:&4@;6%G;FET=61E(&]F(&]I<P 7
  186. XM55  !554  *JJ@!  =$/J    !)T:&4@=F5C=&]R(&9I96QD+@L  +A1  "X
  187. XM40#NB3< CH]<#N4.YP[D   .VP        8          #8   ([=L@ [_G;
  188. XM )2)-P[F#ML.XP[G#J         . !P)           .Z0\#            
  189. XM &0@ 3U-T@' ?.T [_G; 0^#$@[G#ML.Y@  #J         . %54        
  190. XM   .Z@\"             &EN -NWZ "#%-0 %B28  D,,P[H#ML.UP[K#JH 
  191. XM       , '-U -NWZ "#%-0 %B28  D,,P                     "!@  
  192. XM        =&D       %T=  &    !@    ,   !H =(/HOG;  4<"-,<"]L 
  193. XM     'F)-P#O^=L &P  #ND.VP  #P,.Y@        P          'F)-P#O
  194. XM^=L &P                        (1          "JJ@      ! '"  55
  195. XM5  %550  JJJ "@!TP^H8V\ .!P)'AQ4:&4@;W)I9VEN86P@<&AO=&\@;VX@
  196. XM=&AE(&QE9G0@<VAO=W,@;F\@;6]T:6]N(&)L=7(&9FD #E52  555  "JJH 
  197. XM  '3#ZAR;0 _<FEN9R!4:&4@<&AO=&\@;VX@=&AE(')I9VAT('5S97,@=F%R
  198. XM:6%B;&4@;&5N9W1H($Q)0R!T;R!M;W1I;VX@/"AR !=54  %550  JJJ   !
  199. XMTP^H87  /6)L=7(@0F]R:7,@665L='-I;M5S('=A=FEN9R!A<FTL('-I;75L
  200. XM871I;F<@82!S;&]W97(@<VAU='1E<B!U "6JJ  %550  JJJ $ !TP^HCUP 
  201. XM!W-P965D+@O.         /2#$@#O^=L &P  #NH.VP  #P(.YP        P 
  202. XM#N8      /2#$@#O^=L &P                        (2          ")
  203. XM-P       P    555  %550  JJJ "@!U0^H8F\ /!P)(1Q!('1H<F5E+61I
  204. XM;65N<VEO;F%L(#4Q,AL/K1LS&P^H&R!E;&5C=')O<W1A=&EC(&9I96QD(&ES
  205. XM(*JJ  Y54@ %550  JJJ   !U0^H9FD /6EM86=E9"!B>2!V;VQU;65T<FEC
  206. XM86QL>2!R87D@=')A8VEN9R!A('1H<F5E+61I;65N<VEO;F%L('-C80;1 !=5
  207. XM4  %550  JJJ $ !U0^HX!P +VQA<B!F:65L9"!P<F]D=6-E9"!U<VEN9R!,
  208. XM24,@;W9E<B!W:&ET92!N;VES92X+='5D  ,'9 !E9R\ TQ^( $4$60[K#MD.
  209. XMZ   #JH        / '9E#LP.UP     !*1)N ._YVP ;   .[ [;   / 0[C
  210. XM        #          !*1)N ._YVP ;                         A, 
  211. XM                   #    !554  555  "JJH ( ';#ZA1@0 R02!P:&]T
  212. XM;V=R87!H("AT;W I(&]F('1H92!&;&%V:6%N($%M<&AI=&AE871E<B!I<R  
  213. XM#E52  555  "JJH   ';#Z@DF !!=7-E9"!T;R!C<F5A=&4@82!V96-T;W(@
  214. XM9FEE;&0N(%1H92!F:65L9"!I<R!T:&5N('5S960@=&\@:6UA9V4@82!I !=5
  215. XM4  %550  JJJ $@!VP^H    .])P86EN=&5DTR!V97)S:6]N(&]F('1H92!S
  216. XM86UE('!H;W1O9W)A<&@@*&)O='1O;2D<"+\<&P^N&RX+[_G;  !NEP      
  217. XM[QUB -AR"0[M#MD/#0  #L(        /    #OL/#0           ._YVP#O
  218. XM^=L.[@[:#M8   [&        #0!I;@              (       8SQU/CQU
  219. XM/CQC/F%V<SQC/F1A=&$\8SYF:6YA;"YP87!E<CQC/F9I9W5R93,N97!S/%4^
  220. XM+VAO;64O=3 V+V-A8G)A;"]A=G,O9&%T82]F:6YA;"YP87!E<B]F:6=U<F4S
  221. XM+F5P<T,@=            ._YVP#O^=L.[P[:#O0   ZR        #0!B;   
  222. XM            ;       8SQU/CQU/CQC/F%V<SQC/F1A=&$\8SYF:6YA;"YP
  223. XM87!E<CQC/F9I9W5R93$N97!S/%4^+VAO;64O=3 V+V-A8G)A;"]A=G,O9&%T
  224. XM82]F:6YA;"YP87!E<B]F:6=U<F4Q+F5P<^<   $]3=( -@   ._YVP"+B3<.
  225. XM\ [A#I,.L@Z/        #@              #O4/#@                  
  226. XMN%$  +A1 .Z)-P#NB3<.\0[G#MD   [+        !@!D:0!E;G,  +A1  "X
  227. XM40#NB3< <X]<#O(.YP[5   .Q         8 :64 9"!I  "X40  N%$ [HDW
  228. XM 2,8DP[S#N<.X   #MP        & 'D@ &]L=0      ](,2 ._YVP 2   .
  229. XM] [;   .[PZR        # !E;@      ](,2 ._YVP 2                
  230. XM         @4          &0@       "960 !554  555  "JJH * &P#ZAI
  231. XM=  V' C/'%1H92!M87!P:6YG(&]F(&$@=F5C=&]R(&]N=&\@82!$1$$@;&EN
  232. XM92!A;F0@:6YP=70@  Y54@ %550  JJJ $ !L ^H#P$ +G!I>&5L(&9I96QD
  233. XM(&=E;F5R871I;F<@82!S:6YG;&4@;W5T<'5T('!I>&5L+@L      'F)-P#O
  234. XM^=L $@  #O4.VP  #PX.\         P 48$      'F)-P#O^=L $@      
  235. XM                  ('          !H90       G,@  555  %550  JJJ
  236. XM "@!M0^H))@ .!P(XQQ#:7)C=6QA<B!A;F0@='5R8G5L96YT(&9L=6ED(&1Y
  237. XM;F%M:6-S('9E8W1O<B!F:65L9',@(&D #E52  555  "JJH 0 &U#ZBJJ@ O
  238. XM:6UA9V5D('5S:6YG($1$02!C;VYV;VQU=&EO;B!O=F5R('=H:71E(&YO:7-E
  239. XM+@MH("@  &Z7     0 O__\ 2   #OL.V@[=#P0.P@  #P0.[0T <@D     
  240. XM   !#@           %D\=3X\=3X\8SYA=G,\8SYD871A/&,^9FEN86PN<&%P
  241. XM97(\8SYH,#$N<',\53XO:&]M92]U,#8O8V%B<F%L+V%V<R]D871A+V9I;F%L
  242. XM+G!A<&5R+V@P,2YP<V$  +A1  "X40#NB3< <X]<#P .YP[A   .W@      
  243. XM  8 86P 879S  "X40  N%$ [HDW 2,8DP\!#N<.[   #N,        &    
  244. XM        N%$  +A1 .Z)-P#NB3</ @[G#NH   [G        !@       &P 
  245. XM +A1  "X40#NB3< <X]<#P,.YP[I   .Y@        8 <&$ 97(\ "^H*P  
  246. XM    +_[% $@   \$#MH.^P\%#L(.^P\%#NT- '9S         0X    N    
  247. XM  !9/'4^/'4^/&,^879S/&,^9&%T83QC/F9I;F%L+G!A<&5R/&,^:# R+G!S
  248. XM/%4^+VAO;64O=3 V+V-A8G)A;"]A=G,O9&%T82]F:6YA;"YP87!E<B]H,#(N
  249. XM<',W %[@B0    ( +__@ $@   \%#MH/! \&#L(/! \&#NT- (DW        
  250. XM 0X           !9/'4^/'4^/&,^879S/&,^9&%T83QC/F9I;F%L+G!A<&5R
  251. XM/&,^:# S+G!S/%4^+VAO;64O=3 V+V-A8G)A;"]A=G,O9&%T82]F:6YA;"YP
  252. XM87!E<B]H,#,N<',  (\V 0    ( +__\ $@   \&#MH/!0\'#L(/!0\'#NT-
  253. XM             0X   !D      !9/'4^/'4^/&,^879S/&,^9&%T83QC/F9I
  254. XM;F%L+G!A<&5R/&,^:# T+G!S/%4^+VAO;64O=3 V+V-A8G)A;"]A=G,O9&%T
  255. XM82]F:6YA;"YP87!E<B]H,#0N<',! +^+^@    $ +___ $@   \'#MH/!@\(
  256. XM#L(/!@\(#NT- &]U         0X    +      !9/'4^/'4^/&,^879S/&,^
  257. XM9&%T83QC/F9I;F%L+G!A<&5R/&,^:# U+G!S/%4^+VAO;64O=3 V+V-A8G)A
  258. XM;"]A=G,O9&%T82]F:6YA;"YP87!E<B]H,#4N<'-4 %[@B0!'JLL +___ $@ 
  259. XM  \(#MH/!P\)#L(/!P\)#NT- &)U         0X   !D      !A/'4^/'4^
  260. XM/&,^879S/&,^9&%T83QC/F9I;F%L+G!A<&5R/&,^:&%N;FEN9RYP<SQ5/B]H
  261. XM;VUE+W4P-B]C86)R86PO879S+V1A=&$O9FEN86PN<&%P97(O:&%N;FEN9RYP
  262. XM<P$  &Z7 )!R"0 O__P 2   #PD.V@\(#PH.P@\(#PH.[0T  0X        !
  263. XM#@   '4      %D\=3X\=3X\8SYA=G,\8SYD871A/&,^9FEN86PN<&%P97(\
  264. XM8SYH,3$N<',\53XO:&]M92]U,#8O8V%B<F%L+V%V<R]D871A+V9I;F%L+G!A
  265. XM<&5R+V@Q,2YP<U$ +Z@K )!R"0 O_^  2   #PH.V@\)#PL.P@\)#PL.[0T 
  266. XMN%$        !#@   ),      %D\=3X\=3X\8SYA=G,\8SYD871A/&,^9FEN
  267. XM86PN<&%P97(\8SYH,3(N<',\53XO:&]M92]U,#8O8V%B<F%L+V%V<R]D871A
  268. XM+V9I;F%L+G!A<&5R+V@Q,BYP<P  7N") )!R"0 O__P 2   #PL.V@\*#PP.
  269. XMP@\*#PP.[0T #OL        !#@           %D\=3X\=3X\8SYA=G,\8SYD
  270. XM871A/&,^9FEN86PN<&%P97(\8SYH,3,N<',\53XO:&]M92]U,#8O8V%B<F%L
  271. XM+V%V<R]D871A+V9I;F%L+G!A<&5R+V@Q,RYP<W( CS8! )!R"0 O_^  2   
  272. XM#PP.V@\+#PT.P@\+#PT.[0T #P0        !#@           %D\=3X\=3X\
  273. XM8SYA=G,\8SYD871A/&,^9FEN86PN<&%P97(\8SYH,30N<',\53XO:&]M92]U
  274. XM,#8O8V%B<F%L+V%V<R]D871A+V9I;F%L+G!A<&5R+V@Q-"YP<W( OXOZ )!R
  275. XM"  O_L4 2   #PT.V@\,#NT.P@\,   .[0T #P4        !#@          
  276. XM %D\=3X\=3X\8SYA=G,\8SYD871A/&,^9FEN86PN<&%P97(\8SYH,34N<',\
  277. XM53XO:&]M92]U,#8O8V%B<F%L+V%V<R]D871A+V9I;F%L+G!A<&5R+V@Q-2YP
  278. XM<W(  +A1  "X40#NB3< <X]<#PX.YP[U   .\         8 #P8 " X1#B$.
  279. XM%  +   "9    Q@   X1      YU      $ /F$ !$QE9G0   )D   #&   
  280. XM#A(     #G8      0#OP  %4FEG:'1E   P-@)D   #&   #A,     #H  
  281. XM     @ N<  )4F5F97)E;F-E5   X(D"9    Q@   X5   .%@Y]        
  282. XM#P<      F0   ,8   .%@X5#A<.CP !     #QU      )D   #&   #A<.
  283. XM%@X8#I0  @    !P80     "9    Q@   X8#A<.&0Z7  ,     ,#8     
  284. XM F0   ,8   .&0X8#AH.F@ $     &5R      )D   #&   #AH.&0X;#IT 
  285. XM!0             "9    Q@   X;#AH.' Z@  8      0X      F0   ,8
  286. XM   .' X;   .HP '     &%T      \]#XT/10 V                    
  287. XM   $     0    #F9@ !     1F9                      \]#Z$   0#
  288. XM%     (# P4!  #@  P                                    " 0 &
  289. XM!0$             "$%B<W1R86-T                                
  290. XM          (     YF8  0    $9F0                     //@     %
  291. XM 00    " P,% 0  B0 ,                                     @$ 
  292. XM!@4!     0X       A#96QL0F]D>0                              
  293. XM           "     .9F  $    !&9D       (          @  #S\     
  294. XM!0,      @,#!0$  3, #                                     (!
  295. XM  8% 0    \$       +0V5L;$AE861I;F=9                    "0  
  296. XM                  $     YF8  0    $9F0                     /
  297. XM0 ^B   $! 0    " P,% @  <@ )                                
  298. XM     @$ !@4!    #NT       1";V1Y                            
  299. XM   .    !0    (     YF8  0    $9F0                     /00  
  300. XM   % 50    " P,% 0  80 .                                    
  301. XM @$ !@D!   .#N<       =(96%D:6YG              1";V1Y !(    2
  302. XM    $@               @    #F9@ !     1F9                    
  303. XM  ]"      4!!     (# P4!      H                             
  304. XM       " 0 &!0$  &0         #51A8FQE1F]O=&YO=&4'            
  305. XM                                 @    #F9@ !     1F9        
  306. XM              ]#    5 4# 0    (# P4!   V  P                 
  307. XM                   " 0 &"0$  !@         "E1A8FQE5&ET;&4     
  308. XM    #E0Z5&%B;&4@/&XK/CH@                          0    !    
  309. XM .9F  $    !&9D                     #T0     !0,4     @,#!0$ 
  310. XM &8 $@                                    ,!  8) 0   P,%    
  311. XM   %5&ET;&4                                              @  
  312. XM  #F9@ !     1F9       "          (   ]&#Z$   0#      (# P4!
  313. XM  %F  P                                    " 0 &!0$   ,#!0  
  314. XM    "T-E;&Q(96%D:6YG                                        
  315. XM   "     .9F  $    !&9D                     #T</H0  ! $$    
  316. XM @,#!0$     #                                     (!  8% 0  
  317. XM  $S       (0V5L;$)O9'D                                     
  318. XM     @    #F9@ !     1F9                      ](#Z0 5 0# 0  
  319. XM  (# P4!  !F  P                                    " 0 &"0$ 
  320. XM  ,#!0      "E1A8FQE5&ET;&4         #E0Z5&%B;&4@/&XK/CH@    
  321. XM                      0    !     .9F  $    !&9D             
  322. XM        #TD/K   ! ,4     @,#!0$     #                       
  323. XM              (!  8%!0             (06)S=')A8W0             
  324. XM                  ,    #     0    #F9@ !     1F9            
  325. XM          ]*#Z(  00$%0    (# P4!      D                     
  326. XM               " 0 &!0$             "$5Q=6%T:6]N          <H
  327. XM/&X].#XI           )                     0    #F9@ !     1F9
  328. XM                      ],      4$!     (# P4"      D         
  329. XM                           " 0 &!0$             !$)O9'D     
  330. XM                          <          0    #F9@ !     1F9    
  331. XM                  ]-    4P4$%0    (# P4!      D             
  332. XM                       # 0 &"0$             "%-E8W1I;VXS    
  333. XM     !!3.CQN/BX\;CXN/&XK/B @                    !P         !
  334. XM     .9F  $    !&9D                     #TX/JP  ! $4     @,#
  335. XM!0$     "0                                    ,!  8) 0      
  336. XM       &2&5A9&5R              $                             
  337. XM 0    #F9@ !     1F9                      ]/#[D   0$5     (#
  338. XM P4!  "9  @                                    " 0 &!0$   P 
  339. XM        "$9O;W1N;W1E                               '        
  340. XM  $     YF8  0    $9F0                     /4      % 10    "
  341. XM P,% 0  F0 )                                     P$ !@D!   ,
  342. XM          9(96%D97(                                 !P      
  343. XM   !     .9F  $    !&9D                     #U$   !3!0$5    
  344. XM @,#!0$     "0                                    ,!  8) 0  
  345. XM  0#       (4V5C=&EO;C(         #%,Z/&X^+CQN*SX@(           
  346. XM          ,    #     0    #F9@ !     1F9                    
  347. XM  ]2#Z(  00$%0    (# P4!      D                             
  348. XM       " 0 &!0$             "$5Q=6%T:6]N          <H/&X],3XI
  349. XM @                               0    #F9@ !     1F9        
  350. XM              ]4#Z@ 1@0!%0    (# P4!      @                 
  351. XM                   % 0 &!@$             #D9I9W5R92U#87!T:6]N
  352. XM          Y&.D9I9W5R92 \-SXZ(                               
  353. XM 0    #F9@ !     1F9                      ]5      4$5     (#
  354. XM P4!      @                                    " 0 &!0$     
  355. XM        "$9O;W1N;W1E                               #     P  
  356. XM  $     YF8  0    $9F0                     /5@    $%!!4    "
  357. XM P,% 0     )                                     @$ !@4!    
  358. XM          A%<75A=&EO;@         '*#QN/3$^*0                  
  359. XM              $     YF8  0    $9F0                     /5P^H
  360. XM $8$ 14    " P,% 0     (                                    
  361. XM!0$ !@8!              Y&:6=U<F4M0V%P=&EO;@         .1CI&:6=U
  362. XM<F4@/#@^.B                                $     YF8  0    $9
  363. XMF0                     /60^H $8$ 14    " P,% 0  !0 (        
  364. XM                            !0$ !@8!   ) 0        Y&:6=U<F4M
  365. XM0V%P=&EO;@         .1CI&:6=U<F4@/#$^.B                     #
  366. XM     P    $     YF8  0    $9F0                     /6@^B  $$
  367. XM!!4    " P,% 0     )                                     @$ 
  368. XM!@4!   R          A%<75A=&EO;@         '*#QN/3(^*0          
  369. XM           '          $     YF8  0    $9F0                  
  370. XM   /6P^K %,$ 14    " P,% 0     )                            
  371. XM         P$ !@D!  !"          A396-T:6]N,0         )4SH\;BL^
  372. XM+B @                         P    ,    !     .9F  $    !&9D 
  373. XM                    #UP/H@ !! 05     @,#!0$     "0          
  374. XM                          (!  8% 0  0V%P       (17%U871I;VX 
  375. XM        !R@\;CTW/BD                      P    ,    !     .9F
  376. XM  $    !&9D                     #UT/H@ !! 05     @,#!0$     
  377. XM"0                                    (!  8% 0             (
  378. XM17%U871I;VX         !R@\;CTV/BD                      P    , 
  379. XM   !     .9F  $    !&9D                     #UX/H@ !! 05    
  380. XM @,#!0$     "0                                    (!  8% 0  
  381. XM871I       (17%U871I;VX         !R@\;CTT/BD                 
  382. XM     P    ,    !     .9F  $    !&9D                     #U\/
  383. XMH@ !! 05     @,#!0$     "0                                  
  384. XM  (!  8% 0  =7)E       (17%U871I;VX         !R@\;CTU/BD@    
  385. XM                 P    ,    !     .9F  $    !&9D             
  386. XM        #V /H@ !! 05     @,#!0$   4 "0                      
  387. XM              (!  8% 0  "0$        (17%U871I;VX         !R@\
  388. XM;CTS/BDZ    $@                    (    !     .9F  $    !&9D 
  389. XM                    #V(/H@!2! 0%     @,#!0,   0 "0          
  390. XM                          (!  8% 0           0 2   !        
  391. XM"5)E9F5R96YC96X         "%(Z/&XK/BX(                        
  392. XM       !     .9F  $    !&9D                     #V,   !&!0$5
  393. XM     @,#!0$     "                                     4!  8&
  394. XM 0  0@         .1FEG=7)E+4-A<'1I;VX         #D8Z1FEG=7)E(#PS
  395. XM/CH@                    !P         !     .9F  $    !&9D     
  396. XM                #V0   !3!0$5     @,#!0$     "0              
  397. XM                      ,!  8) 0  "$5Q       (4V5C=&EO;C$     
  398. XM    "5,Z/&XK/BX@(                              $     0    # 
  399. XM   !     4                        ]E#Z(   0#%     (" @4!    
  400. XM  D                                    " 0 &!0$  &%T:0    ( 
  401. XMZ@   @ H/  !(#X!U    P     !(   "$%B<W1R86-T                
  402. XM                    !     $     YF8  0    $9F0              
  403. XM       /9@^B   $ Q0    " P,% 0     )                        
  404. XM             @$ !@4!   ]-#X       A!8G-T<F%C=               
  405. XM                !P         !     .9F  $    !&9D             
  406. XM        #V</JP!3! $5     @,#!0$     "0                      
  407. XM              ,!  8) 0             (4V5C=&EO;C(         #%,Z
  408. XM/&X^+CQN*SX@(                     <          0    #F9@ !    
  409. XM 1F9                      ]I#ZL 4P0$%0    (# P4!      D     
  410. XM                               # 0 &"0$   <H/       "%-E8W1I
  411. XM;VXS         !!3.CQN/BX\;CXN/&XK/B @    !@                  
  412. XM       !     .9F  $    !&9D                     #VH/N0  ! %0
  413. XM     @,#!0$     "                                     (!  8%
  414. XM 0             (1F]O=&YO=&4 "$9O;W1N;W1E                    
  415. XM              0    !     .9F  $    !&9D                     
  416. XM#VP/IP  ! ,4     @,#!0$     $@                              
  417. XM      (!  8% 0             (06)S=')A8W0                     
  418. XM                     0    #F9@ !     1F9                    
  419. XM  ]N#Z@ 1@0!%0    (# P4!      @                             
  420. XM       % 0 &!@$   A%<0      #D9I9W5R92U#87!T:6]N         !!&
  421. XM.D9I9W5R92 \;CTR/CH@                            !     $     
  422. XMYF8  0    $9F0                     /< _!   $ Q0    " P,% 0  
  423. XM   2                                     P$ !@D!    *#P     
  424. XM  54:71L90              <F$ $@                    (    !    
  425. XM .9F  $    !&9D                     #W,   !2!00%     @,#!0, 
  426. XM  , "0                                    (!  8% 0          
  427. XM 0 2   ! #0^    "5)E9F5R96YC970         "%(Z/&XK/BX(        
  428. XM                       !     .9F  $    !&9D                 
  429. XM    #W0/J !&! $5     @,#!0$     "                           
  430. XM          4!  8& 0  "%-E       .1FEG=7)E+4-A<'1I;VX         
  431. XM$$8Z1FEG=7)E(#QN/3,^.B                                   0  
  432. XM  #F9@ !     1F9                      ]U#Z@ 1@0!%0    (# P4!
  433. XM      @                                    % 0 &!@$         
  434. XM    #D9I9W5R92U#87!T:6]N         !!&.D9I9W5R92 \;CTT/CH@  #F
  435. XM9@                            $     YF8  0    $9F0          
  436. XM           /=@^H $8$ 14    " P,% 0     (                    
  437. XM                !0$ !@8!              Y&:6=U<F4M0V%P=&EO;@  
  438. XM       01CI&:6=U<F4@/&X]-3XZ(                               
  439. XM   !     .9F  $    !&9D                     #W</J !&! $5    
  440. XM @,#!0$     "                                     4!  8& 0  
  441. XM           .1FEG=7)E+4-A<'1I;VX         $$8Z1FEG=7)E(#QN/38^
  442. XM.B    ,%                             0    #F9@ !     1F9    
  443. XM                  ]X#Z@ 1@0!%0    (# P4!   0  @             
  444. XM                       % 0 &!@$   0         #D9I9W5R92U#87!T
  445. XM:6]N         !!&.D9I9W5R92 \;CTY/CH@   #!0                  
  446. XM          $     YF8  0    $9F0                     />0^H $8$
  447. XM 14    " P,% 0     (                                    !0$ 
  448. XM!@8!              Y&:6=U<F4M0V%P=&EO;@         11CI&:6=U<F4@
  449. XM/&X],3$^.B                                 !     .9F  $    !
  450. XM&9D                     #WH/J !&! $5     @,#!0$     "       
  451. XM                              4!  8& 0             .1FEG=7)E
  452. XM+4-A<'1I;VX         $48Z1FEG=7)E(#QN/3$P/CH@                
  453. XM                 0    #F9@ !     1F9                      ][
  454. XM#Z@ 1@0!%0    (# P4!      @                                 
  455. XM   % 0 &!@$             #D9I9W5R92U#87!T:6]N         !%&.D9I
  456. XM9W5R92 \;CTQ,CXZ(                                 $     YF8 
  457. XM 0    $9F0                     /? ^H $8$ 14    " P,% 0     (
  458. XM                                    !0$ !@8!              Y&
  459. XM:6=U<F4M0V%P=&EO;@         11CI&:6=U<F4@/&X],30^.B          
  460. XM                       !     .9F  $    !&9D                 
  461. XM    #WT/J !&! $5     @,#!0$     "                           
  462. XM          4!  8& 0             .1FEG=7)E+4-A<'1I;VX         
  463. XM$48Z1FEG=7)E(#QN/3$S/CH@    #Z$/T0^R !P #                   
  464. XM     'K95M,/H0    (!  8% 0 $.D9I  !R90 )                    
  465. XM    >ME6TP^B     @$ !@4!  0 YF8       D                  /_,
  466. XMX       #Z,    " 0 &!0$ !0  $  (16UP:&%S:7,    ,            
  467. XM            %>[C^ ^D     @$ !@D!  1U<F4  &%P  D             
  468. XM     /_,X       #Z4    " 0 &!0$ !0     )4W5B<V-R:7!T   )    
  469. XM              #_S.        ^F     @$ !@4!  5&! $ "U-U<&5R<V-R
  470. XM:7!T"    !(                       !ZV5;3#Z<    " 0 &!0$ !   
  471. XM    1FD "                        'NM/HH/J     4!  8& 0 $    
  472. XM       )                        %>[C^ ^I     @$ !@D!  0     
  473. XM      @                        $K;.B#ZH    % 0 &!@0 !       
  474. XM    "0                        >I%:T/JP    ,!  8) 0 $:6]N    
  475. XM   ,                        :+"$L@^L     @$ !@4%  0         
  476. XM  @                       ![K3Z*#ZT    % 0@&!@$ !$8$ 0      
  477. XM"                        /(3>3P/K@    4!  8(!  $!0$    !   )
  478. XM                        :+"$L@^O     @$ !@4%  1U<F4 !DET86QI
  479. XM8P )                        %>[C^ ^P     @$ !@D!  0!&9D !$)O
  480. XM;&0    )                        %>[C^ ^Q     @$(!@D!  0     
  481. XM!$)O;&0    )                        6@1L60^X     @$ !@D%  1I
  482. XM;VX       @                       !ZV5;3#[D    " 0 &!0$ !   
  483. XM        "0                       &BPA+(/N@    (!  8%!0 $1@0!
  484. XM       )                        !$8I>P^[    ! $ !@4!  0% 0  
  485. XM  $   D                       !HL(2R#[T    " 1 &!04 !'5R90 &
  486. XM271A;&EC  D                        $1BE[#[\    $ 0 &!0$ !*$ 
  487. XM   &271A;&EC  @                       !HL(2R#\     " 0 &!04 
  488. XM! (!   &271A;&EC !(                        'J16M#\$    # 0 &
  489. XM"0$ ! 4!     !  "0                       &BPA+(/QP    (!$ 8%
  490. XM!0 $I      !   )                        >ME6TP_(     @$ !@4!
  491. XM  2E    !DET86QI8P )                        >ME6TP_)     @$(
  492. XM!@4!  0     !DET86QI8P[9#OD.X  +  $   [9    #P              
  493. XM 0  #MH #P /                @  .VP /  <,   $ @( $%H'___8.P[<
  494. XM    #P               0  #MT    /      ("        @  .W@    \ 
  495. XM     @(       $   [?  \ !PP   0" @ 06@<  0  #N$ #P '      ("
  496. XM           .X@                    $   [C  \ #P     " @      
  497. XM !F9#N< #@ '#   ! (" !!:!P!E '4 :@ %  "         90     !   $
  498. XM5&AI;@    (         9@     !   &365D:75M  "    "    9P     "
  499. XM   &1&]U8FQE  ,         :      !   %5&AI8VL   !         :0  
  500. XM   !   )5F5R>2!4:&EN             ,D V0#)  ( #     P         
  501. XM      8    &    !@    8    $       % ,H  @ /  \ #P!E &4 90!E
  502. XM &4 90!E &4 90!G#T@%   !     @($     0@   $!!0              
  503. XM  !(   /1@]'#T8       $           !(   /1@]'#T8 ;P    (     
  504. XM      !(   /1@]'#T8       ,           !(   /1@]'#T8       0 
  505. XM          !(   /1@]'#T8 TP (1F]R;6%T($$!   ,    #           
  506. XM    !@    8    &    !@    0       4 RP "  \ #P /    90   &8 
  507. XM90!E &4 90!E &8/2 4  0$    " @0    !"    0$%                
  508. XM $@   ]&#T</1@       0           $@   ]&#T</1@!T     @      
  509. XM     $@   ]&#T</1@       P           $@   ]&#T</1@!T    !   
  510. XM         $@   ]&#T</1@    A&;W)M870@0@                     !
  511. XM+0$] 2T  0       0$N  $" 0    =#;VUM96YT P           !$ !  "
  512. XM  54:6UE<P   P )2&5L=F5T:6-AL  $  93>6UB;VP !0 *079A;G1'87)D
  513. XM90 )  $ !@ '4F5G=6QA<@  #  $  4 !U)E9W5L87(%  8 !$)O;VL "  (
  514. XM1&5M:4)O;&0 "0 $0F]L9  %  ,  0 '4F5G=6QA<@( !  '3V)L:7%U908 
  515. XM!0 &271A;&EC                                                
  516. XM                                                            
  517. XM                                                            
  518. XM                                                            
  519. XM                                                            
  520. XM                                                            
  521. XM                                                            
  522. XM                                                            
  523. XM                                                            
  524. XM                                                            
  525. XM                                                            
  526. XM                                                            
  527. XM                                                            
  528. XM                                                            
  529. XM                                                            
  530. XM                                                            
  531. XM                                                            
  532. XM                                                            
  533. XM                                                            
  534. XM                                                            
  535. XM                                                            
  536. XM                                                            
  537. XM                                                            
  538. XM                                                            
  539. XM                                                            
  540. XM                                                            
  541. XM                                                            
  542. XM                                                            
  543. XM                                                       ! @,$
  544. XMM5%@52EG987Y'1P$1Y_3-=H)S0OA6I*(&46.IKGCJ/"[7]7BT-P9%V_=M%P<
  545. XMV+F[")JOI/F?(<'IC$1ZL*0K7J).L>L%#K6<#.-GMBY@L-PL/C#A]W@P_<&_
  546. XM9$^9;M967O-,+M3",!7+D)5P>%6?'(4DVZG175>EK,ZRW O_KHTQ\#PUF!BP
  547. XM6?!@7K$:8^+I'UGDT01[H2:Z%#<E+5TG.X:5R,&-*7C63A_AM;+D#=VA DJ)
  548. XMD]:[-J8D:#!<L*[V:OP?)D;0Z__TBI2YJQ&M7HM$.NX#D8%/"5X+K?Q0CFD8
  549. XM/X=RO2"^IM&&1]/F);6]9$(1\,U<4HJE*R=]7(S*13GWT2YC;6RN9ZU"K*9J
  550. XMG4:N[T+G\Q'8-C*JKIB0S?S TMAR^Y7]"C#GM5LN]J%B9"EP?[2R(#OVK3W_
  551. XM=0+5$T:[MW4[8%^% 1(2[*>L#Z&B.K^4#"#Q8=Q*#D170MC6ZDB7>UWVLF8 
  552. XMLF0SA(^P **1GN)D1_D>N3FQ^/I4E"2\=2?:7XO'#/E%LU_ON*&>$CI_-N< 
  553. X#    
  554. X 
  555. Xend
  556. END_OF_FILE
  557.   if test 32376 -ne `wc -c <'lic.1.2/doc/siggraph93/p.frame.UU.C'`; then
  558.     echo shar: \"'lic.1.2/doc/siggraph93/p.frame.UU.C'\" unpacked with wrong size!
  559.   elif test -f 'lic.1.2/doc/siggraph93/p.frame.UU.A' && test -f 'lic.1.2/doc/siggraph93/p.frame.UU.B' ; then
  560.     echo shar: Combining  \"'paper.frame.UU'\" \(143959 characters\)
  561.     cat 'lic.1.2/doc/siggraph93/p.frame.UU.A' 'lic.1.2/doc/siggraph93/p.frame.UU.B' 'lic.1.2/doc/siggraph93/p.frame.UU.C' > 'paper.frame.UU'
  562.     if test 143959 -ne `wc -c <'paper.frame.UU'`; then
  563.       echo shar: \"'paper.frame.UU'\" combined with wrong size!
  564.     else 
  565.       rm lic.1.2/doc/siggraph93/p.frame.UU.A lic.1.2/doc/siggraph93/p.frame.UU.B lic.1.2/doc/siggraph93/p.frame.UU.C
  566.       echo shar: Uudecoding \"'lic.1.2/doc/siggraph93/paper.frame'\" \(104448 characters\)
  567.       cat paper.frame.UU | uudecode
  568.       if test 104448 -ne `wc -c <'lic.1.2/doc/siggraph93/paper.frame'`; then
  569.         echo shar: \"'lic.1.2/doc/siggraph93/paper.frame'\" uudecoded with wrong size!
  570.       else
  571.         rm paper.frame.UU
  572.       fi 
  573.       # end of 'paper.frame.UU'
  574.     fi
  575.   fi
  576.   # end of 'lic.1.2/doc/siggraph93/p.frame.UU.B'
  577. fi
  578. if test -f 'lic.1.2/liblic/Filters.c' -a "${1}" != "-c" ; then 
  579.   echo shar: Will not clobber existing file \"'lic.1.2/liblic/Filters.c'\"
  580. else
  581.   echo shar: Extracting \"'lic.1.2/liblic/Filters.c'\" \(6103 characters\)
  582.   sed "s/^X//" >'lic.1.2/liblic/Filters.c' <<'END_OF_FILE'
  583. X/*
  584. X * $Header: /d/sisal/a/casey/tmp/lic/liblic/RCS/Filters.c,v 1.12 1993/08/10 23:35:32 casey Exp $
  585. X */
  586. X
  587. X/*
  588. X * Copyright (c) 1993 The Regents of the University of California.
  589. X * All rights reserved.
  590. X *
  591. X * Redistribution and use in source and binary forms, with or without
  592. X * modification, are permitted provided that the following conditions
  593. X * are met:
  594. X * 1. Redistributions of source code must retain the above copyright
  595. X *    notice, this list of conditions and the following disclaimer.
  596. X * 2. Redistributions in binary form must reproduce the above copyright
  597. X *    notice, this list of conditions and the following disclaimer in the
  598. X *    documentation and/or other materials provided with the distribution.
  599. X * 3. All advertising materials mentioning features or use of this software
  600. X *    must display the following acknowledgement:
  601. X *    This product includes software developed by the University of
  602. X *    California, Lawrence Livermore National Laboratory and its
  603. X *    contributors.
  604. X * 4. Neither the name of the University nor the names of its contributors
  605. X *    may be used to endorse or promote products derived from this software
  606. X *    without specific prior written permission.
  607. X *
  608. X * THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND
  609. X * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
  610. X * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
  611. X * ARE DISCLAIMED.  IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE
  612. X * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
  613. X * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
  614. X * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
  615. X * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
  616. X * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
  617. X * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
  618. X * SUCH DAMAGE.
  619. X */
  620. X
  621. X#ifndef lint
  622. X    static char rcsid[] = "$Header: /d/sisal/a/casey/tmp/lic/liblic/RCS/Filters.c,v 1.12 1993/08/10 23:35:32 casey Exp $";
  623. X    static char copyright[] =
  624. X    "Copyright (c) 1993 The Regents of the University of California.\n"
  625. X    "All rights reserved.\n";
  626. X#endif
  627. X
  628. X
  629. X#include "liblic.h"
  630. X
  631. X
  632. X/*
  633. X *    Build integral tables from filter integrator.
  634. X *    =============================================
  635. X */
  636. X
  637. Xvoid
  638. XLIC_BuildIntegralTables(LIC *This)
  639. X    /*
  640. X     * Build integral tables for LIC instance if needed.
  641. X     */
  642. X{
  643. X    int speed, length;
  644. X
  645. X    if (!This->NeedIntegration)
  646. X    return;
  647. X    This->NeedIntegration = FALSE;
  648. X
  649. X    /*
  650. X     * Build the convolution integral tables.  Build LIC_INTEGRAL_SPEEDS
  651. X     * different versions of them with speed varying between 0 and
  652. X     * LIC_INTEGRAL_SPEEDS.
  653. X     */
  654. X    for (speed = 0; speed < LIC_INTEGRAL_SPEEDS; speed++)
  655. X    {
  656. X    double s;
  657. X
  658. X    for (length = 0; length < LIC_INTEGRAL_LEN; length++)
  659. X    {
  660. X        s = LIC_Length(This) * length / (double)LIC_INTEGRAL_LEN;
  661. X        This->NegIntegralTable[speed][length] = This->Filter(This,  -s, 0.0, speed);
  662. X        This->PosIntegralTable[speed][length] = This->Filter(This, 0.0,   s, speed);
  663. X    }
  664. X    }
  665. X}
  666. X
  667. X
  668. Xdouble
  669. XLIC_Box(LIC *This, double a, double b, int speed)
  670. X    /*
  671. X     * Compute the integral bewteen points a and b, a <= b, of a box
  672. X     * filter:
  673. X     *
  674. X     *    k(s) = 1
  675. X     */
  676. X{
  677. X    return(b - a);
  678. X}
  679. X
  680. X
  681. Xdouble
  682. XLIC_Ripple(LIC *This, double a, double b, int speed)
  683. X    /*
  684. X     * Compute the integral bewteen points a and b, a <= b, of a phase shifted
  685. X     * Hanning low pass filter multiplied by a Hanning window function:
  686. X     *
  687. X     *    k(s) = (cos(d*s + phase) + 1)/2 * (cos(c*s) + 1)/2
  688. X     *
  689. X     * We need to scale the filter to the length of the filter support.  This
  690. X     * requires scaling each occurance of s by Pi/L.  Things fall out in the
  691. X     * integral which allow us to slip the scaling factors into c and d.
  692. X     */
  693. X{
  694. X    REGISTER double m;        /* speed multiple */
  695. X    REGISTER double p;        /* phase of the Hanning ripple function */
  696. X    REGISTER double f;        /* frequency of the Hanning ripple function */
  697. X    REGISTER double c;        /* Hanning window function dilation constant */
  698. X    REGISTER double d;        /* Hanning ripple filter dilation constant */
  699. X    REGISTER double w;        /* integral value */
  700. X    
  701. X    if (LIC_Length(This) == 0 || a == b)
  702. X    return(0.0);
  703. X    
  704. X    /*
  705. X     * Implement ``speed'' variations by scaling the frequency inversely
  706. X     * to speed.  This causes low magnitude vectors to be convolved using
  707. X     * a higher frequency filter which will yield less distinct output
  708. X     * features.  As the frequency increases the filter begins to approximate
  709. X     * the Hanning window filter.
  710. X     */
  711. X
  712. X    p = LIC_Phase(This);
  713. X    f = LIC_Frequency(This) * (6.0 - 5.0*(double)speed/(LIC_INTEGRAL_SPEEDS - 1));
  714. X    c = 1.0/LIC_Length(This) * M_PI;
  715. X    d = 1.0/LIC_Length(This) * M_PI * f;
  716. X    
  717. X    w =   b                          - a
  718. X    + sin(b*c)/c                 - sin(a*c)/c
  719. X    + sin(b*d + p)/d             - sin(a*d + p)/d
  720. X    + sin(b*(c+d) + p)/(2*(c+d)) - sin(a*(c+d) + p)/(2*(c+d));
  721. X    if (c != d)
  722. X    w +=  sin(b*(c-d) - p)/(2*(c-d))
  723. X        - sin(a*(c-d) - p)/(2*(c-d));
  724. X    w /= 4;
  725. X    return(w);
  726. X}
  727. X
  728. X
  729. Xdouble
  730. XLIC_Ramp(LIC *This, double a, double b, int speed)
  731. X    /*
  732. X     * Compute the integral bewteen points a and b, a <= b, of a biased
  733. X     * triangle filter:
  734. X     *
  735. X     *    k(s) = (s + L) / (2 * L)
  736. X     */
  737. X{
  738. X    if (LIC_Length(This) == 0)
  739. X    return(0.0);
  740. X    else
  741. X    return((a*a - b*b)/(4*LIC_Length(This)) + (a - b)/2);
  742. X}
  743. X
  744. X
  745. Xdouble
  746. XLIC_Select(LIC *This, double a, double b, int speed)
  747. X    /*
  748. X     * Compute the integral bewteen points a and b, a <= b, of a 
  749. X     * narrow Gaussian filter which sits at L-2 along the local
  750. X     * streamline.  This filter's delta function like shape and
  751. X     * it's location near the end of the local streamline means that
  752. X     * the filter acts as a pixel selection filter. Hence it's name.
  753. X     *
  754. X     *    k(s) = 1/2*pi * e^(x - L + 2)^2 
  755. X     */
  756. X{
  757. X    if (LIC_Length(This) == 0)
  758. X    return(0.0);
  759. X    else 
  760. X    return(  erf((b - LIC_Length(This) + 2) / sqrt(2.0))
  761. X           - erf((a - LIC_Length(This) + 2) / sqrt(2.0)));
  762. X}
  763. END_OF_FILE
  764.   if test 6103 -ne `wc -c <'lic.1.2/liblic/Filters.c'`; then
  765.     echo shar: \"'lic.1.2/liblic/Filters.c'\" unpacked with wrong size!
  766.   fi
  767.   # end of 'lic.1.2/liblic/Filters.c'
  768. fi
  769. if test -f 'lic.1.2/lic/lic.c' -a "${1}" != "-c" ; then 
  770.   echo shar: Will not clobber existing file \"'lic.1.2/lic/lic.c'\"
  771. else
  772.   echo shar: Extracting \"'lic.1.2/lic/lic.c'\" \(22645 characters\)
  773.   sed "s/^X//" >'lic.1.2/lic/lic.c' <<'END_OF_FILE'
  774. X/*
  775. X * $Header: /d/sisal/a/casey/tmp/lic/lic/RCS/lic.c,v 1.24 1993/08/11 17:25:35 casey Exp $
  776. X */
  777. X
  778. X/*
  779. X * Copyright (c) 1993 The Regents of the University of California.
  780. X * All rights reserved.
  781. X *
  782. X * Redistribution and use in source and binary forms, with or without
  783. X * modification, are permitted provided that the following conditions
  784. X * are met:
  785. X * 1. Redistributions of source code must retain the above copyright
  786. X *    notice, this list of conditions and the following disclaimer.
  787. X * 2. Redistributions in binary form must reproduce the above copyright
  788. X *    notice, this list of conditions and the following disclaimer in the
  789. X *    documentation and/or other materials provided with the distribution.
  790. X * 3. All advertising materials mentioning features or use of this software
  791. X *    must display the following acknowledgement:
  792. X *    This product includes software developed by the University of
  793. X *    California, Lawrence Livermore National Laboratory and its
  794. X *    contributors.
  795. X * 4. Neither the name of the University nor the names of its contributors
  796. X *    may be used to endorse or promote products derived from this software
  797. X *    without specific prior written permission.
  798. X *
  799. X * THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND
  800. X * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
  801. X * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
  802. X * ARE DISCLAIMED.  IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE
  803. X * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
  804. X * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
  805. X * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
  806. X * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
  807. X * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
  808. X * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
  809. X * SUCH DAMAGE.
  810. X */
  811. X
  812. X#ifndef lint
  813. X    static char rcsid[] = "$Header: /d/sisal/a/casey/tmp/lic/lic/RCS/lic.c,v 1.24 1993/08/11 17:25:35 casey Exp $";
  814. X    static char copyright[] =
  815. X    "Copyright (c) 1993 The Regents of the University of California.\n"
  816. X    "All rights reserved.\n";
  817. X#endif
  818. X
  819. X
  820. X/*
  821. X * UNIX command line interface to the Line Integral Convolution library.
  822. X */
  823. X
  824. X
  825. X#include <stdlib.h>
  826. X#include <unistd.h>
  827. X#include <errno.h>
  828. X#include <string.h>
  829. X#include <stdio.h>
  830. X#include <time.h>
  831. X#include <sys/times.h>
  832. X#include <sys/types.h>
  833. X
  834. X#include <fcntl.h>
  835. X#include <math.h>
  836. X#include <sys/stat.h>
  837. X
  838. X#if defined(HAS_MMAP)
  839. X#   include <sys/mman.h>
  840. X#endif
  841. X
  842. X#include <lic.h>
  843. X
  844. X
  845. X#ifdef NEED_EXIT_CODES
  846. X    /*
  847. X     * EXIT_SUCCESS and EXIT_FAILURE are supposed to be defined in <stdlib.h>
  848. X     * according to the ANSI C X3.159-1989 specification, but Sun OS 4.1.1
  849. X     * fails to define them ...
  850. X     */
  851. X#   define EXIT_SUCCESS    0
  852. X#   define EXIT_FAILURE    1
  853. X#endif
  854. X
  855. X
  856. X#ifdef NEED_STRERROR
  857. X    /*
  858. X     * strerror is supposed to be defined in <string.h> and supplied in the
  859. X     * standard C library according to the ANSI C X3.159-1989 specification,
  860. X     * but Sun OS 4.1.1 fails to define or supply it ...  Unfortunately the
  861. X     * only way we can control this is with an externally supplied define.
  862. X     */
  863. X    extern int      errno;        /* system error number */
  864. X    extern char     *sys_errlist[];    /* system error messages */
  865. X    extern int      sys_nerr;        /* number of entries in sys_errlist */
  866. X
  867. X    static char *
  868. X    strerror(int err)
  869. X    {
  870. X    if (err < 0 || err >= sys_nerr) {
  871. X        static char msg[100];
  872. X
  873. X        sprintf(msg, "system error number %d", err);
  874. X        return(msg);
  875. X    }
  876. X    return(sys_errlist[err]);
  877. X    }
  878. X#endif
  879. X
  880. X
  881. X/*
  882. X *    Arguments
  883. X *    =========
  884. X */
  885. Xstatic char *usage =
  886. X"usage: %s options in-image in-vector-field out-image\n"
  887. X"  -x x-extent            - *required*                  (no default)\n"
  888. X"  -y y-extent            - *required*                  (no default)\n"
  889. X"  -z z-extent            -                             (default 1)\n"
  890. X"  -i i-extent            - in-image X extent           (default x-extent)\n"
  891. X"  -j j-extent            - in-image Y extent           (default y-extent)\n"
  892. X"  -k k-extent            - in-image Z extent           (default z-extent)\n"
  893. X"  -f filter              - 'box', 'ripple', 'ramp' or  (default 'box')\n"
  894. X"                           'select'\n"
  895. X"  -n normalization       - 'fixed' or 'variable'       (default 'variable')\n"
  896. X"  -N                     - input vector field already normalized\n"
  897. X"  -l filter-length       -                             (default 10.0)\n"
  898. X"  -d filter-frequency    -                             (default 3.0)\n"
  899. X"  -p filter-phase        -                             (default 0.0)\n"
  900. X"  -L                     - vary filter-length based on vector magnitude\n"
  901. X"  -S                     - vary filter-phase based on vector magnitude\n"
  902. X"  -r default-red         - for zero vectors            (default -1)\n"
  903. X"  -g default-green       - for zero vectors            (default -1)\n"
  904. X"  -b default-blue        - for zero vectors            (default -1)\n"
  905. X"  -a default-alpha       - for zero vectors            (default -1)\n"
  906. X"  -v                     - verbose progress and timing information\n"
  907. X"  -V                     - version information\n"
  908. X    ;
  909. X
  910. Xstatic char *myname;            /* argv[0]: name we were invoked by */
  911. Xstatic int    size_x = -1;        /* -x: X extent */
  912. Xstatic int    size_y = -1;        /* -y: Y extent */
  913. Xstatic int    size_z =  1;        /* -z: Z extent (default 1) */
  914. Xstatic int    size_i = -1;        /* -i: input image X extent */
  915. Xstatic int    size_j = -1;        /* -j: input image Y extent */
  916. Xstatic int    size_k = -1;        /* -k: input image Z extent */
  917. Xstatic LIC_Filter
  918. X              filter = LIC_Box;        /* -f: filter */
  919. Xstatic int    normalization = LIC_VARIABLE;
  920. X                     /* -n: normalization */
  921. Xstatic int    normalized    = FALSE;    /* -N: input vectors prenormalized */
  922. Xstatic double filter_length = 10.0;    /* -l: filter length */
  923. Xstatic double filter_frequency = 3.0;    /* -d: filter frequency */
  924. Xstatic double filter_phase  = 0.0;    /* -p: filter phase */
  925. Xstatic int    vary_length   = FALSE;    /* -L: vary length by vector mag. */
  926. Xstatic int    vary_speed    = FALSE;    /* -S: vary phase by vector mag. */
  927. Xstatic int    default_red   = -1;    /* -r: default red ... */
  928. Xstatic int    default_green = -1;    /* -g: default green ... */
  929. Xstatic int    default_blue  = -1;    /* -b: default blue ... */
  930. Xstatic int    default_alpha = -1;    /* -a: default alpha for zero vector */
  931. Xstatic int    verbose       = FALSE;    /* -v: display progress information */
  932. Xstatic int    version       = FALSE;    /* -V: display version information */
  933. X
  934. Xstatic char  *in_img;            /* input image filename */
  935. Xstatic char  *in_vec;            /* input vector field filename */
  936. Xstatic char  *out_img;            /* output image filename */
  937. X
  938. X
  939. X/*
  940. X *    Local variables
  941. X *    ===============
  942. X */
  943. Xstatic time_t     t0;            /* time we started/ended the LIC */
  944. Xstatic clock_t    t0_clk, tN_clk;    /*   calculation */
  945. Xstatic struct tms t0_tms, tN_tms;
  946. X
  947. Xstatic void      *in_img_addr,        /* virtual addresses of memory */
  948. X                 *in_vec_addr,        /*   mapped input image, input */
  949. X                 *out_img_addr;        /*   vector field and output image */
  950. Xstatic off_t      in_img_size,        /* sizes of image and vector files */
  951. X                  in_vec_size,
  952. X                  out_img_size;
  953. X#if !defined(HAS_MMAP)
  954. X    static int    out_img_fd;        /* output image file descriptor */
  955. X#endif
  956. X
  957. X
  958. X/*
  959. X *    Local routines
  960. X *    ==============
  961. X */
  962. Xint         main(int argc, char *argv[]);
  963. Xstatic void ParseArguments(int argc, char *argv[]);
  964. Xstatic void LoadFiles(void);
  965. Xstatic void UnloadFiles(void);
  966. Xstatic void LicFiles(void);
  967. Xstatic void PrintLicStatus(double PercentDone);
  968. Xstatic void PrintLicError(const char *message);
  969. X
  970. X
  971. Xint
  972. Xmain(int argc, char *argv[])
  973. X{
  974. X    ParseArguments(argc, argv);
  975. X    LoadFiles();
  976. X    LicFiles();
  977. X    UnloadFiles();
  978. X    exit(EXIT_SUCCESS);
  979. X    /*NOTREACHED*/
  980. X}
  981. X
  982. X
  983. Xstatic void
  984. XParseArguments(int argc, char *argv[])
  985. X{
  986. X    int          ch;
  987. X#if defined(HAS_HP_GETOPT)
  988. X    extern int   getopt(int, char * const [], const char *);
  989. X#else
  990. X    extern int   getopt(int, char **, char *);
  991. X#endif
  992. X    extern char *optarg;
  993. X    extern int   optind;
  994. X
  995. X    myname = strrchr(argv[0], '/');
  996. X    if (myname != NULL)
  997. X    myname++;
  998. X    else
  999. X    myname = argv[0];
  1000. X    while ((ch = getopt(argc, argv, "x:y:z:i:j:k:f:n:Nl:d:p:LSr:g:b:a:vV")) != EOF)
  1001. X    switch ((char)ch)
  1002. X    {
  1003. X        default:
  1004. X        case '?':
  1005. X        (void)fprintf(stderr, "%s: unknown option -%c\n",
  1006. X                  myname, ch);
  1007. X        (void)fprintf(stderr, usage, myname);
  1008. X        exit(EXIT_FAILURE);
  1009. X        /*NOTREACHED*/
  1010. X        case 'x':
  1011. X        size_x = atoi(optarg);
  1012. X        break;
  1013. X        case 'y':
  1014. X        size_y = atoi(optarg);
  1015. X        break;
  1016. X        case 'z':
  1017. X        size_z = atoi(optarg);
  1018. X        break;
  1019. X        case 'i':
  1020. X        size_i = atoi(optarg);
  1021. X        break;
  1022. X        case 'j':
  1023. X        size_j = atoi(optarg);
  1024. X        break;
  1025. X        case 'k':
  1026. X        size_k = atoi(optarg);
  1027. X        break;
  1028. X        case 'f':
  1029. X        if (strcasecmp(optarg, "box") == 0)
  1030. X            filter = LIC_Box;
  1031. X        else if (strcasecmp(optarg, "ripple") == 0)
  1032. X            filter = LIC_Ripple;
  1033. X        else if (strcasecmp(optarg, "ramp") == 0)
  1034. X            filter = LIC_Ramp;
  1035. X        else if (strcasecmp(optarg, "select") == 0)
  1036. X            filter = LIC_Select;
  1037. X        else
  1038. X        {
  1039. X            (void)fprintf(stderr, "%: invalid filter type %s:"
  1040. X                  " use 'box', 'ripple' or 'ramp'\n",
  1041. X                  myname, optarg);
  1042. X            exit(EXIT_FAILURE);
  1043. X            /*NOTREACHED*/
  1044. X        }
  1045. X        break;
  1046. X        case 'n':
  1047. X        if (strcasecmp(optarg, "fixed") == 0)
  1048. X            normalization = LIC_FIXED;
  1049. X        else if (strcasecmp(optarg, "variable") == 0)
  1050. X            normalization = LIC_VARIABLE;
  1051. X        else
  1052. X        {
  1053. X            (void)fprintf(stderr, "%: invalid normalization type %s:"
  1054. X                  " use 'fixed' or 'variable'\n",
  1055. X                  myname, optarg);
  1056. X            exit(EXIT_FAILURE);
  1057. X            /*NOTREACHED*/
  1058. X        }
  1059. X        break;
  1060. X        case 'N':
  1061. X        normalized = TRUE;
  1062. X        break;
  1063. X        case 'l':
  1064. X        filter_length = atof(optarg);
  1065. X        break;
  1066. X        case 'd':
  1067. X        filter_frequency = atof(optarg);
  1068. X        break;
  1069. X        case 'p':
  1070. X        filter_phase = atof(optarg);
  1071. X        break;
  1072. X        case 'L':
  1073. X        vary_length = TRUE;
  1074. X        break;
  1075. X        case 'S':
  1076. X        vary_speed = TRUE;
  1077. X        break;
  1078. X        case 'r':
  1079. X        default_red = (*optarg == '-')
  1080. X            ? -atoi(optarg+1)
  1081. X            :  atoi(optarg);
  1082. X        break;
  1083. X        case 'g':
  1084. X        default_green = (*optarg == '-')
  1085. X            ? -atoi(optarg+1)
  1086. X            :  atoi(optarg);
  1087. X        break;
  1088. X        case 'b':
  1089. X        default_blue = (*optarg == '-')
  1090. X            ? -atoi(optarg+1)
  1091. X            :  atoi(optarg);
  1092. X        break;
  1093. X        case 'a':
  1094. X        default_alpha = (*optarg == '-')
  1095. X            ? -atoi(optarg+1)
  1096. X            :  atoi(optarg);
  1097. X        break;
  1098. X        case 'v':
  1099. X        verbose = TRUE;
  1100. X        break;
  1101. X        case 'V':
  1102. X        version = TRUE;
  1103. X        break;
  1104. X    }
  1105. X
  1106. X    if (version)
  1107. X    (void)printf("%s: compiled with LIC library using %d byte %s pixels\n",
  1108. X             myname, LIC_ConfiguredPixelSize(),
  1109. X             LIC_ConfiguredPixelType());
  1110. X
  1111. X    if (argc != optind + 3)
  1112. X    {
  1113. X    (void)fprintf(stderr, "%s: need three file options\n", myname);
  1114. X    (void)fprintf(stderr, usage, myname);
  1115. X    exit(EXIT_FAILURE);
  1116. X    /*NOTREACHED*/
  1117. X    }
  1118. X    in_img = argv[optind + 0];
  1119. X    in_vec = argv[optind + 1];
  1120. X    out_img = argv[optind + 2];
  1121. X    if (size_x == -1 || size_y == -1)
  1122. X    {
  1123. X    (void)fprintf(stderr, "%s: must specify both -x and -y options ...\n",
  1124. X              myname);
  1125. X    exit(EXIT_FAILURE);
  1126. X    /*NOTREACHED*/
  1127. X    }
  1128. X    if (size_i == -1)
  1129. X    size_i = size_x;
  1130. X    if (size_j == -1)
  1131. X    size_j = size_y;
  1132. X    if (size_k == -1)
  1133. X    size_k = size_z;
  1134. X}
  1135. X
  1136. X
  1137. X#if defined(HAS_MMAP)
  1138. X
  1139. X/*
  1140. X * Versions of LoadFiles and UnloadFiles for systems that support mmap(2).
  1141. X */
  1142. X
  1143. Xstatic void
  1144. XLoadFiles(void)
  1145. X    /*
  1146. X     * Map input and output files into virtual memory.  Make sure that the
  1147. X     * input files match the user specified sizes and pre-extend the
  1148. X     * output file to its final length.
  1149. X     */
  1150. X{
  1151. X    char        c;
  1152. X    int         fd;
  1153. X    struct stat stbuf;
  1154. X
  1155. X    /*
  1156. X     * Map input vector field into virtual memory.  (We always map the
  1157. X     * input vector field first because it will always be larger than
  1158. X     * either of the other two files.  This should, hopefully, give the
  1159. X     * operating system the best possible chance to find virtual address
  1160. X     * space for the three files.)
  1161. X     */
  1162. X    fd = open(in_vec, O_RDONLY);
  1163. X    if (fd < 0)
  1164. X    {
  1165. X    (void)fprintf(stderr, "%s: unable to open %s: %s\n",
  1166. X              myname, in_vec, strerror(errno));
  1167. X    exit(EXIT_FAILURE);
  1168. X    /*NOTREACHED*/
  1169. X    }
  1170. X    if (fstat(fd, &stbuf) < 0)
  1171. X    {
  1172. X    (void)fprintf(stderr, "%s: unable to stat %s: %s\n",
  1173. X              myname, in_vec, strerror(errno));
  1174. X    exit(EXIT_FAILURE);
  1175. X    /*NOTREACHED*/
  1176. X    }
  1177. X    in_vec_size = stbuf.st_size;
  1178. X    if (size_x * size_y * size_z * sizeof(float) * (size_z == 1 ? 2 : 3)
  1179. X    != in_vec_size)
  1180. X    {
  1181. X    (void)fprintf(stderr, "%s: input vector field size doesn't"
  1182. X              " match %d * x=%d * y=%d * z=%d\n",
  1183. X              myname, sizeof(float) * (size_z == 1 ? 2 : 3),
  1184. X              size_x, size_y, size_z);
  1185. X    exit(EXIT_FAILURE);
  1186. X    /*NOTREACHED*/
  1187. X    }
  1188. X    in_vec_addr = mmap((void *)0, (int)in_vec_size, PROT_READ,
  1189. X               MAP_SHARED, fd, (off_t)0);
  1190. X    if ((int)in_vec_addr == -1)
  1191. X    {
  1192. X    (void)fprintf(stderr, "%s: unable to map %s: %s\n",
  1193. X              myname, in_vec, strerror(errno));
  1194. X    exit(EXIT_FAILURE);
  1195. X    /*NOTREACHED*/
  1196. X    }
  1197. X    (void)close(fd);
  1198. X
  1199. X    /*
  1200. X     * Map input image into virtual memory.
  1201. X     */
  1202. X    fd = open(in_img, O_RDONLY);
  1203. X    if (fd < 0)
  1204. X    {
  1205. X    (void)fprintf(stderr, "%s: unable to open %s: %s\n",
  1206. X              myname, in_img, strerror(errno));
  1207. X    exit(EXIT_FAILURE);
  1208. X    /*NOTREACHED*/
  1209. X    }
  1210. X    if (fstat(fd, &stbuf) < 0)
  1211. X    {
  1212. X    (void)fprintf(stderr, "%s: unable to stat %s: %s\n",
  1213. X              myname, in_img, strerror(errno));
  1214. X    exit(EXIT_FAILURE);
  1215. X    /*NOTREACHED*/
  1216. X    }
  1217. X    in_img_size = stbuf.st_size;
  1218. X    if (size_i * size_j * size_k * LIC_ConfiguredPixelSize() != in_img_size)
  1219. X    {
  1220. X    (void)fprintf(stderr, "%s: input image size doesn't"
  1221. X              " match %d * x=%d * y=%d * z=%d\n",
  1222. X              myname, LIC_ConfiguredPixelSize(),
  1223. X              size_i, size_j, size_k);
  1224. X    exit(EXIT_FAILURE);
  1225. X    /*NOTREACHED*/
  1226. X    }
  1227. X    in_img_addr = mmap((void *)0, (int)in_img_size, PROT_READ,
  1228. X               MAP_SHARED, fd, (off_t)0);
  1229. X    if ((int)in_img_addr == -1)
  1230. X    {
  1231. X    (void)fprintf(stderr, "%s: unable to map %s: %s\n",
  1232. X              myname, in_img, strerror(errno));
  1233. X    exit(EXIT_FAILURE);
  1234. X    /*NOTREACHED*/
  1235. X    }
  1236. X    (void)close(fd);
  1237. X
  1238. X    /*
  1239. X     * Map output image into virtual memory.
  1240. X     */
  1241. X    fd = open(out_img, O_CREAT|O_RDWR|O_EXCL, 0666);
  1242. X    if (fd < 0)
  1243. X    {
  1244. X    (void)fprintf(stderr, "%s: unable to open %s: %s\n",
  1245. X              myname, out_img, strerror(errno));
  1246. X    exit(EXIT_FAILURE);
  1247. X    /*NOTREACHED*/
  1248. X    }
  1249. X    out_img_size = size_x * size_y * size_z * LIC_ConfiguredPixelSize();
  1250. X    (void)lseek(fd, out_img_size - 1, SEEK_SET);
  1251. X    c = '\0';
  1252. X    if (write(fd, &c, 1) != 1)
  1253. X    {
  1254. X    (void)fprintf(stderr, "%s: unable to preextend %s to %ld bytes: %s\n",
  1255. X              myname, out_img, strerror(errno));
  1256. X    exit(EXIT_FAILURE);
  1257. X    /*NOTREACHED*/
  1258. X    }
  1259. X    out_img_addr = mmap((void *)0, (int)out_img_size, PROT_WRITE,
  1260. X            MAP_SHARED, fd, (off_t)0);
  1261. X    if ((int)out_img_addr == -1)
  1262. X    {
  1263. X    (void)fprintf(stderr, "%s: unable to map %s: %s\n",
  1264. X              myname, out_img, strerror(errno));
  1265. X    exit(EXIT_FAILURE);
  1266. X    /*NOTREACHED*/
  1267. X    }
  1268. X    (void)close(fd);
  1269. X}
  1270. X
  1271. X
  1272. Xstatic void
  1273. XUnloadFiles(void)
  1274. X    /*
  1275. X     * Unmap input and output files from virtual memory.
  1276. X     */
  1277. X{
  1278. X    (void)munmap(in_img_addr, (int)in_img_size);
  1279. X    (void)munmap(in_vec_addr, (int)in_vec_size);
  1280. X    (void)munmap(out_img_addr, (int)out_img_size);
  1281. X}
  1282. X
  1283. X#else /* HAS_MMAP */
  1284. X
  1285. X/*
  1286. X * Versions of LoadFiles and UnloadFiles for systems that *don't* support
  1287. X * mmap(2).
  1288. X */
  1289. X
  1290. Xstatic void
  1291. XLoadFiles(void)
  1292. X    /*
  1293. X     * Load input files into malloc'ed memory and malloc memory for
  1294. X     * output image.  Make sure that the input files match the user
  1295. X     * specified sizes.  Preallocate output image on disk in order to
  1296. X     * make sure we have the space to output the image later in
  1297. X     * UnloadFiles.
  1298. X     */
  1299. X{
  1300. X    int         fd;
  1301. X    struct stat stbuf;
  1302. X
  1303. X    /*
  1304. X     * Load input image into malloc'ed memory.
  1305. X     */
  1306. X    fd = open(in_img, O_RDONLY);
  1307. X    if (fd < 0)
  1308. X    {
  1309. X    (void)fprintf(stderr, "%s: unable to open %s: %s\n",
  1310. X              myname, in_img, strerror(errno));
  1311. X    exit(EXIT_FAILURE);
  1312. X    /*NOTREACHED*/
  1313. X    }
  1314. X    if (fstat(fd, &stbuf) < 0)
  1315. X    {
  1316. X    (void)fprintf(stderr, "%s: unable to stat %s: %s\n",
  1317. X              myname, in_img, strerror(errno));
  1318. X    exit(EXIT_FAILURE);
  1319. X    /*NOTREACHED*/
  1320. X    }
  1321. X    in_img_size = stbuf.st_size;
  1322. X    if (size_i * size_j * size_k * LIC_ConfiguredPixelSize()
  1323. X    != in_img_size)
  1324. X    {
  1325. X    (void)fprintf(stderr, "%s: input image size doesn't"
  1326. X              " match %d * x=%d * y=%d * z=%d\n",
  1327. X              myname, LIC_ConfiguredPixelSize(),
  1328. X              size_i, size_j, size_k);
  1329. X    exit(EXIT_FAILURE);
  1330. X    /*NOTREACHED*/
  1331. X    }
  1332. X    in_img_addr = malloc(in_img_size);
  1333. X    if (in_img_addr == NULL)
  1334. X    {
  1335. X    (void)fprintf(stderr, "%s: unable to allocate %u bytes for"
  1336. X              " input image\n", myname, in_img_size);
  1337. X    exit(EXIT_FAILURE);
  1338. X    /*NOTREACHED*/
  1339. X    }
  1340. X    if (read(fd, in_img_addr, in_img_size) != in_img_size)
  1341. X    {
  1342. X    (void)fprintf(stderr, "%s: unable to read %u bytes from %s\n",
  1343. X              myname, in_img_size, in_img);
  1344. X    exit(EXIT_FAILURE);
  1345. X    /*NOTREACHED*/
  1346. X    }
  1347. X    (void)close(fd);
  1348. X
  1349. X    /*
  1350. X     * Load input vector field into malloc'ed memory.
  1351. X     */
  1352. X    fd = open(in_vec, O_RDONLY);
  1353. X    if (fd < 0)
  1354. X    {
  1355. X    (void)fprintf(stderr, "%s: unable to open %s: %s\n",
  1356. X              myname, in_vec, strerror(errno));
  1357. X    exit(EXIT_FAILURE);
  1358. X    /*NOTREACHED*/
  1359. X    }
  1360. X    if (fstat(fd, &stbuf) < 0)
  1361. X    {
  1362. X    (void)fprintf(stderr, "%s: unable to stat %s: %s\n",
  1363. X              myname, in_vec, strerror(errno));
  1364. X    exit(EXIT_FAILURE);
  1365. X    /*NOTREACHED*/
  1366. X    }
  1367. X    in_vec_size = stbuf.st_size;
  1368. X    if (size_x * size_y * size_z * sizeof(float) * (size_z == 1 ? 2 : 3)
  1369. X    != in_vec_size)
  1370. X    {
  1371. X    (void)fprintf(stderr, "%s: input vector field size doesn't"
  1372. X              " match %d * x=%d * y=%d * z=%d\n",
  1373. X              myname, (size_z == 1 ? 2 : 3),
  1374. X              size_x, size_y, size_z);
  1375. X    exit(EXIT_FAILURE);
  1376. X    /*NOTREACHED*/
  1377. X    }
  1378. X    in_vec_addr = malloc(in_vec_size);
  1379. X    if (in_vec_addr == NULL)
  1380. X    {
  1381. X    (void)fprintf(stderr, "%s: unable to allocate %u bytes for"
  1382. X              " input vector field\n", myname, in_vec_size);
  1383. X    exit(EXIT_FAILURE);
  1384. X    /*NOTREACHED*/
  1385. X    }
  1386. X    if (read(fd, in_vec_addr, in_vec_size) != in_vec_size)
  1387. X    {
  1388. X    (void)fprintf(stderr, "%s: unable to read %u bytes from %s\n",
  1389. X              myname, in_vec_size, in_vec);
  1390. X    exit(EXIT_FAILURE);
  1391. X    /*NOTREACHED*/
  1392. X    }
  1393. X    (void)close(fd);
  1394. X
  1395. X    /*
  1396. X     * Open output file and allocate space for output image (both in
  1397. X     * memory and on disk).
  1398. X     */
  1399. X    out_img_fd = open(out_img, O_CREAT|O_WRONLY|O_EXCL, 0666);
  1400. X    if (out_img_fd < 0)
  1401. X    {
  1402. X    (void)fprintf(stderr, "%s: unable to open %s: %s\n",
  1403. X              myname, out_img, strerror(errno));
  1404. X    exit(EXIT_FAILURE);
  1405. X    /*NOTREACHED*/
  1406. X    }
  1407. X    out_img_size = size_x * size_y * size_z * LIC_ConfiguredPixelSize();
  1408. X    out_img_addr = malloc(out_img_size);
  1409. X    if (out_img_addr == NULL)
  1410. X    {
  1411. X    (void)fprintf(stderr, "%s: unable to allocate %u bytes for"
  1412. X              " output image\n", myname, out_img_size);
  1413. X    exit(EXIT_FAILURE);
  1414. X    /*NOTREACHED*/
  1415. X    }
  1416. X    if (write(out_img_fd, out_img_addr, out_img_size) != out_img_size)
  1417. X    {
  1418. X    (void)fprintf(stderr, "%s: unable to extend %s to %u bytes\n",
  1419. X              myname, out_img, out_img_size);
  1420. X    exit(EXIT_FAILURE);
  1421. X    /*NOTREACHED*/
  1422. X    }
  1423. X    (void)lseek(out_img_fd, (off_t)0, SEEK_SET);
  1424. X}
  1425. X
  1426. X
  1427. Xstatic void
  1428. XUnloadFiles(void)
  1429. X    /*
  1430. X     * Free malloc'ed space used by input files, write output image and
  1431. X     * free up it's malloc'ed memory.
  1432. X     */
  1433. X{
  1434. X    if (write(out_img_fd, out_img_addr, out_img_size) != out_img_size)
  1435. X    (void)fprintf(stderr, "%s: unable to write %u bytes to %s!\n",
  1436. X              myname, out_img_size, out_img);
  1437. X    close(out_img_fd);
  1438. X    free(in_img_addr);
  1439. X    free(in_vec_addr);
  1440. X    free(out_img_addr);
  1441. X}
  1442. X
  1443. X#endif /* HAS_MMAP */
  1444. X
  1445. X
  1446. Xstatic void
  1447. XLicFiles()
  1448. X    /*
  1449. X     * Execute LIC algorithm on arguments.
  1450. X     */
  1451. X{
  1452. X    LIC *lic;
  1453. X
  1454. X    lic = LIC_Create((unsigned char *)in_img_addr, size_i, size_j, size_k,
  1455. X             (float *)        in_vec_addr, size_x, size_y, size_z,
  1456. X             (unsigned char *)out_img_addr,
  1457. X             filter, normalization, normalized,
  1458. X             filter_length,
  1459. X             filter_frequency,
  1460. X             vary_length,
  1461. X             vary_speed && (filter == LIC_Ripple),
  1462. X             default_red, default_green, default_blue, default_alpha,
  1463. X             verbose ? PrintLicStatus : (void (*)(double))NULL,
  1464. X             PrintLicError);
  1465. X    if (lic == (LIC *)NULL)
  1466. X    {
  1467. X    (void)fprintf(stderr, "%s: LIC_Create returned NULL!\n", myname);
  1468. X    return;
  1469. X    }
  1470. X    LIC_ChangePhase(lic, filter_phase);
  1471. X
  1472. X    /*
  1473. X     * Build integral tables here so the build isn't computed as part of
  1474. X     * time to compute the image.
  1475. X     */
  1476. X    LIC_BuildIntegralTables(lic);
  1477. X
  1478. X    t0 = time(NULL);
  1479. X    t0_clk = times(&t0_tms);
  1480. X    LIC_ComputeImage(lic);
  1481. X    tN_clk = times(&tN_tms);
  1482. X
  1483. X    LIC_Destroy(lic);
  1484. X
  1485. X    if (verbose)
  1486. X    {
  1487. X    long   clk_tck = sysconf(_SC_CLK_TCK);
  1488. X#if defined(HAS_OLD_TIMES)
  1489. X        double wall   = (double)(time(NULL)       - t0);
  1490. X#else
  1491. X    double wall   = (double)(tN_clk           - t0_clk          )/clk_tck;
  1492. X#endif
  1493. X    double cpu    = (double)(tN_tms.tms_utime - t0_tms.tms_utime)/clk_tck;
  1494. X    double ncells = (double)(size_x * size_y * size_z);
  1495. X
  1496. X    if (wall == 0.0)
  1497. X        wall = 1.0e-6;
  1498. X    if (cpu == 0.0)
  1499. X        cpu  = 1.0e-6;
  1500. X    (void)printf("\n"
  1501. X             "Performance statistics\n"
  1502. X             "--------------------------------\n"
  1503. X             "CPU utilization             = %d%%\n"
  1504. X             "CPU time (seconds)          = %.2f\n"
  1505. X             "Cells processed per second  = %.2f\n"
  1506. X             "Ave loop count              = %.2f\n"
  1507. X             "Ave length                  = %.2f\n",
  1508. X             (int)(cpu/wall*100),
  1509. X             cpu,
  1510. X             ncells/cpu,
  1511. X             lic->TotalLoopCount/ncells/2,
  1512. X             lic->TotalLength/ncells/2);
  1513. X    }
  1514. X}
  1515. X
  1516. X
  1517. Xstatic void
  1518. XPrintLicStatus(double PercentDone)
  1519. X    /*
  1520. X     * Output the current LIC progress statistics.  (Only called if verbose
  1521. X     * is TRUE.)
  1522. X     */
  1523. X{
  1524. X    if (PercentDone == 0.0)
  1525. X    (void)printf("%s: %6.2f%% done ...\r", myname, PercentDone);
  1526. X    else
  1527. X    {
  1528. X    /*
  1529. X     * Should really use a smoothed exponentially decaying sample series
  1530. X     * of CPU utilization to predict future CPU utilization to get better
  1531. X     * estimates of completion times.  But, I mean, get real.  This is
  1532. X     * really just a convenience! :-)  So, we cheat and just use CPU
  1533. X     * utilization since the start of the computation to predict future
  1534. X     * CPU ultilization.
  1535. X     */
  1536. X    time_t t  = time(NULL);
  1537. X    time_t tN = t0 + (unsigned int)((double)(t-t0) * 100/PercentDone);
  1538. X    char      *cp = ctime(&tN);
  1539. X    char       cbuf[26];
  1540. X
  1541. X    (void)strcpy(cbuf, cp);
  1542. X    cbuf[24] = '\0';
  1543. X    (void)printf("%s: %6.2f%% done, estimated completion: %s\r",
  1544. X             myname, PercentDone, cbuf);
  1545. X    }
  1546. X    fflush(stdout);
  1547. X}
  1548. X
  1549. X
  1550. Xstatic void
  1551. XPrintLicError(const char *message)
  1552. X    /*
  1553. X     * Print an error reported by the LIC subroutine library.
  1554. X     */
  1555. X{
  1556. X    (void)fprintf(stderr, "%s: LIC library error: %s\n", myname);
  1557. X}
  1558. END_OF_FILE
  1559.   if test 22645 -ne `wc -c <'lic.1.2/lic/lic.c'`; then
  1560.     echo shar: \"'lic.1.2/lic/lic.c'\" unpacked with wrong size!
  1561.   fi
  1562.   # end of 'lic.1.2/lic/lic.c'
  1563. fi
  1564. echo shar: End of archive 6 \(of 10\).
  1565. cp /dev/null ark6isdone
  1566. MISSING=""
  1567. for I in 1 2 3 4 5 6 7 8 9 10 ; do
  1568.     if test ! -f ark${I}isdone ; then
  1569.     MISSING="${MISSING} ${I}"
  1570.     fi
  1571. done
  1572. if test "${MISSING}" = "" ; then
  1573.     echo You have unpacked all 10 archives.
  1574.     rm -f ark[1-9]isdone ark[1-9][0-9]isdone
  1575. else
  1576.     echo You still must unpack the following archives:
  1577.     echo "        " ${MISSING}
  1578. fi
  1579. exit 0
  1580. exit 0 # Just in case...
  1581.