psblas3-type-indexed
Salvatore Filippone 17 years ago
parent 0c4a981595
commit 9a98261efb

@ -26,7 +26,8 @@ lib: $(OBJS)
/bin/cp -p $(HERE)/$(LIBNAME) $(LIBDIR)
/bin/cp -p $(LIBMOD) $(LOCAL_MODS) $(LIBDIR)
$(F90OBJS): $(MODOBJS)
$(OBJS): $(LIBDIR)/psb_base_mod$(.mod)
$(F90OBJS): psb_prec_mod$(.mod)
psb_prec_mod.o: psb_prec_type.o
veryclean: clean

@ -113,7 +113,7 @@ subroutine psb_zbjac_aply(alpha,prec,x,beta,y,desc_data,trans,work,info)
& trans=trans_,unit='U',choice=psb_none_, work=aux)
if(info /=0) goto 9999
case('T','C')
case('T')
call psb_spsm(zone,prec%av(u_pr_),x,zzero,ww,desc_data,info,&
& trans=trans_,unit='L',diag=prec%d,choice=psb_none_, work=aux)
if(info /=0) goto 9999
@ -121,6 +121,14 @@ subroutine psb_zbjac_aply(alpha,prec,x,beta,y,desc_data,trans,work,info)
& trans=trans_,unit='U',choice=psb_none_,work=aux)
if(info /=0) goto 9999
case('C')
call psb_spsm(zone,prec%av(u_pr_),x,zzero,ww,desc_data,info,&
& trans=trans_,unit='L',diag=conjg(prec%d),choice=psb_none_, work=aux)
if(info /=0) goto 9999
call psb_spsm(alpha,prec%av(l_pr_),ww,beta,y,desc_data,info,&
& trans=trans_,unit='U',choice=psb_none_,work=aux)
if(info /=0) goto 9999
end select

Loading…
Cancel
Save