From 8970a140e08276bd1c8677d7a1919e64b113ad79 Mon Sep 17 00:00:00 2001
From: Joe Taylor <k1jt@arrl.org>
Date: Tue, 2 Oct 2012 18:08:44 +0000
Subject: [PATCH] Working toward decode capability.

git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/wsjtx@2628 ab8295b8-cf94-4d9e-aec4-7959e3be5d79
---
 libm65/genjt9.f90 |   5 +-
 libm65/jt9.f90    |   4 +-
 libm65/jt9sim.f90 |   9 ++-
 libm65/spec9.f90  | 142 ++++++++++++++++++++++++++++++++++++++++++++++
 libm65/sync9.f90  |  23 +++-----
 mainwindow.cpp    |   2 +-
 6 files changed, 164 insertions(+), 21 deletions(-)
 create mode 100644 libm65/spec9.f90

diff --git a/libm65/genjt9.f90 b/libm65/genjt9.f90
index 8af75ba18..b401ee479 100644
--- a/libm65/genjt9.f90
+++ b/libm65/genjt9.f90
@@ -31,7 +31,10 @@ subroutine genjt9(message,minutes,msgsent,d6)
   call encode232(d1,nsym2,d2)        !Convolutional code, K=32, r=1/2
   call interleave9(d2,1,d3)          !Interleave the single bits
   call packbits(d3,nsym2,3,d4)       !Pack 3-bit groups into words
-  call graycode(d4,69,1,d5)          !Apply Gray code
+  d5=d4
+  print*,d5
+  
+!  call graycode(d4,69,1,d5)          !Apply Gray code
 
 ! Insert sync symbols (ntone=0) and add 1 to the data-tone numbers.
   j=0
diff --git a/libm65/jt9.f90 b/libm65/jt9.f90
index 2ab66b5d0..4924a39ba 100644
--- a/libm65/jt9.f90
+++ b/libm65/jt9.f90
@@ -108,10 +108,10 @@ program jt9
 3001    format(i8,2f12.3)
      enddo
 
-     print*,npts8,npts8/1500.0
      nutc=nutc0
      nstandalone=1
-     call sync9(ss,tstep,f0a,df3)
+     call sync9(ss,tstep,f0a,df3,lagpk,fpk)
+     call spec9(c0,npts8,nsps,f0a,lagpk,fpk)
 !     call decode0(dd,ss,savg,nstandalone,nfsample)
   enddo
 
diff --git a/libm65/jt9sim.f90 b/libm65/jt9sim.f90
index 86466763d..870ed436e 100644
--- a/libm65/jt9sim.f90
+++ b/libm65/jt9sim.f90
@@ -155,9 +155,11 @@ program jt9sim
         t5(j)=itone(i)-1
      enddo
 
-     call graycode(t5,69,-1,t4)
+!     call graycode(t5,69,-1,t4)
+     t4=t5
      call unpackbits(t4,69,3,t3)
-     call interleave9(t3,-1,t2)
+!     call interleave9(t3,-1,t2)
+     t2=t3
 
      nbits=72
      ndelta=17
@@ -168,9 +170,12 @@ program jt9sim
         i4=i4+128
         t2(i)=i1
      enddo
+
+     print*,t2
   
      call fano232(t2,nbits+31,mettab,ndelta,limit,t1,ncycles,metric,ierr,   &
           maxmetric,maxnp)
+     print*,ncycles
 
      nbytes=(nbits+7)/8
      do i=1,nbytes
diff --git a/libm65/spec9.f90 b/libm65/spec9.f90
new file mode 100644
index 000000000..0c0c8c70a
--- /dev/null
+++ b/libm65/spec9.f90
@@ -0,0 +1,142 @@
+subroutine spec9(c0,npts8,nsps,f0a,lagpk,fpk)
+
+  parameter (MAXFFT=31500)
+  complex c0(0:npts8-1)
+  real s(0:MAXFFT-1)
+  real ssym(0:9,184)
+  complex c(0:MAXFFT-1)
+  integer*1 softbit(207)
+  integer ibit(207)
+
+  integer*1 t1(13)              !72 bits and zero tail as 8-bit bytes
+  integer*4 t4(69)              !Symbols from t5, values 0-7
+  integer*4 mettab(0:255,0:1)
+  integer*1 tmp(72)
+  character*22 msg
+
+  integer isync(85)
+  integer ii(16)                       !Locations of sync symbols
+  data ii/1,6,11,16,21,26,31,39,45,51,57,63,69,75,81,85/
+
+  isync=0
+  do i=1,16
+     isync(ii(i))=1
+  enddo
+
+  idt=-400
+  idf=0.
+  fshift=fpk-f0a + 0.1*idf
+  twopi=8.0*atan(1.0)
+  dphi=twopi*fshift/1500.0
+  nsps8=nsps/8
+  nfft=nsps8
+  df=1500.0/nfft
+  s=0.
+  istart=lagpk*nsps8 + idt
+  nsym=min((npts8-istart)/nsps8,85)
+  print*,f0a,lagpk,fpk,fshift,nsym
+
+  do j=0,nsym-1
+     ia=j*nsps8 + istart
+     ib=ia+nsps8-1
+     c(0:nfft-1)=c0(ia:ib)
+
+     phi=0.
+     do i=0,nfft-1
+        phi=phi + dphi
+        c(i)=c(i) * cmplx(cos(phi),-sin(phi))
+     enddo
+
+     call four2a(c,nfft,1,-1,1)
+     do i=0,nfft-1
+        sx=real(c(i))**2 + aimag(c(i))**2
+        if(i.le.8) ssym(i,1+j)=sx
+        s(i)=s(i) + sx
+     enddo
+  enddo
+
+!  sync=0.
+!  do i=1,16
+!     sync=sync + ssym(0,ii(i))
+!  enddo
+
+  do i=0,nfft-1
+     freq = f0a + i*df
+     write(71,3001) i,i*df,freq,s(i)
+3001 format(i8,3f12.3)
+  enddo
+
+  do j=1,nsym
+     write(72,3002) j,(ssym(i,j),i=0,8)
+3002 format(i3,9f7.2)
+  enddo
+
+  m0=3
+  ntones=8
+  k=0
+  do j=1,nsym
+!     print*,j,k,nsym
+     if(isync(j).eq.1) cycle
+     do m=m0-1,0,-1                   !Get bit-wise soft symbols
+        n=2**m
+        r1=0.
+        r2=0.
+        do ig=0,ntones-1
+!           i=igray(ig,-1)
+           i=ig
+           if(iand(i,n).ne.0) then
+              r1=max(r1,ssym(i+1,j))
+           else
+              r2=max(r2,ssym(i+1,j))
+           endif
+        enddo
+        k=k+1
+        softbit(k)=min(127,max(-127,nint(10.0*(r1-r2)))) + 128
+     enddo
+  enddo
+
+  ibit=0
+  do i=1,207
+     if(softbit(i).lt.0) ibit(i)=1
+  enddo
+  print*,ibit
+
+! Get the metric table
+  bias=0.37                          !To be optimized, in decoder program
+  scale=10                           !  ... ditto ...
+  open(19,file='met8.21',status='old')
+
+  do i=0,255
+     read(19,*) x00,x0,x1
+     mettab(i,0)=nint(scale*(x0-bias))
+     mettab(i,1)=nint(scale*(x1-bias))    !### Check range, etc.  ###
+  enddo
+  close(19)
+  nbits=72
+  ndelta=17
+  limit=1000
+
+  print*,softbit
+  call fano232(softbit,nbits+31,mettab,ndelta,limit,t1,ncycles,metric,ierr,   &
+          maxmetric,maxnp)
+  print*,ncycles
+
+  nbytes=(nbits+7)/8
+  do i=1,nbytes
+     n=t1(i)
+     t4(i)=iand(n,255)
+  enddo
+  call unpackbits(t4,nbytes,8,tmp)
+  call packbits(tmp,12,6,t4)
+  do i=1,12
+     if(t4(i).lt.128) t1(i)=t4(i)
+     if(t4(i).ge.128) t1(i)=t4(i)-256
+  enddo
+  do i=1,12
+     t4(i)=t1(i)
+  enddo
+  call unpackmsg(t4,msg)         !Unpack decoded msg
+  print*,msg
+
+  return
+end subroutine spec9
diff --git a/libm65/sync9.f90 b/libm65/sync9.f90
index 4c8a5c725..b989c1035 100644
--- a/libm65/sync9.f90
+++ b/libm65/sync9.f90
@@ -1,25 +1,18 @@
-subroutine sync9(ss,tstep,f0a,df3)
+subroutine sync9(ss,tstep,f0a,df3,lagpk,fpk)
 
   parameter (NSMAX=22000)            !Max length of saved spectra
   real ss(184,NSMAX)
 
-  integer ii(16)
-  integer isync(85)             !Sync vector
+  integer ii0(16)
+  integer ii(16)                     !Locations of sync half-symbols
+  data ii/1,11,21,31,41,51,61,77,89,101,113,125,137,149,161,169/
+  integer isync(85)                  !Sync vector for half-symbols
   data isync/                                    &
        1,0,0,0,0,1,0,0,0,0,1,0,0,0,0,1,0,0,0,0,  &
        1,0,0,0,0,1,0,0,0,0,1,0,0,0,0,0,0,0,1,0,  &
        0,0,0,0,1,0,0,0,0,0,1,0,0,0,0,0,1,0,0,0,  &
        0,0,1,0,0,0,0,0,1,0,0,0,0,0,1,0,0,0,0,0,  &
        1,0,0,0,1/
-  
-  ii=0
-  k=0
-  do i=1,85
-     if(isync(i).eq.1) then
-        k=k+1
-        ii(k)=2*i-1
-     endif
-  enddo
 
   nz=1000.0/df3
 
@@ -40,9 +33,9 @@ subroutine sync9(ss,tstep,f0a,df3)
      enddo
   enddo
 
-  freq=f0a + (npk-1)*df3
-  write(*,1010) lagpk,npk,freq
-1010 format('lagpk:',i4,'   npk:',i6,'   freq:',f8.2)
+  fpk=f0a + (npk-1)*df3
+  write(*,1010) lagpk,npk,fpk
+1010 format('lagpk:',i4,'   npk:',i6,'   fpk:',f8.2)
 
   do lag=-lagmax,lagmax
      sum=0.
diff --git a/mainwindow.cpp b/mainwindow.cpp
index a120db1b3..48d16902e 100644
--- a/mainwindow.cpp
+++ b/mainwindow.cpp
@@ -1,4 +1,4 @@
-//------------------------------------------------------------- MainWindow
+//-------------------------------------------------------------- MainWindow
 #include "mainwindow.h"
 #include "ui_mainwindow.h"
 #include "devsetup.h"