From b35062cf29f1154dcf6788353c5a739a383e4c9d Mon Sep 17 00:00:00 2001 From: William Dawson Date: Thu, 18 Apr 2024 11:42:06 +0900 Subject: [PATCH] Seems to be working --- Source/Fortran/FermiOperatorModule.F90 | 1 + Source/Fortran/distributed_includes/FilterMatrix.f90 | 12 ++---------- 2 files changed, 3 insertions(+), 10 deletions(-) diff --git a/Source/Fortran/FermiOperatorModule.F90 b/Source/Fortran/FermiOperatorModule.F90 index 768478a4..b8f88a03 100644 --- a/Source/Fortran/FermiOperatorModule.F90 +++ b/Source/Fortran/FermiOperatorModule.F90 @@ -185,6 +185,7 @@ SUBROUTINE ComputeDenseFOE(H, ISQ, trace, K, inv_temp_in, & !! Scale the eigenvectors CALL MatrixDiagonalScale(vecs, gathered_list) + CALL FilterMatrix(vecs, params%threshold) !! Multiply Back Together CALL TransposeMatrix(vecs, vecsT) diff --git a/Source/Fortran/distributed_includes/FilterMatrix.f90 b/Source/Fortran/distributed_includes/FilterMatrix.f90 index 0373f513..2530b267 100644 --- a/Source/Fortran/distributed_includes/FilterMatrix.f90 +++ b/Source/Fortran/distributed_includes/FilterMatrix.f90 @@ -1,7 +1,5 @@ !! Local Data - TYPE(ProcessGrid_t) :: grid_temp - LOGICAL :: is_complex_temp - INTEGER :: II, size_temp + INTEGER :: II CALL GetMatrixTripletList(this, tlist) CALL ConstructTripletList(new_list) @@ -13,10 +11,4 @@ END IF END DO - size_temp = this%actual_matrix_dimension - grid_temp = this%process_grid - is_complex_temp = this%is_complex - CALL DestructMatrix(this) - CALL ConstructEmptyMatrix(this, size_temp, grid_temp, is_complex_temp) - - CALL FillMatrixFromTripletList(this, new_list) + CALL FillMatrixFromTripletList(this, new_list, preduplicated_in=.TRUE.)