Changed sp2run in examples (#143)

Author: cnegre

examples/sp2run/src/Makefile

@@ -7,7 +7,7 @@ F90      = gfortran
 COMPFLAGS += -fopenmp -I/$(bml_PATH)/install/include -I/$(progress_PATH)/install/include
 
 LINKFLAGS+= -L$(progress_PATH)/install/lib -lprogress
-LINKFLAGS+= -L$(bml_PATH)/install/lib -lbml
+LINKFLAGS+= -L$(bml_PATH)/install/lib  -lbml_fortran -lbml
 LINKFLAGS+= -llapack -lblas
 
 all:

examples/sp2run/src/main.F90

@@ -27,8 +27,19 @@ program sp2run
   call bml_read_matrix(ham_bml,"ham.mtx")
 
   !> Do sp2.
-  call prg_sp2_alg1(ham_bml,rho_bml,sp2%threshold,sp2%bndfil,sp2%minsp2iter,sp2%maxsp2iter &
-      ,sp2%sp2conv,sp2%sp2tol,sp2%verbose)
+  write(*,*)sp2%flavor
+  if(sp2%flavor .eq. "Basic") then 
+    call prg_sp2_basic(ham_bml,rho_bml,sp2%threshold,sp2%bndfil,sp2%minsp2iter,sp2%maxsp2iter &
+       ,sp2%sp2conv,sp2%sp2tol,sp2%verbose)
+  elseif(sp2%flavor .eq. "Alg1") then 
+    call prg_sp2_alg1(ham_bml,rho_bml,sp2%threshold,sp2%bndfil,sp2%minsp2iter,sp2%maxsp2iter &
+       ,sp2%sp2conv,sp2%sp2tol,sp2%verbose)
+  elseif(sp2%flavor .eq. "Alg2") then
+    call prg_sp2_alg2(ham_bml,rho_bml,sp2%threshold,sp2%bndfil,sp2%minsp2iter,sp2%maxsp2iter &
+       ,sp2%sp2conv,sp2%sp2tol,sp2%verbose)
+  else 
+    write(*,*)"SP2 algorithm variant/flavor not recognized"
+  endif 
 
   call bml_print_matrix("rho_bml",rho_bml,0,6,0,6)