psblas3:

base/serial/impl/psb_c_base_mat_impl.f90
 base/serial/impl/psb_d_base_mat_impl.f90
 base/serial/impl/psb_s_base_mat_impl.f90
 base/serial/impl/psb_z_base_mat_impl.f90
 docs/html/footnode.html
 docs/html/img100.png
 docs/html/img101.png
 docs/html/img102.png
 docs/html/img103.png
 docs/html/img104.png
 docs/html/img105.png
 docs/html/img106.png
 docs/html/img107.png
 docs/html/img108.png
 docs/html/img109.png
 docs/html/img110.png
 docs/html/img111.png
 docs/html/img112.png
 docs/html/img113.png
 docs/html/img114.png
 docs/html/img115.png
 docs/html/img116.png
 docs/html/img117.png
 docs/html/img118.png
 docs/html/img119.png
 docs/html/img12.png
 docs/html/img120.png
 docs/html/img121.png
 docs/html/img122.png
 docs/html/img123.png
 docs/html/img124.png
 docs/html/img125.png
 docs/html/img126.png
 docs/html/img127.png
 docs/html/img128.png
 docs/html/img129.png
 docs/html/img13.png
 docs/html/img130.png
 docs/html/img131.png
 docs/html/img132.png
 docs/html/img133.png
 docs/html/img134.png
 docs/html/img135.png
 docs/html/img136.png
 docs/html/img137.png
 docs/html/img138.png
 docs/html/img139.png
 docs/html/img140.png
 docs/html/img141.png
 docs/html/img142.png
 docs/html/img143.png
 docs/html/img144.png
 docs/html/img19.png
 docs/html/img2.png
 docs/html/img20.png
 docs/html/img21.png
 docs/html/img22.png
 docs/html/img23.png
 docs/html/img24.png
 docs/html/img25.png
 docs/html/img26.png
 docs/html/img27.png
 docs/html/img28.png
 docs/html/img29.png
 docs/html/img3.png
 docs/html/img30.png
 docs/html/img31.png
 docs/html/img32.png
 docs/html/img33.png
 docs/html/img34.png
 docs/html/img35.png
 docs/html/img36.png
 docs/html/img37.png
 docs/html/img38.png
 docs/html/img39.png
 docs/html/img40.png
 docs/html/img41.png
 docs/html/img42.png
 docs/html/img43.png
 docs/html/img44.png
 docs/html/img45.png
 docs/html/img46.png
 docs/html/img47.png
 docs/html/img48.png
 docs/html/img49.png
 docs/html/img50.png
 docs/html/img51.png
 docs/html/img52.png
 docs/html/img53.png
 docs/html/img54.png
 docs/html/img55.png
 docs/html/img56.png
 docs/html/img57.png
 docs/html/img58.png
 docs/html/img59.png
 docs/html/img60.png
 docs/html/img61.png
 docs/html/img62.png
 docs/html/img63.png
 docs/html/img64.png
 docs/html/img65.png
 docs/html/img66.png
 docs/html/img67.png
 docs/html/img68.png
 docs/html/img69.png
 docs/html/img70.png
 docs/html/img71.png
 docs/html/img72.png
 docs/html/img73.png
 docs/html/img74.png
 docs/html/img75.png
 docs/html/img76.png
 docs/html/img77.png
 docs/html/img78.png
 docs/html/img79.png
 docs/html/img80.png
 docs/html/img81.png
 docs/html/img82.png
 docs/html/img83.png
 docs/html/img84.png
 docs/html/img85.png
 docs/html/img86.png
 docs/html/img87.png
 docs/html/img88.png
 docs/html/img89.png
 docs/html/img90.png
 docs/html/img91.png
 docs/html/img92.png
 docs/html/img93.png
 docs/html/img94.png
 docs/html/img95.png
 docs/html/img96.png
 docs/html/img97.png
 docs/html/img98.png
 docs/html/img99.png
 docs/html/node100.html
 docs/html/node112.html
 docs/html/node117.html
 docs/html/node119.html
 docs/html/node19.html
 docs/html/node29.html
 docs/html/node35.html
 docs/html/node37.html
 docs/html/node38.html
 docs/html/node39.html
 docs/html/node40.html
 docs/html/node41.html
 docs/html/node42.html
 docs/html/node43.html
 docs/html/node44.html
 docs/html/node45.html
 docs/html/node46.html
 docs/html/node47.html
 docs/html/node48.html
 docs/html/node49.html
 docs/html/node51.html
 docs/html/node52.html
 docs/html/node53.html
 docs/html/node54.html
 docs/html/node56.html
 docs/html/node57.html
 docs/html/node61.html
 docs/html/node62.html
 docs/html/node63.html
 docs/html/node67.html
 docs/html/node68.html
 docs/html/node69.html
 docs/html/node71.html
 docs/html/node72.html
 docs/html/node73.html
 docs/html/node74.html
 docs/html/node75.html
 docs/html/node76.html
 docs/html/node77.html
 docs/html/node80.html
 docs/html/node82.html
 docs/html/node84.html
 docs/html/node85.html
 docs/html/node86.html
 docs/html/node88.html
 docs/html/node9.html
 docs/html/node92.html
 docs/html/node93.html
 docs/html/node94.html
 docs/html/node95.html
 docs/html/node96.html
 docs/html/node97.html
 docs/html/node98.html
 docs/html/node99.html
 docs/psblas-3.0.pdf
 docs/src/datastruct.tex

Provided default implementations of MV_TO_COO and MV_FROM_COO based on
CP.
Fix for documentation.

base/serial/impl/psb_c_base_mat_impl.f90

@@ -86,9 +86,28 @@ subroutine psb_c_base_cp_to_fmt(a,b,info)
   ! Default implementation
   ! 
   info = psb_success_
+  call psb_erractionsave(err_act)
+
   call a%cp_to_coo(tmp,info)
   if (info == psb_success_) call b%mv_from_coo(tmp,info)
 
+  if (info /= psb_success_) then 
+    info = psb_err_from_subroutine_ 
+    call psb_errpush(info,name, a_err='to/from coo')
+    goto 9999
+  end if
+    
+  call psb_erractionrestore(err_act)
+  return
+
+9999 continue
+  call psb_erractionrestore(err_act)
+
+  if (err_act == psb_act_abort_) then
+    call psb_error()
+    return
+  end if
+  return
   return
 
 end subroutine psb_c_base_cp_to_fmt
@@ -111,10 +130,29 @@ subroutine psb_c_base_cp_from_fmt(a,b,info)
   !
   ! Default implementation
   ! 
-  info = psb_success_
+  info  = psb_success_
+  call psb_erractionsave(err_act)
+
   call b%cp_to_coo(tmp,info)
   if (info == psb_success_) call a%mv_from_coo(tmp,info)
 
+  if (info /= psb_success_) then 
+    info = psb_err_from_subroutine_ 
+    call psb_errpush(info,name, a_err='to/from coo')
+    goto 9999
+  end if
+
+  call psb_erractionrestore(err_act)
+  return
+
+
+9999 continue
+  call psb_erractionrestore(err_act)
+
+  if (err_act == psb_act_abort_) then
+    call psb_error()
+    return
+  end if
   return
 
 end subroutine psb_c_base_cp_from_fmt
@@ -134,15 +172,30 @@ subroutine psb_c_base_mv_to_coo(a,b,info)
   character(len=20)  :: name='to_coo'
   logical, parameter :: debug=.false.
 
-  call psb_get_erraction(err_act)
-  ! This is the base version. If we get here
-  ! it means the derived class is incomplete,
-  ! so we throw an error.
-  info = psb_err_missing_override_method_
-  call psb_errpush(info,name,a_err=a%get_fmt())
 
-  if (err_act /= psb_act_ret_) then
+  info  = psb_success_
+  call psb_erractionsave(err_act)
+  
+  call a%cp_to_coo(b,info) 
+
+  if (info /= psb_success_) then 
+    info = psb_err_from_subroutine_ 
+    call psb_errpush(info,name, a_err='to coo')
+    goto 9999
+  end if
+
+  call a%free()
+  
+  call psb_erractionrestore(err_act)
+  return
+
+
+9999 continue
+  call psb_erractionrestore(err_act)
+
+  if (err_act == psb_act_abort_) then
     call psb_error()
+    return
   end if
   return
 
@@ -162,15 +215,29 @@ subroutine psb_c_base_mv_from_coo(a,b,info)
   character(len=20)  :: name='from_coo'
   logical, parameter :: debug=.false.
 
-  call psb_get_erraction(err_act)
+  info  = psb_success_
-  ! This is the base version. If we get here
+  call psb_erractionsave(err_act)
-  ! it means the derived class is incomplete,
+  
-  ! so we throw an error.
+  call a%cp_from_coo(b,info) 
-  info = psb_err_missing_override_method_
-  call psb_errpush(info,name,a_err=a%get_fmt())
 
-  if (err_act /= psb_act_ret_) then
+  if (info /= psb_success_) then 
+    info = psb_err_from_subroutine_ 
+    call psb_errpush(info,name, a_err='from coo')
+    goto 9999
+  end if
+
+  call b%free()
+  
+  call psb_erractionrestore(err_act)
+  return
+
+
+9999 continue
+  call psb_erractionrestore(err_act)
+
+  if (err_act == psb_act_abort_) then
     call psb_error()
+    return
   end if
   return
 
@@ -1550,11 +1617,8 @@ subroutine psb_c_base_inner_vect_sv(alpha,a,x,beta,y,info,trans)
   character(len=20)  :: name='c_base_inner_vect_sv'
   logical, parameter :: debug=.false.
 
-  call psb_get_erraction(err_act)
+  info  = psb_success_
-  ! This is the base version. If we get here
+  call psb_erractionsave(err_act)
-  ! it means the derived class is incomplete,
-  ! so we throw an error.
-  info = psb_success_
 
   call a%inner_cssm(alpha,x%v,beta,y%v,info,trans)  
 

base/serial/impl/psb_d_base_mat_impl.f90

@@ -86,9 +86,28 @@ subroutine psb_d_base_cp_to_fmt(a,b,info)
   ! Default implementation
   ! 
   info = psb_success_
+  call psb_erractionsave(err_act)
+
   call a%cp_to_coo(tmp,info)
   if (info == psb_success_) call b%mv_from_coo(tmp,info)
 
+  if (info /= psb_success_) then 
+    info = psb_err_from_subroutine_ 
+    call psb_errpush(info,name, a_err='to/from coo')
+    goto 9999
+  end if
+    
+  call psb_erractionrestore(err_act)
+  return
+
+9999 continue
+  call psb_erractionrestore(err_act)
+
+  if (err_act == psb_act_abort_) then
+    call psb_error()
+    return
+  end if
+  return
   return
 
 end subroutine psb_d_base_cp_to_fmt
@@ -111,10 +130,29 @@ subroutine psb_d_base_cp_from_fmt(a,b,info)
   !
   ! Default implementation
   ! 
-  info = psb_success_
+  info  = psb_success_
+  call psb_erractionsave(err_act)
+
   call b%cp_to_coo(tmp,info)
   if (info == psb_success_) call a%mv_from_coo(tmp,info)
 
+  if (info /= psb_success_) then 
+    info = psb_err_from_subroutine_ 
+    call psb_errpush(info,name, a_err='to/from coo')
+    goto 9999
+  end if
+
+  call psb_erractionrestore(err_act)
+  return
+
+
+9999 continue
+  call psb_erractionrestore(err_act)
+
+  if (err_act == psb_act_abort_) then
+    call psb_error()
+    return
+  end if
   return
 
 end subroutine psb_d_base_cp_from_fmt
@@ -134,15 +172,30 @@ subroutine psb_d_base_mv_to_coo(a,b,info)
   character(len=20)  :: name='to_coo'
   logical, parameter :: debug=.false.
 
-  call psb_get_erraction(err_act)
-  ! This is the base version. If we get here
-  ! it means the derived class is incomplete,
-  ! so we throw an error.
-  info = psb_err_missing_override_method_
-  call psb_errpush(info,name,a_err=a%get_fmt())
 
-  if (err_act /= psb_act_ret_) then
+  info  = psb_success_
+  call psb_erractionsave(err_act)
+  
+  call a%cp_to_coo(b,info) 
+
+  if (info /= psb_success_) then 
+    info = psb_err_from_subroutine_ 
+    call psb_errpush(info,name, a_err='to coo')
+    goto 9999
+  end if
+
+  call a%free()
+  
+  call psb_erractionrestore(err_act)
+  return
+
+
+9999 continue
+  call psb_erractionrestore(err_act)
+
+  if (err_act == psb_act_abort_) then
     call psb_error()
+    return
   end if
   return
 
@@ -162,15 +215,29 @@ subroutine psb_d_base_mv_from_coo(a,b,info)
   character(len=20)  :: name='from_coo'
   logical, parameter :: debug=.false.
 
-  call psb_get_erraction(err_act)
+  info  = psb_success_
-  ! This is the base version. If we get here
+  call psb_erractionsave(err_act)
-  ! it means the derived class is incomplete,
+  
-  ! so we throw an error.
+  call a%cp_from_coo(b,info) 
-  info = psb_err_missing_override_method_
-  call psb_errpush(info,name,a_err=a%get_fmt())
 
-  if (err_act /= psb_act_ret_) then
+  if (info /= psb_success_) then 
+    info = psb_err_from_subroutine_ 
+    call psb_errpush(info,name, a_err='from coo')
+    goto 9999
+  end if
+
+  call b%free()
+  
+  call psb_erractionrestore(err_act)
+  return
+
+
+9999 continue
+  call psb_erractionrestore(err_act)
+
+  if (err_act == psb_act_abort_) then
     call psb_error()
+    return
   end if
   return
 
@@ -1550,11 +1617,8 @@ subroutine psb_d_base_inner_vect_sv(alpha,a,x,beta,y,info,trans)
   character(len=20)  :: name='d_base_inner_vect_sv'
   logical, parameter :: debug=.false.
 
-  call psb_get_erraction(err_act)
+  info  = psb_success_
-  ! This is the base version. If we get here
+  call psb_erractionsave(err_act)
-  ! it means the derived class is incomplete,
-  ! so we throw an error.
-  info = psb_success_
 
   call a%inner_cssm(alpha,x%v,beta,y%v,info,trans)  
 

base/serial/impl/psb_s_base_mat_impl.f90

@@ -86,9 +86,28 @@ subroutine psb_s_base_cp_to_fmt(a,b,info)
   ! Default implementation
   ! 
   info = psb_success_
+  call psb_erractionsave(err_act)
+
   call a%cp_to_coo(tmp,info)
   if (info == psb_success_) call b%mv_from_coo(tmp,info)
 
+  if (info /= psb_success_) then 
+    info = psb_err_from_subroutine_ 
+    call psb_errpush(info,name, a_err='to/from coo')
+    goto 9999
+  end if
+    
+  call psb_erractionrestore(err_act)
+  return
+
+9999 continue
+  call psb_erractionrestore(err_act)
+
+  if (err_act == psb_act_abort_) then
+    call psb_error()
+    return
+  end if
+  return
   return
 
 end subroutine psb_s_base_cp_to_fmt
@@ -111,10 +130,29 @@ subroutine psb_s_base_cp_from_fmt(a,b,info)
   !
   ! Default implementation
   ! 
-  info = psb_success_
+  info  = psb_success_
+  call psb_erractionsave(err_act)
+
   call b%cp_to_coo(tmp,info)
   if (info == psb_success_) call a%mv_from_coo(tmp,info)
 
+  if (info /= psb_success_) then 
+    info = psb_err_from_subroutine_ 
+    call psb_errpush(info,name, a_err='to/from coo')
+    goto 9999
+  end if
+
+  call psb_erractionrestore(err_act)
+  return
+
+
+9999 continue
+  call psb_erractionrestore(err_act)
+
+  if (err_act == psb_act_abort_) then
+    call psb_error()
+    return
+  end if
   return
 
 end subroutine psb_s_base_cp_from_fmt
@@ -134,15 +172,30 @@ subroutine psb_s_base_mv_to_coo(a,b,info)
   character(len=20)  :: name='to_coo'
   logical, parameter :: debug=.false.
 
-  call psb_get_erraction(err_act)
-  ! This is the base version. If we get here
-  ! it means the derived class is incomplete,
-  ! so we throw an error.
-  info = psb_err_missing_override_method_
-  call psb_errpush(info,name,a_err=a%get_fmt())
 
-  if (err_act /= psb_act_ret_) then
+  info  = psb_success_
+  call psb_erractionsave(err_act)
+  
+  call a%cp_to_coo(b,info) 
+
+  if (info /= psb_success_) then 
+    info = psb_err_from_subroutine_ 
+    call psb_errpush(info,name, a_err='to coo')
+    goto 9999
+  end if
+
+  call a%free()
+  
+  call psb_erractionrestore(err_act)
+  return
+
+
+9999 continue
+  call psb_erractionrestore(err_act)
+
+  if (err_act == psb_act_abort_) then
     call psb_error()
+    return
   end if
   return
 
@@ -162,15 +215,29 @@ subroutine psb_s_base_mv_from_coo(a,b,info)
   character(len=20)  :: name='from_coo'
   logical, parameter :: debug=.false.
 
-  call psb_get_erraction(err_act)
+  info  = psb_success_
-  ! This is the base version. If we get here
+  call psb_erractionsave(err_act)
-  ! it means the derived class is incomplete,
+  
-  ! so we throw an error.
+  call a%cp_from_coo(b,info) 
-  info = psb_err_missing_override_method_
-  call psb_errpush(info,name,a_err=a%get_fmt())
 
-  if (err_act /= psb_act_ret_) then
+  if (info /= psb_success_) then 
+    info = psb_err_from_subroutine_ 
+    call psb_errpush(info,name, a_err='from coo')
+    goto 9999
+  end if
+
+  call b%free()
+  
+  call psb_erractionrestore(err_act)
+  return
+
+
+9999 continue
+  call psb_erractionrestore(err_act)
+
+  if (err_act == psb_act_abort_) then
     call psb_error()
+    return
   end if
   return
 
@@ -1550,11 +1617,8 @@ subroutine psb_s_base_inner_vect_sv(alpha,a,x,beta,y,info,trans)
   character(len=20)  :: name='s_base_inner_vect_sv'
   logical, parameter :: debug=.false.
 
-  call psb_get_erraction(err_act)
+  info  = psb_success_
-  ! This is the base version. If we get here
+  call psb_erractionsave(err_act)
-  ! it means the derived class is incomplete,
-  ! so we throw an error.
-  info = psb_success_
 
   call a%inner_cssm(alpha,x%v,beta,y%v,info,trans)  
 

base/serial/impl/psb_z_base_mat_impl.f90

@@ -86,9 +86,28 @@ subroutine psb_z_base_cp_to_fmt(a,b,info)
   ! Default implementation
   ! 
   info = psb_success_
+  call psb_erractionsave(err_act)
+
   call a%cp_to_coo(tmp,info)
   if (info == psb_success_) call b%mv_from_coo(tmp,info)
 
+  if (info /= psb_success_) then 
+    info = psb_err_from_subroutine_ 
+    call psb_errpush(info,name, a_err='to/from coo')
+    goto 9999
+  end if
+    
+  call psb_erractionrestore(err_act)
+  return
+
+9999 continue
+  call psb_erractionrestore(err_act)
+
+  if (err_act == psb_act_abort_) then
+    call psb_error()
+    return
+  end if
+  return
   return
 
 end subroutine psb_z_base_cp_to_fmt
@@ -111,10 +130,29 @@ subroutine psb_z_base_cp_from_fmt(a,b,info)
   !
   ! Default implementation
   ! 
-  info = psb_success_
+  info  = psb_success_
+  call psb_erractionsave(err_act)
+
   call b%cp_to_coo(tmp,info)
   if (info == psb_success_) call a%mv_from_coo(tmp,info)
 
+  if (info /= psb_success_) then 
+    info = psb_err_from_subroutine_ 
+    call psb_errpush(info,name, a_err='to/from coo')
+    goto 9999
+  end if
+
+  call psb_erractionrestore(err_act)
+  return
+
+
+9999 continue
+  call psb_erractionrestore(err_act)
+
+  if (err_act == psb_act_abort_) then
+    call psb_error()
+    return
+  end if
   return
 
 end subroutine psb_z_base_cp_from_fmt
@@ -134,15 +172,30 @@ subroutine psb_z_base_mv_to_coo(a,b,info)
   character(len=20)  :: name='to_coo'
   logical, parameter :: debug=.false.
 
-  call psb_get_erraction(err_act)
-  ! This is the base version. If we get here
-  ! it means the derived class is incomplete,
-  ! so we throw an error.
-  info = psb_err_missing_override_method_
-  call psb_errpush(info,name,a_err=a%get_fmt())
 
-  if (err_act /= psb_act_ret_) then
+  info  = psb_success_
+  call psb_erractionsave(err_act)
+  
+  call a%cp_to_coo(b,info) 
+
+  if (info /= psb_success_) then 
+    info = psb_err_from_subroutine_ 
+    call psb_errpush(info,name, a_err='to coo')
+    goto 9999
+  end if
+
+  call a%free()
+  
+  call psb_erractionrestore(err_act)
+  return
+
+
+9999 continue
+  call psb_erractionrestore(err_act)
+
+  if (err_act == psb_act_abort_) then
     call psb_error()
+    return
   end if
   return
 
@@ -162,15 +215,29 @@ subroutine psb_z_base_mv_from_coo(a,b,info)
   character(len=20)  :: name='from_coo'
   logical, parameter :: debug=.false.
 
-  call psb_get_erraction(err_act)
+  info  = psb_success_
-  ! This is the base version. If we get here
+  call psb_erractionsave(err_act)
-  ! it means the derived class is incomplete,
+  
-  ! so we throw an error.
+  call a%cp_from_coo(b,info) 
-  info = psb_err_missing_override_method_
-  call psb_errpush(info,name,a_err=a%get_fmt())
 
-  if (err_act /= psb_act_ret_) then
+  if (info /= psb_success_) then 
+    info = psb_err_from_subroutine_ 
+    call psb_errpush(info,name, a_err='from coo')
+    goto 9999
+  end if
+
+  call b%free()
+  
+  call psb_erractionrestore(err_act)
+  return
+
+
+9999 continue
+  call psb_erractionrestore(err_act)
+
+  if (err_act == psb_act_abort_) then
     call psb_error()
+    return
   end if
   return
 
@@ -1550,11 +1617,8 @@ subroutine psb_z_base_inner_vect_sv(alpha,a,x,beta,y,info,trans)
   character(len=20)  :: name='z_base_inner_vect_sv'
   logical, parameter :: debug=.false.
 
-  call psb_get_erraction(err_act)
+  info  = psb_success_
-  ! This is the base version. If we get here
+  call psb_erractionsave(err_act)
-  ! it means the derived class is incomplete,
-  ! so we throw an error.
-  info = psb_success_
 
   call a%inner_cssm(alpha,x%v,beta,y%v,info,trans)  
 

docs/html/footnode.html

@@ -104,7 +104,7 @@ sample scatter/gather routines.
 .
 </PRE>
 </DD>
-<DT><A NAME="foot7091">... follows</A><A
+<DT><A NAME="foot7055">... follows</A><A
  HREF="node112.html#tex2html30"><SUP>3</SUP></A></DT>
 <DD>The string is case-insensitive