diff --git a/LICENSE b/LICENSE
index f30830cd..89db0950 100644
--- a/LICENSE
+++ b/LICENSE
@@ -1,4 +1,4 @@
-              Parallel Sparse BLAS  version 3.2
+              Parallel Sparse BLAS  version 3.3
     (C) Copyright 2006, 2007, 2008, 2009, 2010, 2012, 2013, 2014
 
                        Salvatore Filippone    University of Rome Tor Vergata
diff --git a/README b/README
index 22a4e996..5c00929b 100644
--- a/README
+++ b/README
@@ -1,4 +1,9 @@
-This directory contains  the PSBLAS library, version 3.2.
+This directory contains  the PSBLAS library, version 3.3.
+
+
+Existing codes will work untouched with version 3.3, but they will
+have to be recompiled; new methods have been added for the sake of GPU
+computations, plus the usual bunch of bug fixes.
 
 
 This version requires a working Fortran 2003 compiler; we do not use
diff --git a/test/fileread/df_sample.f90 b/test/fileread/df_sample.f90
index c182d2b1..993422cf 100644
--- a/test/fileread/df_sample.f90
+++ b/test/fileread/df_sample.f90
@@ -142,7 +142,7 @@ program df_sample
     
     m_problem = aux_a%get_nrows()
     call psb_bcast(ictxt,m_problem)
-    call psb_mat_renum(psb_mat_renum_identity_,aux_a,info,perm) 
+    call psb_mat_renum(psb_mat_renum_gps_,aux_a,info,perm) 
 
     ! At this point aux_b may still be unallocated
     if (size(aux_b,dim=1) == m_problem) then
@@ -151,6 +151,10 @@ program df_sample
       b_col_glob =>aux_b(:,1)
       call psb_gelp('N',perm(1:m_problem),&
            & b_col_glob(1:m_problem),info)
+      write(fnout,'(a,i3.3,a)') 'amat-',iam,'-gps.mtx'
+      call aux_a%print(fname=fnout)
+      write(fnout,'(a,i3.3,a)') 'rhs-',iam,'-gps.mtx'
+      call mm_array_write(b_col_glob(1:m_problem),'GPS RHS',info,filename=fnout)
     else
       write(psb_out_unit,'("Generating an rhs...")')
       write(psb_out_unit,'(" ")')
diff --git a/test/kernel/runs/spmv.inp b/test/kernel/runs/spmv.inp
index 3b4e9323..91844ce2 100644
--- a/test/kernel/runs/spmv.inp
+++ b/test/kernel/runs/spmv.inp
@@ -1,3 +1,3 @@
 CSR 
-10
+50
 
diff --git a/test/pargen/runs/ppde.inp b/test/pargen/runs/ppde.inp
index 1a8722e1..87d51541 100644
--- a/test/pargen/runs/ppde.inp
+++ b/test/pargen/runs/ppde.inp
@@ -2,7 +2,7 @@
 BICGSTAB               Iterative method BICGSTAB CGS  BICG BICGSTABL RGMRES
 BJAC                   Preconditioner NONE  DIAG  BJAC 
 CSR                    Storage format for matrix A:  CSR COO JAD 
-050                    Domain size (acutal system is this**3)
+160                    Domain size (acutal system is this**3)
 2                      Stopping criterion
 1000                   MAXIT
 -1                     ITRACE