diff --git a/slv.mk b/slv.mk index 3a33c72..1e2ebe1 100644 --- a/slv.mk +++ b/slv.mk @@ -61,7 +61,7 @@ FINC = $(foreach i,$(INCDIR),-I$(i)/include) -I./ $(MKLINC) $(OTHERINCDIR) FLIB = $(foreach i,$(LIBDIR),-L$(i)/lib) $(PBSDIR) $(OTHERLIBDIR) FLIB += $(MKLLIB) $(METISLIB) $(CGNSLIB) $(HDF5LIB) $(SZIPLIB) $(PBSLIB) $(OTHERLIB) -FCBASE = $(REALLOCFLAG) $(TRACEFLAG) -g +FCBASE = $(REALLOCFLAG) $(TRACEFLAG) ADDFLAGS = $(MKLFLAG) $(MOD) $(FINC) $(FCADD) $(CMPFLAGS) diff --git a/typesgen_mod.f90 b/typesgen_mod.f90 index d141546..61460ec 100644 --- a/typesgen_mod.f90 +++ b/typesgen_mod.f90 @@ -834,6 +834,7 @@ pure function dot_mat_sparse(this,idx,array) result(return_value) do n = 1,size(array,dim=2) #ifdef USE_INTEL_MKL return_value(n) = doti( this%row_val(n1:n2) , & + this%row_idx(n1:n2) , & array( this%row_idx(n1:n2) , n ) ) #else return_value(n) = dot( this%row_val(n1:n2) , & @@ -1255,10 +1256,11 @@ pure function dot_mat2_sparse(this,idx,array) result(return_value) do n = 1,size(array,dim=2) #ifdef USE_INTEL_MKL return_value(n,m) = doti( this%row_val(n1:n2) , & - array( this%row_idx(n1:n2) , n,m ) ) + this%row_idx(n1:n2) , & + array( this%row_idx(n1:n2) , n,m ) ) #else return_value(n,m) = dot( this%row_val(n1:n2) , & - array( this%row_idx(n1:n2) , n,m ) ) + array( this%row_idx(n1:n2) , n,m ) ) #endif end do end do