[cig-commits] r3870 - short/3D/lithomop/trunk/lithomop3d/liblithomop3d

[willic3 at geodynamics.org]

Author: willic3
Date: 2006-06-23 10:10:55 -0700 (Fri, 23 Jun 2006)
New Revision: 3870

Modified:
   short/3D/lithomop/trunk/lithomop3d/liblithomop3d/formf_ss.f
   short/3D/lithomop/trunk/lithomop3d/liblithomop3d/iterate.F
   short/3D/lithomop/trunk/lithomop3d/liblithomop3d/residu.f
Log:
Commented out some debugging statements.


Modified: short/3D/lithomop/trunk/lithomop3d/liblithomop3d/formf_ss.f
===================================================================
--- short/3D/lithomop/trunk/lithomop3d/liblithomop3d/formf_ss.f	2006-06-23 15:54:47 UTC (rev 3869)
+++ short/3D/lithomop/trunk/lithomop3d/liblithomop3d/formf_ss.f	2006-06-23 17:10:55 UTC (rev 3870)
@@ -76,6 +76,7 @@
 c
       integer i,ielg
       double precision p(60),dl(60)
+cdebug      integer idb
 c
 cdebug      write(6,*) "Hello from formf_ss_f!"
 c
@@ -102,6 +103,7 @@
 	call dsymv("u",nee,one,s,nee,dl,ione,zero,p,ione)
         call addfor(bintern,p,lm(1,ielg),lmx(1,ielg),neq,nee)
       end do
+cdebug      write(6,*) "bintern:",(bintern(idb),idb=1,neq)
       return
       end
 c

Modified: short/3D/lithomop/trunk/lithomop3d/liblithomop3d/iterate.F
===================================================================
--- short/3D/lithomop/trunk/lithomop3d/liblithomop3d/iterate.F	2006-06-23 15:54:47 UTC (rev 3869)
+++ short/3D/lithomop/trunk/lithomop3d/liblithomop3d/iterate.F	2006-06-23 17:10:55 UTC (rev 3870)
@@ -225,6 +225,7 @@
      &                               ierr)
         call VecCreateSeqWithArray(MPI_COMM_WORLD, neq, dispvec,sol,
      &                               ierr)
+cdebug        call VecView(rhs, PETSC_VIEWER_STDOUT_WORLD, ierr)
         call KSPCreate(MPI_COMM_WORLD, ksp, ierr)
         call KSPSetfromOptions(ksp, ierr)
         if (used) then
@@ -233,6 +234,7 @@
         end if
         call KSPSetOperators(ksp,A,A,DIFFERENT_NONZERO_PATTERN,ierr)
         call KSPSolve(ksp, rhs, sol, ierr)
+cdebug        call VecView(sol, PETSC_VIEWER_STDOUT_WORLD, ierr)
         call KSPDestroy(ksp, ierr)
         call VecDestroy(rhs, ierr)
         call VecDestroy(sol, ierr)

Modified: short/3D/lithomop/trunk/lithomop3d/liblithomop3d/residu.f
===================================================================
--- short/3D/lithomop/trunk/lithomop3d/liblithomop3d/residu.f	2006-06-23 15:54:47 UTC (rev 3869)
+++ short/3D/lithomop/trunk/lithomop3d/liblithomop3d/residu.f	2006-06-23 17:10:55 UTC (rev 3870)
@@ -84,12 +84,15 @@
         en=dispvec(i)*bresid(i)
         tmp(3)=tmp(3)+en*en
       end do
+cdebug      write(6,*) "Initial bresid:",(bresid(idb),idb=1,neq)
+cdebug      call flush(6)
       if(nextflag.eq.izero) then
         call fill(bresid,zero,neq)
       else
         call dcopy(neq,bextern,ione,bresid,ione)
       end if
       call daxpy(neq,alpha,bintern,ione,bresid,ione)
+cdebug      write(6,*) "Final bresid:",(bresid(idb),idb=1,neq)
       tmp(1)=dnrm2(neq,dispvec,ione)
       tmp(2)=dnrm2(neq,bresid,ione)
       tmp(3)=sqrt(tmp(3))