diff --git a/fvn_sparse/fvn_sparse.f90 b/fvn_sparse/fvn_sparse.f90 index cc7af45..ea60eaf 100644 --- a/fvn_sparse/fvn_sparse.f90 +++ b/fvn_sparse/fvn_sparse.f90 @@ -112,7 +112,7 @@ allocate(Ap(n+1),Ai(nz)) call umfpack_zl_triplet_to_col(n,n,nz,wTi,wTj,Tx,Tz,Ap,Ai,Ax,Az,status) ! if status is not zero a problem has occured if (status /= 0) then - write(*,*) "Problem during umfpack_zl_triplet_to_col" + write(*,*) "Problem during umfpack_zl_triplet_to_col : ",trim(umfpack_return_code(int(status,kind=sp_kind))) endif ! Define defaults control values @@ -122,7 +122,7 @@ call umfpack_zl_defaults(control) call umfpack_zl_symbolic(n,n,Ap,Ai,Ax,Az,symbolic, control, info) ! info(1) should be zero if (info(1) /= 0) then - write(*,*) "Problem during symbolic analysis : ",trim(umfpack_return_code(int(info(1)))) + write(*,*) "Problem during symbolic analysis : ",trim(umfpack_return_code(int(info(1),kind=sp_kind))) status=info(1) endif @@ -130,7 +130,7 @@ endif call umfpack_zl_numeric (Ap, Ai, Ax, Az, symbolic, numeric, control, info) ! info(1) should be zero if (info(1) /= 0) then - write(*,*) "Problem during numerical factorization : ",trim(umfpack_return_code(int(info(1)))) + write(*,*) "Problem during numerical factorization : ",trim(umfpack_return_code(int(info(1),kind=sp_kind))) status=info(1) endif @@ -143,7 +143,7 @@ if (present(det) ) then ! info(1) should be zero if (info(1) /= 0) then if ( (info(1) < 1) .or. (info(1) >3) ) then ! not a warning - write(*,*) "Problem during sparse determinant : ",trim(umfpack_return_code(int(info(1)))) + write(*,*) "Problem during sparse determinant : ",trim(umfpack_return_code(int(info(1),kind=sp_kind))) endif status=info(1) endif @@ -161,7 +161,7 @@ sys=0 call umfpack_zl_solve (sys, Ap, Ai, Ax,Az, xx,xz, bx,bz, numeric, control, info) ! info(1) should be zero if (info(1) /= 0) then - write(*,*) "Problem during solving : ",trim(umfpack_return_code(int(info(1)))) + write(*,*) "Problem during solving : ",trim(umfpack_return_code(int(info(1),kind=sp_kind))) status=info(1) endif @@ -221,7 +221,7 @@ allocate(Ap(n+1),Ai(nz)) call umfpack_zi_triplet_to_col(n,n,nz,wTi,wTj,Tx,Tz,Ap,Ai,Ax,Az,status) ! if status is not zero a problem has occured if (status /= 0) then - write(*,*) "Problem during umfpack_zl_triplet_to_col" + write(*,*) "Problem during umfpack_zl_triplet_to_col : ",trim(umfpack_return_code(status)) endif ! Define defaults control values @@ -231,7 +231,7 @@ call umfpack_zi_defaults(control) call umfpack_zi_symbolic(n,n,Ap,Ai,Ax,Az,symbolic, control, info) ! info(1) should be zero if (info(1) /= 0) then - write(*,*) "Problem during symbolic analysis : ",trim(umfpack_return_code(int(info(1)))) + write(*,*) "Problem during symbolic analysis : ",trim(umfpack_return_code(int(info(1),kind=sp_kind))) status=info(1) endif @@ -239,7 +239,7 @@ endif call umfpack_zi_numeric (Ap, Ai, Ax, Az, symbolic, numeric, control, info) ! info(1) should be zero if (info(1) /= 0) then - write(*,*) "Problem during numerical factorization : ",trim(umfpack_return_code(int(info(1)))) + write(*,*) "Problem during numerical factorization : ",trim(umfpack_return_code(int(info(1),kind=sp_kind))) status=info(1) endif @@ -252,7 +252,7 @@ if (present(det) ) then ! info(1) should be zero if (info(1) /= 0) then if ( (info(1) < 1) .or. (info(1) >3) ) then ! not a warning - write(*,*) "Problem during sparse determinant : ",trim(umfpack_return_code(int(info(1)))) + write(*,*) "Problem during sparse determinant : ",trim(umfpack_return_code(int(info(1),kind=sp_kind))) endif status=info(1) endif @@ -271,7 +271,7 @@ sys=0 call umfpack_zi_solve (sys, Ap, Ai, Ax,Az, xx,xz, bx,bz, numeric, control, info) ! info(1) should be zero if (info(1) /= 0) then - write(*,*) "Problem during solving : ",trim(umfpack_return_code(int(info(1)))) + write(*,*) "Problem during solving : ",trim(umfpack_return_code(int(info(1),kind=sp_kind))) status=info(1) endif @@ -322,7 +322,7 @@ allocate(Ap(n+1),Ai(nz)) call umfpack_dl_triplet_to_col(n,n,nz,wTi,wTj,T,Ap,Ai,A,status) ! if status is not zero a problem has occured if (status /= 0) then - write(*,*) "Problem during umfpack_dl_triplet_to_col" + write(*,*) "Problem during umfpack_dl_triplet_to_col : ",trim(umfpack_return_code(int(status,kind=sp_kind))) endif ! Define defaults control values @@ -332,7 +332,7 @@ call umfpack_dl_defaults(control) call umfpack_dl_symbolic(n,n,Ap,Ai,A,symbolic, control, info) ! info(1) should be zero if (info(1) /= 0) then - write(*,*) "Problem during symbolic analysis : ",trim(umfpack_return_code(int(info(1)))) + write(*,*) "Problem during symbolic analysis : ",trim(umfpack_return_code(int(info(1),kind=sp_kind))) status=info(1) endif @@ -340,7 +340,7 @@ endif call umfpack_dl_numeric (Ap, Ai, A, symbolic, numeric, control, info) ! info(1) should be zero if (info(1) /= 0) then - write(*,*) "Problem during numerical factorization : ",trim(umfpack_return_code(int(info(1)))) + write(*,*) "Problem during numerical factorization : ",trim(umfpack_return_code(int(info(1),kind=sp_kind))) status=info(1) endif @@ -353,7 +353,7 @@ if (present(det) ) then ! info(1) should be zero if (info(1) /= 0) then if ( (info(1) < 1) .or. (info(1) >3) ) then ! not a warning - write(*,*) "Problem during sparse determinant : ",trim(umfpack_return_code(int(info(1)))) + write(*,*) "Problem during sparse determinant : ",trim(umfpack_return_code(int(info(1),kind=sp_kind))) endif status=info(1) endif @@ -367,7 +367,7 @@ sys=0 call umfpack_dl_solve (sys, Ap, Ai, A, x, B, numeric, control, info) ! info(1) should be zero if (info(1) /= 0) then - write(*,*) "Problem during solving : ",trim(umfpack_return_code(int(info(1)))) + write(*,*) "Problem during solving : ",trim(umfpack_return_code(int(info(1),kind=sp_kind))) status=info(1) endif @@ -417,7 +417,7 @@ allocate(Ap(n+1),Ai(nz)) call umfpack_di_triplet_to_col(n,n,nz,wTi,wTj,T,Ap,Ai,A,status) ! if status is not zero a problem has occured if (status /= 0) then - write(*,*) "Problem during umfpack_di_triplet_to_col" + write(*,*) "Problem during umfpack_di_triplet_to_col : ",trim(umfpack_return_code(status)) endif ! Define defaults control values @@ -427,7 +427,7 @@ call umfpack_di_defaults(control) call umfpack_di_symbolic(n,n,Ap,Ai,A,symbolic, control, info) ! info(1) should be zero if (info(1) /= 0) then - write(*,*) "Problem during symbolic analysis : ",trim(umfpack_return_code(int(info(1)))) + write(*,*) "Problem during symbolic analysis : ",trim(umfpack_return_code(int(info(1),kind=sp_kind))) status=info(1) endif @@ -435,7 +435,7 @@ endif call umfpack_di_numeric (Ap, Ai, A, symbolic, numeric, control, info) ! info(1) should be zero if (info(1) /= 0) then - write(*,*) "Problem during numerical factorization : ",trim(umfpack_return_code(int(info(1)))) + write(*,*) "Problem during numerical factorization : ",trim(umfpack_return_code(int(info(1),kind=sp_kind))) status=info(1) endif @@ -448,7 +448,7 @@ if (present(det) ) then ! info(1) should be zero if (info(1) /= 0) then if ( (info(1) < 1) .or. (info(1) >3) ) then ! not a warning - write(*,*) "Problem during sparse determinant : ",trim(umfpack_return_code(int(info(1)))) + write(*,*) "Problem during sparse determinant : ",trim(umfpack_return_code(int(info(1),kind=sp_kind))) endif status=info(1) endif @@ -460,7 +460,7 @@ sys=0 call umfpack_di_solve (sys, Ap, Ai, A, x, B, numeric, control, info) ! info(1) should be zero if (info(1) /= 0) then - write(*,*) "Problem during solving : ",trim(umfpack_return_code(int(info(1)))) + write(*,*) "Problem during solving : ",trim(umfpack_return_code(int(info(1),kind=sp_kind))) status=info(1) endif @@ -512,7 +512,7 @@ allocate(Ap(n+1),Ai(nz)) call umfpack_zl_triplet_to_col(n,n,nz,wTi,wTj,Tx,Tz,Ap,Ai,Ax,Az,status) ! if status is not zero a problem has occured if (status /= 0) then - write(*,*) "Problem during umfpack_zl_triplet_to_col" + write(*,*) "Problem during umfpack_zl_triplet_to_col ",trim(umfpack_return_code(int(status,kind=sp_kind))) endif ! Define defaults control values @@ -522,7 +522,7 @@ call umfpack_zl_defaults(control) call umfpack_zl_symbolic(n,n,Ap,Ai,Ax,Az,symbolic, control, info) ! info(1) should be zero if (info(1) /= 0) then - write(*,*) "Problem during symbolic analysis : ",trim(umfpack_return_code(int(info(1)))) + write(*,*) "Problem during symbolic analysis : ",trim(umfpack_return_code(int(info(1),kind=sp_kind))) status=info(1) endif @@ -530,7 +530,7 @@ endif call umfpack_zl_numeric (Ap, Ai, Ax, Az, symbolic, numeric, control, info) ! info(1) should be zero if (info(1) /= 0) then - write(*,*) "Problem during numerical factorization : ",trim(umfpack_return_code(int(info(1)))) + write(*,*) "Problem during numerical factorization : ",trim(umfpack_return_code(int(info(1),kind=sp_kind))) status=info(1) endif @@ -541,7 +541,7 @@ call umfpack_zl_get_determinant(det(1),det(2),det(3),numeric,info,status) ! info(1) should be zero if (info(1) /= 0) then if ( (info(1) < 1) .or. (info(1) >3) ) then ! not a warning - write(*,*) "Problem during sparse determinant : ",trim(umfpack_return_code(int(info(1)))) + write(*,*) "Problem during sparse determinant : ",trim(umfpack_return_code(int(info(1),kind=sp_kind))) endif status=info(1) endif @@ -591,7 +591,7 @@ allocate(Ap(n+1),Ai(nz)) call umfpack_zi_triplet_to_col(n,n,nz,wTi,wTj,Tx,Tz,Ap,Ai,Ax,Az,status) ! if status is not zero a problem has occured if (status /= 0) then - write(*,*) "Problem during umfpack_zl_triplet_to_col" + write(*,*) "Problem during umfpack_zl_triplet_to_col : ",trim(umfpack_return_code(status)) endif ! Define defaults control values @@ -601,7 +601,7 @@ call umfpack_zi_defaults(control) call umfpack_zi_symbolic(n,n,Ap,Ai,Ax,Az,symbolic, control, info) ! info(1) should be zero if (info(1) /= 0) then - write(*,*) "Problem during symbolic analysis : ",trim(umfpack_return_code(int(info(1)))) + write(*,*) "Problem during symbolic analysis : ",trim(umfpack_return_code(int(info(1),kind=sp_kind))) status=info(1) endif @@ -609,7 +609,7 @@ endif call umfpack_zi_numeric (Ap, Ai, Ax, Az, symbolic, numeric, control, info) ! info(1) should be zero if (info(1) /= 0) then - write(*,*) "Problem during numerical factorization : ",trim(umfpack_return_code(int(info(1)))) + write(*,*) "Problem during numerical factorization : ",trim(umfpack_return_code(int(info(1),kind=sp_kind))) status=info(1) endif @@ -620,7 +620,7 @@ call umfpack_zi_get_determinant(det(1),det(2),det(3),numeric,info,status) ! info(1) should be zero if (info(1) /= 0) then if ( (info(1) < 1) .or. (info(1) >3) ) then ! not a warning - write(*,*) "Problem during sparse determinant : ",trim(umfpack_return_code(int(info(1)))) + write(*,*) "Problem during sparse determinant : ",trim(umfpack_return_code(int(info(1),kind=sp_kind))) endif status=info(1) endif @@ -661,7 +661,7 @@ allocate(Ap(n+1),Ai(nz)) call umfpack_dl_triplet_to_col(n,n,nz,wTi,wTj,T,Ap,Ai,A,status) ! if status is not zero a problem has occured if (status /= 0) then - write(*,*) "Problem during umfpack_dl_triplet_to_col" + write(*,*) "Problem during umfpack_dl_triplet_to_col : ",trim(umfpack_return_code(int(status,kind=sp_kind))) endif ! Define defaults control values @@ -671,7 +671,7 @@ call umfpack_dl_defaults(control) call umfpack_dl_symbolic(n,n,Ap,Ai,A,symbolic, control, info) ! info(1) should be zero if (info(1) /= 0) then - write(*,*) "Problem during symbolic analysis : ",trim(umfpack_return_code(int(info(1)))) + write(*,*) "Problem during symbolic analysis : ",trim(umfpack_return_code(int(info(1),kind=sp_kind))) status=info(1) endif @@ -679,7 +679,7 @@ endif call umfpack_dl_numeric (Ap, Ai, A, symbolic, numeric, control, info) ! info(1) should be zero if (info(1) /= 0) then - write(*,*) "Problem during numerical factorization : ",trim(umfpack_return_code(int(info(1)))) + write(*,*) "Problem during numerical factorization : ",trim(umfpack_return_code(int(info(1),kind=sp_kind))) status=info(1) endif @@ -690,7 +690,7 @@ call umfpack_dl_get_determinant(det(1),det(2),numeric,info,status) ! info(1) should be zero if (info(1) /= 0) then if ( (info(1) < 1) .or. (info(1) >3) ) then ! not a warning - write(*,*) "Problem during sparse determinant : ",trim(umfpack_return_code(int(info(1)))) + write(*,*) "Problem during sparse determinant : ",trim(umfpack_return_code(int(info(1),kind=sp_kind))) endif status=info(1) endif @@ -736,7 +736,7 @@ allocate(Ap(n+1),Ai(nz)) call umfpack_di_triplet_to_col(n,n,nz,wTi,wTj,T,Ap,Ai,A,status) ! if status is not zero a problem has occured if (status /= 0) then - write(*,*) "Problem during umfpack_di_triplet_to_col" + write(*,*) "Problem during umfpack_di_triplet_to_col : ",trim(umfpack_return_code(status)) endif ! Define defaults control values @@ -746,7 +746,7 @@ call umfpack_di_defaults(control) call umfpack_di_symbolic(n,n,Ap,Ai,A,symbolic, control, info) ! info(1) should be zero if (info(1) /= 0) then - write(*,*) "Problem during symbolic analysis : ",trim(umfpack_return_code(int(info(1)))) + write(*,*) "Problem during symbolic analysis : ",trim(umfpack_return_code(int(info(1),kind=sp_kind))) status=info(1) endif @@ -754,7 +754,7 @@ endif call umfpack_di_numeric (Ap, Ai, A, symbolic, numeric, control, info) ! info(1) should be zero if (info(1) /= 0) then - write(*,*) "Problem during numerical factorization : ",trim(umfpack_return_code(int(info(1)))) + write(*,*) "Problem during numerical factorization : ",trim(umfpack_return_code(int(info(1),kind=sp_kind))) status=info(1) endif @@ -766,7 +766,7 @@ call umfpack_di_get_determinant(det(1),det(2),numeric,info,status) ! info(1) should be zero if (info(1) /= 0) then if ( (info(1) < 1) .or. (info(1) >3) ) then ! not a warning - write(*,*) "Problem during sparse determinant : ",trim(umfpack_return_code(int(info(1)))) + write(*,*) "Problem during sparse determinant : ",trim(umfpack_return_code(int(info(1),kind=sp_kind))) endif status=info(1) endif