230908.100959.HKT remove unwanted writings

libprima · Sep 8, 2023 · eb3d198 · eb3d198
1 parent 125a5c5
commit eb3d198
Show file tree

Hide file tree

Showing 3 changed files with 3 additions and 47 deletions.
diff --git a/fortran/lincoa/getact.f90 b/fortran/lincoa/getact.f90
@@ -12,7 +12,7 @@ module getact_mod
 !
 ! Started: February 2022
 !
-! Last Modified: Thursday, September 07, 2023 PM12:01:27
+! Last Modified: Friday, September 08, 2023 AM10:09:46
 !--------------------------------------------------------------------------------------------------!
 
 implicit none
@@ -392,8 +392,6 @@ subroutine getact(amat, delta, g, iact, nact, qfac, resact, resnew, rfac, psd)
  call assert(all(is_finite(psd)) .or. nact == 0, 'PSD is finite unless NACT == 0', srname)
  ! In theory, ||PSD||^2 <= GG and -GG <= PSD^T*G <= 0.
  ! N.B. 1. Do not use DD, which may not be up to date. 2. PSD^T*G can be NaN if G is huge.
- write (*, *) 'GG = ', GG, 'PSD*PSD = ', inprod(psd, psd), 'PSD*G = ', inprod(psd, g)
- write (*, *) 'PSD = ', psd
  call assert(inprod(psd, psd) <= TWO * gg, '||PSD||^2 <= 2*GG', srname)
  call assert(.not. (inprod(psd, g) > 1.0E2_RP * EPS * gg .or. inprod(psd, g) < -TWO * gg), '-2*GG <= PSD^T*G <= 0', srname)
 end if

diff --git a/fortran/tests/testsuite/freeform.m b/fortran/tests/testsuite/freeform.m
diff --git a/matlab/setup_tools/freeform.m b/matlab/setup_tools/freeform.m
@@ -8,6 +8,8 @@ function freeform(files)
  files = {listing.name};
 elseif ischarstr(files)
  files = {files};
+else
+ error('Invalid input.')
 end
 
 for ifile = 1 : length(files)