base/internals/psi_crea_index.f90
 base/serial/psb_znumbmm.f90

Silly bug in psb_znumbmm calling cnumbmm. 
crea_index was trying to get the size of an unallocated array.
psblas3-type-indexed
Salvatore Filippone 14 years ago
parent d4d408f950
commit 8360e7478e

@ -132,8 +132,7 @@ subroutine psi_crea_index(desc_a,index_in,index_out,glob_idx,nxch,nsnd,nrcv,info
end if
if(debug_level >= psb_debug_inner_)&
& write(debug_unit,*) me,' ',trim(name),': calling psi_desc_index',&
& size(index_out)
& write(debug_unit,*) me,' ',trim(name),': calling psi_desc_index'
! Do the actual format conversion.
call psi_desc_index(desc_a,index_in,dep_list(1:,me),&
& length_dl(me),nsnd,nrcv, index_out,glob_idx,info)

@ -169,7 +169,7 @@ contains
mb = b%get_nrows()
nb = b%get_ncols()
call cnumbmm(ma,na,nb,a%irp,a%ja,0,a%val,&
call znumbmm(ma,na,nb,a%irp,a%ja,0,a%val,&
& b%irp,b%ja,0,b%val,&
& c%irp,c%ja,0,c%val,temp)

Loading…
Cancel
Save