diff --git a/lib/bpdecode144.f90 b/lib/bpdecode144.f90
index b67e5f944..bcbd6c16c 100644
--- a/lib/bpdecode144.f90
+++ b/lib/bpdecode144.f90
@@ -1,5 +1,6 @@
 subroutine pltanh(x,y)
   isign=+1
+  z=x
   if( x.lt.0 ) then
     isign=-1
     z=abs(x)
@@ -24,6 +25,7 @@ end subroutine pltanh
 
 subroutine platanh(x,y)
   isign=+1
+  z=x
   if( x.lt.0 ) then
     isign=-1
     z=abs(x)
diff --git a/lib/bpdecode40.f90 b/lib/bpdecode40.f90
index 7839980b3..3a35bab49 100644
--- a/lib/bpdecode40.f90
+++ b/lib/bpdecode40.f90
@@ -136,9 +136,9 @@ do iter=0,maxiterations
   do j=1,N
     do i=1,ncw
       ichk=Mn(i,j)  ! Mn(:,j) are the checks that include bit j
-      Tmn=product(tanhtoc(1:nrw(i),ichk),mask=Nm(1:nrw(i),ichk).ne.j)
+      Tmn=product(tanhtoc(1:nrw(ichk),ichk),mask=Nm(1:nrw(ichk),ichk).ne.j)
       call platanh(-Tmn,y)
-      tov(i,j)=y
+      tov(i,j)=2*y
     enddo
   enddo