mycompsub_real.f

FFT Dual Recycling Code Results on Advanced LIGO and the Caltech 40m IFO

c**********************************************************************c
c********************** sum two real arrays ***************************c

        subroutine myreal_somme (c1,c2,c,nvecs)
	
        implicit  none
        integer*4 nvecs,i
        real*8    c1(nvecs),c2(nvecs),c(nvecs)

        do 100 i = 1, nvecs
           c(i) = c1(i) + c2(i)
 100    continue

        return
        end

c**********************************************************************c
c********************** subtract two real arrays **********************c

        subroutine myreal_diffr (c1,c2,c,nvecs)

        implicit  none
        integer*4 nvecs,i
        real*8    c1(nvecs),c2(nvecs),c(nvecs)

        do 100 i = 1, nvecs
           c(i) = c1(i) - c2(i)
 100    continue

        return
        end

c**********************************************************************c
c************* multiply an real array with a real constant ************c

        subroutine myreal_mult (in,const,out,nvecs)

        implicit  none
        integer*4 nvecs,i
        real*8    in(nvecs),const,out(nvecs)

        do 100 i = 1, nvecs
           out(i) = const*in(i)
 100    continue

        return
        end

c**********************************************************************c
c**************** scalar sum product of two real arrays ***************c

        subroutine myrescal (c1,c2,ps,nvecs)

        implicit  none
        integer*4 nvecs,i
        real*8    c1(nvecs),c2(nvecs),ps

        ps = 0.0

        do 100 i = 1, nvecs
           ps = ps + c1(i)*c2(i)
 100    continue

        return
        end

c**********************************************************************c
c**********************************************************************c