CTITLESAFXINB -- CORRELATION FOR SPATIAL INTERPOLATION 00010006 C*********************************************************************** C COPYRIGHT ATLANTIC RICHFIELD COMPANY 1991 * C*********************************************************************** CA AUTHOR D CORRIGAN 00020003 CA DESIGNER D CORRIGAN 00030003 CA LANGUAGE FORTRAN 77 00040003 CA SYSTEM IBM / CRAY 00050005 CA WRITTEN AUGUST 1991 00060003 C REVISED 02-10-92 JJC RENAMED SAFXIN2 TO SAFXINB. 00070003 CA 00080003 CA CALL SAFXINB( LX,R,C,A,B,OP,LOP ) 00090004 CA 00100004 CA IN/OUT ARGUMENT TYPE DESCRIPTION 00110004 CA 00120004 CA INPUT LX I4 NUMBER OF TRACES IN SPATIAL FILTER 00130004 CA INPUT R C4 AUTOCORRELLATION 00140004 CA INPUT OP R4 INTERPOLATION OPERATOR 00150004 CA INPUT LOP I4 LENGTH OF OP 00160004 CA OUTPUT C C4 CROSS CORRELATION 00170004 CA OUTPUT A R4 AUXILIARY ARRAYS 00180004 CA OUTPUT B R4 AUXILIARY ARRAYS 00190004 CA 00200003 CA CORRELATION FOR SPATIAL INTERPOLATION USING SPATIAL PREDICTION 00210006 CA METHOD. 00220006 CA 00230003 C 00240000 SUBROUTINE SAFXINB( LX,R,C,A,B,OP,LOP ) 00250002 C 00260000 IMPLICIT INTEGER (A-Z) 00270007 C 00280000 COMPLEX R(1),C(1) 00290000 C 00300007 REAL A(1),B(1) 00310000 REAL OP(1) 00320000 C 00330007 REAL AIMAG 00340007 REAL REAL 00350007 REAL CI 00360007 REAL CR 00370007 REAL SDOT 00380009 C 00390000 C INTERPOLATE INPUT AUTOCORRELATION FUNCTION 00400000 C 00410000 LH = LX/2 00420000 LOH = LOP/2 00430000 IR = LOH 00440000 DO 100 I = 1,LOH-1 00450000 A(I) = REAL(R(IR)) 00460000 B(I) = -AIMAG(R(IR)) 00470000 100 IR = IR - 1 00480000 C 00490000 IR = 1 00500000 DO 200 I = LOH,LOH+LX-1 00510000 A(I) = REAL(R(IR)) 00520000 B(I) = AIMAG(R(IR)) 00530000 200 IR = IR + 1 00540000 C 00550000 IC = LH + 1 00560000 JC = LH 00570000 DO 300 I = 1,LH 00580000 CR = SDOT( LOP,OP,1,A(I),1 ) 00590000 CI = SDOT( LOP,OP,1,B(I),1 ) 00600000 C(IC) = CMPLX( CR,CI ) 00610000 C(JC) = CMPLX( CR,-CI) 00620000 IC = IC + 1 00630000 300 JC = JC - 1 00640000 C 00650000 C 00660000 RETURN 00670000 END 00680000