Merge lp:~nickpapior/siesta/kovalp-fixes into lp:~kovalp/siesta/siesta

[Nick Papior]

There were a few missing items that you have missed.

Please, read the commit log and see what I have changed. I have runned your example with the committed code without problems. But perhaps you should run it again.

[Petr Koval (kovalp)]

review: Approve

[Petr Koval (kovalp)]

It works as it should, thank you. I will try later to update other
Makefiles as well. Thank you.

Peter

On Tue, Nov 22, 2016 at 9:34 AM, <email address hidden> wrote:

> The proposal to merge lp:~nickpapior/siesta/kovalp-fixes into
> lp:~kovalp/siesta/siesta has been updated.
>
> Status: Needs review => Merged
>
> For more details, see:
> https://code.launchpad.net/~nickpapior/siesta/kovalp-fixes/+merge/311468
> --
> You are reviewing the proposed merge of lp:~nickpapior/siesta/kovalp-fixes
> into lp:~kovalp/siesta/siesta.
>

=== modified file 'Src/Makefile'
--- Src/Makefile	2016-11-17 13:54:18 +0000
+++ Src/Makefile	2016-11-22 08:16:45 +0000
@@ -675,9 +675,7 @@
 compute_rhog.o: siesta_options.o sparse_matrices.o
 conjgr.o: precision.o
 conjgr_old.o: precision.o
-coor.o: alloc.o parallel.o precision.o siesta_geom.o sys.o units.o zmatrix.o \
-  m_poisson_fdf.o 
- 
+coor.o: alloc.o parallel.o precision.o siesta_geom.o sys.o units.o zmatrix.o
 coxmol.o: files.o periodic_table.o precision.o
 create_Sparsity_SC.o: class_Sparsity.o geom_helper.o intrinsic_missing.o
 create_Sparsity_Union.o: class_OrbitalDistribution.o class_Sparsity.o
@@ -692,31 +690,14 @@
 detover.o: alloc.o linpack.o parallel.o parallelsubs.o precision.o
 dfscf.o: alloc.o atm_types.o atmfuncs.o atomlist.o listsc.o m_spin.o mesh.o
 dfscf.o: meshdscf.o meshphi.o parallel.o parallelsubs.o precision.o sys.o
-dhscf.o: alloc.o atmfuncs.o bsc_xcmod.o cellxc_mod.o delk.o dfscf.o \
-  m_poisson_fdf.o m_poisson_switch.o
-
-m_poisson_fdf.o : $(FDF) m_char.o
-
-m_poisson_switch.o: poison.o m_poisson_obc_fastconv.o
-   
-m_poisson_obc_fastconv.o: m_plot_1d_arr3d.o m_fastconv_dims.o \
-  m_poisson_obc_fastconv_r2c.o m_mpi_stop.o
-    
-m_poisson_obc_fastconv_r2c.o: m_plot_1d_arr3d.o m_fastconv_dims.o \
-  m_poisson_obc_comm.o  cellsubs.o
-  
-m_poisson_obc_comm.o: atm_types.o siesta_geom.o atmfuncs.o \
-  $(FDF) parallel.o  m_mpi_stop.o
-
-m_plot_1d_arr3d.o: precision.o
-
+dhscf.o: alloc.o atmfuncs.o bsc_xcmod.o cellxc_mod.o delk.o dfscf.o
 dhscf.o: doping_uniform.o files.o forhar.o iogrid_netcdf.o m_charge_add.o
 dhscf.o: m_efield.o m_hartree_add.o m_iorho.o m_iotddft.o m_mesh_node.o
-dhscf.o: m_ncdf_io.o m_ncdf_siesta.o m_partial_charges.o m_rhog.o m_spin.o
-dhscf.o: m_ts_global_vars.o m_ts_hartree.o m_ts_options.o m_ts_tdir.o
-dhscf.o: m_ts_voltage.o mesh.o meshdscf.o meshsubs.o moremeshsubs.o parallel.o
-dhscf.o: parsing.o precision.o rhofft.o rhoofd.o siesta_options.o sys.o units.o
-dhscf.o: vmat.o
+dhscf.o: m_ncdf_io.o m_ncdf_siesta.o m_partial_charges.o m_poisson_fdf.o
+dhscf.o: m_poisson_switch.o m_rhog.o m_spin.o m_ts_global_vars.o m_ts_hartree.o
+dhscf.o: m_ts_options.o m_ts_tdir.o m_ts_voltage.o mesh.o meshdscf.o meshsubs.o
+dhscf.o: moremeshsubs.o parallel.o parsing.o precision.o rhofft.o rhoofd.o
+dhscf.o: siesta_options.o sys.o units.o vmat.o
 diag2g.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
 diag2k.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
 diag3g.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
@@ -797,9 +778,9 @@
 initatom.o: atmparams.o atom.o atom_options.o basis_io.o basis_specs.o
 initatom.o: basis_types.o chemical.o electrostatic.o ldau_specs.o m_spin.o
 initatom.o: old_atmfuncs.o precision.o pseudopotential.o
-initparallel.o: alloc.o class_OrbitalDistribution.o domain_decom.o parallel.o
-initparallel.o: parallelsubs.o precision.o siesta_options.o sparse_matrices.o
-initparallel.o: spatial.o sys.o m_poisson_fdf.o
+initparallel.o: alloc.o class_OrbitalDistribution.o domain_decom.o
+initparallel.o: m_poisson_fdf.o parallel.o parallelsubs.o precision.o
+initparallel.o: siesta_options.o sparse_matrices.o spatial.o sys.o
 io.o: sys.o
 iocg.o: files.o precision.o sys.o
 iodm.o: alloc.o files.o parallel.o parallelsubs.o precision.o sys.o
@@ -860,6 +841,7 @@
 m_eo.o: precision.o
 m_evolve.o: alloc.o m_memory.o parallel.o parallelsubs.o precision.o sys.o
 m_exp_coord.o: files.o m_os.o parallel.o precision.o units.o
+m_fastconv_dims.o: precision.o
 m_filter.o: bessph.o precision.o radfft.o sorting.o sys.o
 m_fire.o: parallel.o precision.o
 m_fire_mixing.o: alloc.o m_fire_para.o m_mpi_utils.o parallel.o precision.o
@@ -950,6 +932,15 @@
 m_planewavematrix.o: atomlist.o dhscf.o m_planewavematrixvar.o precision.o
 m_planewavematrix.o: siesta_geom.o sparse_matrices.o
 m_planewavematrixvar.o: precision.o
+m_plot_1d_arr3d.o: precision.o
+m_poisson_fdf.o: m_char.o precision.o
+m_poisson_obc_comm.o: atm_types.o atmfuncs.o m_mpi_stop.o parallel.o
+m_poisson_obc_comm.o: precision.o siesta_geom.o sys.o
+m_poisson_obc_fastconv.o: cellsubs.o m_mpi_stop.o m_poisson_obc_fastconv_r2c.o
+m_poisson_obc_fastconv.o: parallel.o precision.o sys.o
+m_poisson_obc_fastconv_r2c.o: cellsubs.o m_fastconv_dims.o m_plot_1d_arr3d.o
+m_poisson_obc_fastconv_r2c.o: m_poisson_obc_comm.o precision.o sys.o
+m_poisson_switch.o: m_poisson_obc_fastconv.o precision.o sys.o
 m_pulay.o: alloc.o atomlist.o m_mpi_utils.o m_spin.o m_svd.o parallel.o
 m_pulay.o: precision.o siesta_options.o sparse_matrices.o sys.o
 m_redist_spmatrix.o: alloc.o class_Distribution.o
@@ -1312,12 +1303,12 @@
 siesta_init.o: class_Fstack_Pair_Geometry_SpData2D.o densematrix.o files.o
 siesta_init.o: flook_siesta.o ioxv.o kpoint_grid.o kpoint_pdos.o ksvinit.o
 siesta_init.o: m_cite.o m_energies.o m_eo.o m_fixed.o m_forces.o m_gamma.o
-siesta_init.o: m_iostruct.o m_mpi_utils.o m_new_dm.o m_rmaxh.o m_spin.o
-siesta_init.o: m_steps.o m_supercell.o m_timer.o m_wallclock.o metaforce.o
-siesta_init.o: molecularmechanics.o object_debug.o parallel.o parallelsubs.o
-siesta_init.o: projected_DOS.o siesta_cmlsubs.o siesta_dicts.o siesta_geom.o
-siesta_init.o: siesta_options.o sparse_matrices.o struct_init.o sys.o timer.o
-siesta_init.o: timestamp.o ts_init.o writewave.o zmatrix.o
+siesta_init.o: m_iostruct.o m_mpi_utils.o m_new_dm.o m_poisson_fdf.o m_rmaxh.o
+siesta_init.o: m_spin.o m_steps.o m_supercell.o m_timer.o m_wallclock.o
+siesta_init.o: metaforce.o molecularmechanics.o object_debug.o parallel.o
+siesta_init.o: parallelsubs.o projected_DOS.o siesta_cmlsubs.o siesta_dicts.o
+siesta_init.o: siesta_geom.o siesta_options.o sparse_matrices.o struct_init.o
+siesta_init.o: sys.o timer.o timestamp.o ts_init.o writewave.o zmatrix.o
 siesta_master.o: iopipes.o iosockets.o precision.o sys.o
 siesta_move.o: atomlist.o broyden_optim.o cell_broyden_optim.o
 siesta_move.o: cell_fire_optim.o dynamics.o fire_optim.o flook_siesta.o ioxv.o
@@ -1348,11 +1339,10 @@
 state_init.o: m_ts_tri_init.o normalize_dm.o overlap.o parallel.o
 state_init.o: proximity_check.o siesta_cml.o siesta_geom.o siesta_options.o
 state_init.o: sparse_matrices.o sys.o units.o write_subs.o zmatrix.o
-struct_init.o: alloc.o atmfuncs.o atomlist.o files.o ioxv.o m_exp_coord.o \
-  m_poisson_fdf.o m_poisson_obc_comm.o
-struct_init.o: m_iostruct.o m_mpi_utils.o m_steps.o parallel.o periodic_table.o
-struct_init.o: siesta_cml.o siesta_geom.o siesta_master.o siesta_options.o
-struct_init.o: units.o zmatrix.o
+struct_init.o: alloc.o atmfuncs.o atomlist.o files.o ioxv.o m_exp_coord.o
+struct_init.o: m_iostruct.o m_mpi_utils.o m_poisson_fdf.o m_poisson_obc_comm.o
+struct_init.o: m_steps.o parallel.o periodic_table.o siesta_cml.o siesta_geom.o
+struct_init.o: siesta_master.o siesta_options.o units.o zmatrix.o
 sys.o: parallel.o
 timer.o: extrae_eventllist.o extrae_module.o m_timer.o parallel.o sys.o
 timer.o: timer_tree.o

=== modified file 'Src/elpa2_kernels.f90'
--- Src/elpa2_kernels.f90	2016-10-27 08:15:19 +0000
+++ Src/elpa2_kernels.f90	2016-11-22 08:16:45 +0000
@@ -18,7 +18,7 @@
    integer, parameter :: dprec = kind(1.d0)
 
    integer, intent(in) :: nb, nq, ldq, ldh
-   real(dprec), intent(inout) :: q(ldq,*) 
+   real(dprec), intent(inout) :: q(ldq,*)
    real(dprec), intent(in) :: hh(ldh,*)
 
    real(dprec) s
@@ -51,7 +51,7 @@
       call hh_trafo_kernel_4(q(i,1),hh, nb, ldq, ldh, s)
    endif
 
-end subroutine
+end
 
 ! --------------------------------------------------------------------------------------------------
 ! The following kernels perform the Householder transformation on Q for 12/8/4 rows.
@@ -572,4 +572,5 @@
    enddo
 
 end
+
 ! --------------------------------------------------------------------------------------------------

=== modified file 'Src/m_fastconv_dims.F90'
--- Src/m_fastconv_dims.F90	2016-11-16 09:58:12 +0000
+++ Src/m_fastconv_dims.F90	2016-11-22 08:16:45 +0000
@@ -9,12 +9,6 @@
   use precision, only :  dp
   
   implicit none
-  
-  complex(dp), private, allocatable :: coul_r2c(:,:,:)
-  integer(8), private :: p_r2c, p_c2r
-  real(dp), allocatable, private    :: io_dreal(:,:,:)
-  complex(dp), allocatable, private :: dens_r2c(:,:,:)
-  integer, private :: icall = 0
 
 contains
 
@@ -22,7 +16,7 @@
 ! Computes 6 numbers out of limits for a general 1d convolution
 !
 subroutine conv_lims(s1,f1,s2,f2,s3,f3,n1,n2,ncaux,saux,salloc_aux,falloc_aux)
-  implicit none
+
   ! external
   integer, intent(in) :: s1, f1, s2, f2, s3, f3 ! start/finish indices for functions in a convolution  f3 = f1*f2
   integer, intent(out) :: n1,n2,ncaux,saux,salloc_aux,falloc_aux
@@ -40,7 +34,6 @@
 ! Calculate the limits for FFTW r2c transform
 !
 subroutine conv_calc_r2c_dims(nff_in_1, nff_in_2, nff_complex, nff_real)
-  implicit none
   integer, intent(in)  :: nff_in_1, nff_in_2
   integer, intent(out) :: nff_complex, nff_real
   nff_real = (nff_in_1+nff_in_2)

=== modified file 'Src/m_mpi_stop.F90'
--- Src/m_mpi_stop.F90	2016-11-16 13:02:00 +0000
+++ Src/m_mpi_stop.F90	2016-11-22 08:16:45 +0000
@@ -6,15 +6,12 @@
 ! See Docs/Contributors.txt for a list of contributors.
 !
 module m_mpi_stop
+
   implicit none
-
+  
 contains
 
-!
-!
-!
 subroutine mpi_stop(cfile, iline, ierr)
-  implicit none
   character(*), intent(in) :: cfile
   integer, intent(in) :: iline
   integer, intent(in) :: ierr

=== modified file 'Src/m_plot_1d_arr3d.F90'
--- Src/m_plot_1d_arr3d.F90	2016-10-18 13:15:38 +0000
+++ Src/m_plot_1d_arr3d.F90	2016-11-22 08:16:45 +0000
@@ -6,6 +6,7 @@
 ! See Docs/Contributors.txt for a list of contributors.
 !
 module m_plot_1d_arr3d
+  
   implicit none
   
 contains
@@ -15,7 +16,6 @@
 ! 
 subroutine plot_1d_arr3d(nnn, arr, ialong,i1,i2, fname)
   use precision, only : grid_p
-  implicit none
   !! external
   integer, intent(in) :: nnn(:) ! 3 dimensions of array arr(*)
   real(grid_p), intent(in) :: arr(*) ! data array
@@ -67,7 +67,6 @@
 !! 
 !!
 integer function get_free_handle()
-  implicit none
   ! internal
   integer :: ihandle
   logical :: logical_opened

=== modified file 'Src/m_poisson_fdf.F90'
--- Src/m_poisson_fdf.F90	2016-10-25 15:07:07 +0000
+++ Src/m_poisson_fdf.F90	2016-11-22 08:16:45 +0000
@@ -12,7 +12,7 @@
   
   implicit none 
   
-  character(80) :: PoissonMethod = ''
+  character(32) :: PoissonMethod = ''
   real(dp)      :: PoissonPBCCutoff0D = -999
   
 contains
@@ -23,16 +23,20 @@
 subroutine poisson_fdf()
   use fdf, only : fdf_get
   use m_char, only : ucase
-  implicit none
 
-  ! internal
-     
   !
   ! Choice of Poison solver  
   !
  
   PoissonMethod = fdf_get('Poisson.Method', 'PBC')
   PoissonMethod = ucase(PoissonMethod)
+  select case (PoissonMethod)
+  case('PBC', 'OBC.FASTCONV')
+  case default
+     write(*,*) 'The Poisson.Method is not recognized, &
+          &should be one of PBC/OBC.FastConv'
+     call die('The Poisson.Method is not recognized')
+  end select
   
   PoissonPBCCutoff0D = fdf_get('Poisson.PBC.Cutoff0D', 1.0_dp)
 

=== modified file 'Src/m_poisson_obc_comm.F90'
--- Src/m_poisson_obc_comm.F90	2016-11-16 13:02:00 +0000
+++ Src/m_poisson_obc_comm.F90	2016-11-22 08:16:45 +0000
@@ -72,7 +72,9 @@
         box(ix,2) = max(xa(ix,ia)+rc, box(ix,2))
       enddo  
     enddo
-    
+
+    deallocate(sp2rcut)
+
     box_c(1:3) = (box(:,2)+box(:,1))/2
     
 !    write(6,'(a)') 'box'
@@ -95,8 +97,6 @@
 !    write(6,*) origin
 !    write(6,'(3g12.5)') xa
   endif
-
-  if(allocated(sp2rcut)) deallocate(sp2rcut)
   
 end subroutine ! obc_coord_shift
 
@@ -184,6 +184,7 @@
   
   j3min = n2s3(node)
   j3max = n2s3(node+1)-1
+  deallocate(n2n3,n2s3)
 
   nnn(1:3) = [j1max-j1min+1,j2max-j2min+1,j3max-j3min+1]
   if(any(nnn/=[n1,n2,n3])) then

=== modified file 'Src/m_poisson_obc_fastconv.F90'
--- Src/m_poisson_obc_fastconv.F90	2016-11-16 13:02:00 +0000
+++ Src/m_poisson_obc_fastconv.F90	2016-11-22 08:16:45 +0000
@@ -106,15 +106,21 @@
       allocate(v_gth(np_gth))                                                             
       call poisson_obc_fastconv_r2c( cell, mesh(1),mesh(2),mesh(3),mesh, &
         rho_gth, U, v_gth, STRESS)
-        
+      
 !      write(6,*) 'sum(v_gth)', sum(v_gth)
 
     endif
-
+    deallocate(rho_gth)
+    
     !! only potential must be scattered  
     call MPI_Scatterv(v_gth, n2c, n2d, MPI_grid_real, v,np, MPI_grid_real, root, comm, ierr)
-    call mpi_stop(__FILE__, __LINE__, ierr)  
-
+    call mpi_stop(__FILE__, __LINE__, ierr)
+    if(root==node) then
+       deallocate(v_gth)
+    end if
+   
+    deallocate(n2c, n2d)
+    
   endif
   
 end subroutine !poisson_obc_fastconv

=== modified file 'Src/m_poisson_obc_fastconv_r2c.F90'
--- Src/m_poisson_obc_fastconv_r2c.F90	2016-11-17 14:11:46 +0000
+++ Src/m_poisson_obc_fastconv_r2c.F90	2016-11-22 08:16:45 +0000
@@ -165,10 +165,10 @@
   do i3=1,nr(3)
   do i2=1,nr(2)
   do i1=1,nc(1)
-    z1 = cmplx(coul_r2c((i1-1)*2+1,i2,i3), coul_r2c((i1-1)*2+2,i2,i3), 8)
-    z2 = cmplx(dens_fft((i1-1)*2+1,i2,i3), dens_fft((i1-1)*2+2,i2,i3), 8)
+    z1 = cmplx(coul_r2c((i1-1)*2+1,i2,i3), coul_r2c((i1-1)*2+2,i2,i3), dp)
+    z2 = cmplx(dens_fft((i1-1)*2+1,i2,i3), dens_fft((i1-1)*2+2,i2,i3), dp)
     z3 = z1*z2
-    dens_fft((i1-1)*2+1,i2,i3) = real(z3,8)              ! Product density and Coulomb interaction in momentum space
+    dens_fft((i1-1)*2+1,i2,i3) = real(z3,dp)             ! Product density and Coulomb interaction in momentum space
     dens_fft((i1-1)*2+2,i2,i3) = aimag(z3)               ! Product density and Coulomb interaction in momentum space
   enddo
   enddo 
@@ -220,7 +220,7 @@
 !
 !
 subroutine comp_coul_r2c_inplace(ic, nn, dcell, n, coul_r2c)
-  use precision,   only : dp
+  use precision,   only : dp, i8b
   use cellsubs,    only : volcel  ! Finds unit cell volume
   implicit none
   !! external
@@ -230,7 +230,7 @@
   integer, intent(in) :: n(:)        ! 3
   real(dp), intent(inout) :: coul_r2c(:,:,:)
   !! internal
-  integer(8) :: p_r2c
+  integer(i8b) :: p_r2c
 
   integer :: i1, i2, i3
   real(dp) :: rp(3), dcell2(3,3), dvol

=== modified file 'Util/COOP/Makefile'
--- Util/COOP/Makefile	2016-08-30 14:02:33 +0000
+++ Util/COOP/Makefile	2016-11-22 08:16:45 +0000
@@ -108,7 +108,7 @@
 bsc_cellxc.o: alloc.o bsc_xcmod.o cellxc_mod.o mesh.o moremeshsubs.o parallel.o
 bsc_cellxc.o: parallelsubs.o precision.o sys_local.o
 bsc_xcmod.o: parallel.o precision.o sys_local.o
-cart2frac.o: sys_local.o
+cart2frac.o: precision.o sys_local.o
 cdiag.o: alloc.o m_diagon_opt.o memoryinfo.o parallel.o precision.o sys_local.o
 cell_broyden_optim.o: m_broyddj_nocomm.o parallel.o precision.o sys_local.o
 cell_broyden_optim.o: units.o zmatrix.o
@@ -121,6 +121,8 @@
 cgvc_zmatrix.o: alloc.o conjgr.o m_mpi_utils.o parallel.o precision.o
 cgvc_zmatrix.o: sys_local.o units.o zmatrix.o
 cgwf.o: egandd.o onmod.o parallel.o precision.o sys_local.o
+changebasis.o: alloc.o m_matdiag.o parallel.o parallelsubs.o precision.o
+changebasis.o: sparse_matrices.o sys_local.o wavefunctions.o
 chemical.o: parallel.o precision.o sys_local.o
 chempot.o: alloc.o mneighb.o parallel.o parallelsubs.o precision.o sys_local.o
 chkdim.o: sys_local.o
@@ -173,17 +175,17 @@
 denmat.o: alloc.o globalise.o m_mpi_utils.o onmod.o precision.o
 denmatlomem.o: alloc.o globalise.o onmod.o precision.o
 densematrix.o: alloc.o precision.o
-detover.o: alloc.o parallel.o parallelsubs.o precision.o
+detover.o: alloc.o linpack.o parallel.o parallelsubs.o precision.o
 dfscf.o: alloc.o atm_types.o atmfuncs.o atomlist.o listsc.o m_spin.o mesh.o
 dfscf.o: meshdscf.o meshphi.o parallel.o parallelsubs.o precision.o sys_local.o
 dhscf.o: alloc.o atmfuncs.o bsc_xcmod.o cellxc_mod.o delk.o dfscf.o
 dhscf.o: doping_uniform.o files.o forhar.o iogrid_netcdf.o m_charge_add.o
-dhscf.o: m_efield.o m_hartree_add.o m_iorho.o m_mesh_node.o m_ncdf_io.o
-dhscf.o: m_ncdf_siesta.o m_partial_charges.o m_rhog.o m_spin.o
-dhscf.o: m_ts_global_vars.o m_ts_hartree.o m_ts_options.o m_ts_tdir.o
-dhscf.o: m_ts_voltage.o mesh.o meshdscf.o meshsubs.o moremeshsubs.o parallel.o
-dhscf.o: parsing.o precision.o rhofft.o rhoofd.o siesta_options.o sys_local.o
-dhscf.o: units.o vmat.o
+dhscf.o: m_efield.o m_hartree_add.o m_iorho.o m_iotddft.o m_mesh_node.o
+dhscf.o: m_ncdf_io.o m_ncdf_siesta.o m_partial_charges.o m_poisson_fdf.o
+dhscf.o: m_poisson_switch.o m_rhog.o m_spin.o m_ts_global_vars.o m_ts_hartree.o
+dhscf.o: m_ts_options.o m_ts_tdir.o m_ts_voltage.o mesh.o meshdscf.o meshsubs.o
+dhscf.o: moremeshsubs.o parallel.o parsing.o precision.o rhofft.o rhoofd.o
+dhscf.o: siesta_options.o sys_local.o units.o vmat.o
 diag2g.o: fermid.o parallel.o parallelsubs.o precision.o sys_local.o
 diag2k.o: fermid.o parallel.o parallelsubs.o precision.o sys_local.o
 diag3g.o: fermid.o parallel.o parallelsubs.o precision.o sys_local.o
@@ -211,7 +213,7 @@
 domain_decom.o: sparse_matrices.o sys_local.o
 doping_uniform.o: alloc.o m_ntm.o mesh.o parallel.o precision.o sys_local.o
 dynamics.o: alloc.o atomlist.o files.o ioxv.o m_mpi_utils.o parallel.o
-dynamics.o: precision.o sys_local.o units.o
+dynamics.o: precision.o siesta_options.o sys_local.o units.o
 egandd.o: alloc.o denmat.o ener3.o globalise.o gradient.o on_subs.o onmod.o
 egandd.o: onmod.o precision.o sys_local.o
 eggbox.o: parallel.o precision.o
@@ -220,6 +222,11 @@
 elpa2.o: elpa1.o
 ener3.o: alloc.o globalise.o m_mpi_utils.o onmod.o precision.o
 ener3lomem.o: alloc.o globalise.o m_mpi_utils.o onmod.o precision.o
+evolg.o: alloc.o m_eo.o m_matswinvers.o m_steps.o parallel.o parallelsubs.o
+evolg.o: precision.o siesta_options.o sparse_matrices.o sys_local.o
+evolg.o: wavefunctions.o
+evolk.o: m_diagon_opt.o m_matswinvers.o parallel.o precision.o siesta_options.o
+evolk.o: sys_local.o wavefunctions.o
 extrapolateSpData2D.o: class_Data2D.o class_OrbitalDistribution.o
 extrapolateSpData2D.o: class_SpData2D.o class_Sparsity.o
 extrapolon.o: parallel.o precision.o sys_local.o
@@ -260,9 +267,9 @@
 initatom.o: atmparams.o atom.o atom_options.o basis_io.o basis_specs.o
 initatom.o: basis_types.o chemical.o electrostatic.o ldau_specs.o m_spin.o
 initatom.o: old_atmfuncs.o precision.o pseudopotential.o
-initparallel.o: alloc.o class_OrbitalDistribution.o domain_decom.o parallel.o
-initparallel.o: parallelsubs.o precision.o siesta_options.o sparse_matrices.o
-initparallel.o: spatial.o sys_local.o
+initparallel.o: alloc.o class_OrbitalDistribution.o domain_decom.o
+initparallel.o: m_poisson_fdf.o parallel.o parallelsubs.o precision.o
+initparallel.o: siesta_options.o sparse_matrices.o spatial.o sys_local.o
 io.o: sys_local.o
 iocg.o: files.o precision.o sys_local.o
 iodm.o: alloc.o files.o parallel.o parallelsubs.o precision.o sys_local.o
@@ -325,7 +332,10 @@
 m_efield.o: siesta_geom.o sys_local.o units.o
 m_energies.o: precision.o
 m_eo.o: precision.o
+m_evolve.o: alloc.o m_memory.o parallel.o parallelsubs.o precision.o
+m_evolve.o: sys_local.o
 m_exp_coord.o: files.o m_os.o parallel.o precision.o units.o
+m_fastconv_dims.o: precision.o
 m_filter.o: bessph.o precision.o radfft.o sorting.o sys_local.o
 m_fire.o: parallel.o precision.o
 m_fire_mixing.o: alloc.o m_fire_para.o m_mpi_utils.o parallel.o precision.o
@@ -347,10 +357,14 @@
 m_hs_matrix.o: alloc.o cellsubs.o geom_helper.o precision.o sys_local.o
 m_hsx.o: atm_types.o atmfuncs.o atomlist.o files.o parallel.o parallelsubs.o
 m_hsx.o: precision.o siesta_geom.o sys_local.o
+m_initwf.o: alloc.o densematrix.o fermid.o m_diagon_opt.o m_eo.o m_memory.o
+m_initwf.o: parallel.o parallelsubs.o precision.o sparse_matrices.o sys_local.o
+m_initwf.o: wavefunctions.o
 m_integrate.o: precision.o
 m_io.o: sys_local.o
 m_io_s.o: class_OrbitalDistribution.o class_SpData1D.o class_SpData2D.o
 m_io_s.o: class_Sparsity.o parallel.o precision.o
+m_iodipol.o: files.o m_dipol.o m_steps.o
 m_iodm.o: class_OrbitalDistribution.o class_SpData2D.o class_Sparsity.o
 m_iodm.o: m_io_s.o parallel.o
 m_iodm_old.o: alloc.o files.o m_matio.o m_mpi_utils.o parallel.o parallelsubs.o
@@ -358,9 +372,13 @@
 m_iorho.o: alloc.o parallel.o parallelsubs.o precision.o sys_local.o
 m_iostruct.o: alloc.o files.o m_mpi_utils.o parallel.o precision.o
 m_iostruct.o: siesta_geom.o sys_local.o units.o
+m_iotddft.o: files.o m_dipol.o m_steps.o parallel.o siesta_options.o
+m_iotddft.o: wavefunctions.o
 m_kinetic.o: precision.o
 m_mat_invert.o: intrinsic_missing.o m_pivot_array.o precision.o
+m_matdiag.o: precision.o
 m_matio.o: alloc.o parallel.o
+m_matswinvers.o: precision.o
 m_memory.o: memoryinfo.o parallel.o sys_local.o
 m_mesh_node.o: intrinsic_missing.o parallel.o precision.o
 m_mixing.o: class_Data1D.o class_Fstack_Data1D.o parallel.o precision.o
@@ -409,6 +427,15 @@
 m_planewavematrix.o: atomlist.o dhscf.o m_planewavematrixvar.o precision.o
 m_planewavematrix.o: siesta_geom.o sparse_matrices.o
 m_planewavematrixvar.o: precision.o
+m_plot_1d_arr3d.o: precision.o
+m_poisson_fdf.o: m_char.o precision.o
+m_poisson_obc_comm.o: atm_types.o atmfuncs.o m_mpi_stop.o parallel.o
+m_poisson_obc_comm.o: precision.o siesta_geom.o sys_local.o
+m_poisson_obc_fastconv.o: cellsubs.o m_mpi_stop.o m_poisson_obc_fastconv_r2c.o
+m_poisson_obc_fastconv.o: parallel.o precision.o sys_local.o
+m_poisson_obc_fastconv_r2c.o: cellsubs.o m_fastconv_dims.o m_plot_1d_arr3d.o
+m_poisson_obc_fastconv_r2c.o: m_poisson_obc_comm.o precision.o sys_local.o
+m_poisson_switch.o: m_poisson_obc_fastconv.o precision.o sys_local.o
 m_pulay.o: alloc.o atomlist.o m_mpi_utils.o m_spin.o m_svd.o parallel.o
 m_pulay.o: precision.o siesta_options.o sparse_matrices.o sys_local.o
 m_redist_spmatrix.o: alloc.o class_Distribution.o
@@ -424,7 +451,8 @@
 m_sparsity_handling.o: class_Data1D.o class_Data2D.o
 m_sparsity_handling.o: class_OrbitalDistribution.o class_SpData1D.o
 m_sparsity_handling.o: class_SpData2D.o class_Sparsity.o geom_helper.o
-m_sparsity_handling.o: m_interpolate.o m_region.o precision.o
+m_sparsity_handling.o: intrinsic_missing.o m_interpolate.o m_region.o
+m_sparsity_handling.o: precision.o
 m_spin.o: alloc.o parallel.o precision.o sys_local.o units.o
 m_stress.o: precision.o
 m_supercell.o: atom_graph.o class_OrbitalDistribution.o class_SpData2D.o
@@ -538,9 +566,10 @@
 m_ts_options.o: m_ts_hartree.o m_ts_method.o m_ts_mumps_init.o m_ts_tdir.o
 m_ts_options.o: m_ts_weight.o parallel.o precision.o siesta_options.o units.o
 m_ts_pivot.o: alloc.o class_OrbitalDistribution.o class_Sparsity.o
-m_ts_pivot.o: create_Sparsity_SC.o create_Sparsity_Union.o intrinsic_missing.o
-m_ts_pivot.o: m_char.o m_pivot.o m_region.o m_sparsity_handling.o m_ts_debug.o
-m_ts_pivot.o: m_ts_electype.o parallel.o precision.o
+m_ts_pivot.o: create_Sparsity_SC.o create_Sparsity_Union.o geom_helper.o
+m_ts_pivot.o: intrinsic_missing.o m_char.o m_pivot.o m_region.o
+m_ts_pivot.o: m_sparsity_handling.o m_ts_debug.o m_ts_electype.o parallel.o
+m_ts_pivot.o: precision.o
 m_ts_rgn2trimat.o: alloc.o class_OrbitalDistribution.o class_Sparsity.o
 m_ts_rgn2trimat.o: create_Sparsity_Union.o geom_helper.o m_region.o
 m_ts_rgn2trimat.o: m_ts_electype.o m_ts_method.o m_ts_tri_common.o parallel.o
@@ -685,8 +714,8 @@
 rdiag.o: alloc.o m_diagon_opt.o memoryinfo.o parallel.o precision.o sys_local.o
 rdiag_elpa.o: alloc.o elpa1.o elpa2.o parallel.o precision.o sys_local.o
 rdiag_mrrr.o: alloc.o parallel.o precision.o sys_local.o
-read_options.o: files.o m_charge_add.o m_hartree_add.o m_mixing_scf.o m_spin.o
-read_options.o: m_target_stress.o memoryinfo.o parallel.o precision.o
+read_options.o: files.o m_charge_add.o m_cite.o m_hartree_add.o m_mixing_scf.o
+read_options.o: m_spin.o m_target_stress.o memoryinfo.o parallel.o precision.o
 read_options.o: siesta_cml.o siesta_options.o sys_local.o units.o
 read_xc_info.o: parallel.o precision.o sys_local.o
 readsp.o: parallel.o precision.o sys_local.o
@@ -756,30 +785,32 @@
 siesta_end.o: meshdscf.o meshphi.o moremeshsubs.o parallel.o siesta_cml.o
 siesta_end.o: siesta_cmlsubs.o siesta_master.o siesta_options.o
 siesta_end.o: sparse_matrices.o timestamp.o
-siesta_forces.o: atomlist.o class_Fstack_Data1D.o class_SpData2D.o compute_dm.o
-siesta_forces.o: compute_energies.o compute_max_diff.o files.o
-siesta_forces.o: final_H_f_stress.o flook_siesta.o kpoint_grid.o
-siesta_forces.o: m_convergence.o m_energies.o m_forces.o m_iodm_old.o
-siesta_forces.o: m_mixing.o m_mixing_scf.o m_mpi_utils.o m_ncdf_siesta.o
-siesta_forces.o: m_pexsi.o m_pexsi_driver.o m_rhog.o m_spin.o m_steps.o
-siesta_forces.o: m_stress.o m_transiesta.o m_ts_charge.o m_ts_electype.o
-siesta_forces.o: m_ts_global_vars.o m_ts_method.o m_ts_options.o mixer.o
-siesta_forces.o: parallel.o post_scf_work.o precision.o save_density_matrix.o
-siesta_forces.o: scfconvergence_test.o setup_H0.o setup_hamiltonian.o
-siesta_forces.o: siesta_cml.o siesta_dicts.o siesta_geom.o siesta_master.o
-siesta_forces.o: siesta_options.o sparse_matrices.o state_analysis.o
-siesta_forces.o: state_init.o sys_local.o units.o write_subs.o
+siesta_forces.o: alloc.o atomlist.o class_Fstack_Data1D.o class_SpData2D.o
+siesta_forces.o: compute_dm.o compute_energies.o compute_max_diff.o
+siesta_forces.o: densematrix.o files.o final_H_f_stress.o flook_siesta.o
+siesta_forces.o: kpoint_grid.o m_convergence.o m_energies.o m_eo.o m_evolve.o
+siesta_forces.o: m_forces.o m_gamma.o m_initwf.o m_iodm.o m_iodm_old.o
+siesta_forces.o: m_iotddft.o m_mixing.o m_mixing_scf.o m_mpi_utils.o
+siesta_forces.o: m_ncdf_siesta.o m_pexsi.o m_pexsi_driver.o m_rhog.o m_spin.o
+siesta_forces.o: m_steps.o m_stress.o m_transiesta.o m_ts_charge.o
+siesta_forces.o: m_ts_electype.o m_ts_global_vars.o m_ts_method.o
+siesta_forces.o: m_ts_options.o mixer.o overfsm.o parallel.o post_scf_work.o
+siesta_forces.o: precision.o save_density_matrix.o scfconvergence_test.o
+siesta_forces.o: setup_H0.o setup_hamiltonian.o siesta_cml.o siesta_dicts.o
+siesta_forces.o: siesta_geom.o siesta_master.o siesta_options.o
+siesta_forces.o: sparse_matrices.o state_analysis.o state_init.o sys_local.o
+siesta_forces.o: units.o wavefunctions.o write_subs.o
 siesta_geom.o: precision.o
 siesta_init.o: alloc.o atomlist.o bands.o bsc_xcmod.o
 siesta_init.o: class_Fstack_Pair_Geometry_SpData2D.o densematrix.o files.o
 siesta_init.o: flook_siesta.o ioxv.o kpoint_grid.o kpoint_pdos.o ksvinit.o
 siesta_init.o: m_cite.o m_energies.o m_eo.o m_fixed.o m_forces.o m_gamma.o
-siesta_init.o: m_iostruct.o m_mpi_utils.o m_new_dm.o m_rmaxh.o m_spin.o
-siesta_init.o: m_steps.o m_supercell.o m_timer.o m_wallclock.o metaforce.o
-siesta_init.o: molecularmechanics.o object_debug.o parallel.o parallelsubs.o
-siesta_init.o: projected_DOS.o siesta_cmlsubs.o siesta_dicts.o siesta_geom.o
-siesta_init.o: siesta_options.o sparse_matrices.o struct_init.o sys_local.o
-siesta_init.o: timer.o timestamp.o ts_init.o writewave.o zmatrix.o
+siesta_init.o: m_iostruct.o m_mpi_utils.o m_new_dm.o m_poisson_fdf.o m_rmaxh.o
+siesta_init.o: m_spin.o m_steps.o m_supercell.o m_timer.o m_wallclock.o
+siesta_init.o: metaforce.o molecularmechanics.o object_debug.o parallel.o
+siesta_init.o: parallelsubs.o projected_DOS.o siesta_cmlsubs.o siesta_dicts.o
+siesta_init.o: siesta_geom.o siesta_options.o sparse_matrices.o struct_init.o
+siesta_init.o: sys_local.o timer.o timestamp.o ts_init.o writewave.o zmatrix.o
 siesta_master.o: iopipes.o iosockets.o precision.o sys_local.o
 siesta_move.o: atomlist.o broyden_optim.o cell_broyden_optim.o
 siesta_move.o: cell_fire_optim.o dynamics.o fire_optim.o flook_siesta.o ioxv.o
@@ -813,19 +844,18 @@
 state_init.o: proximity_check.o siesta_cml.o siesta_geom.o siesta_options.o
 state_init.o: sparse_matrices.o sys_local.o units.o write_subs.o zmatrix.o
 struct_init.o: alloc.o atmfuncs.o atomlist.o files.o ioxv.o m_exp_coord.o
-struct_init.o: m_iostruct.o m_mpi_utils.o m_steps.o parallel.o periodic_table.o
-struct_init.o: siesta_cml.o siesta_geom.o siesta_master.o siesta_options.o
-struct_init.o: units.o zmatrix.o
+struct_init.o: m_iostruct.o m_mpi_utils.o m_poisson_fdf.o m_poisson_obc_comm.o
+struct_init.o: m_steps.o parallel.o periodic_table.o siesta_cml.o siesta_geom.o
+struct_init.o: siesta_master.o siesta_options.o units.o zmatrix.o
 sys.o: parallel.o
 timer.o: extrae_eventllist.o extrae_module.o m_timer.o parallel.o sys_local.o
 timer.o: timer_tree.o
 timer_tree.o: m_walltime.o
 transition_rate.o: alloc.o fermid.o parallel.o parallelsubs.o precision.o
 transition_rate.o: sys_local.o
-ts_init.o: m_cite.o m_fixed.o m_os.o m_ts_cctype.o m_ts_charge.o
-ts_init.o: m_ts_electrode.o m_ts_electype.o m_ts_gf.o m_ts_global_vars.o
-ts_init.o: m_ts_kpoints.o m_ts_method.o m_ts_options.o parallel.o
-ts_init.o: siesta_options.o
+ts_init.o: m_fixed.o m_os.o m_ts_cctype.o m_ts_charge.o m_ts_electrode.o
+ts_init.o: m_ts_electype.o m_ts_gf.o m_ts_global_vars.o m_ts_kpoints.o
+ts_init.o: m_ts_method.o m_ts_options.o parallel.o siesta_options.o
 ts_show_regions.o: m_region.o m_ts_electype.o m_ts_method.o parallel.o
 ts_show_regions.o: precision.o units.o
 typecell.o: precision.o
@@ -837,6 +867,7 @@
 vmatsp.o: alloc.o atm_types.o atmfuncs.o atomlist.o listsc.o mesh.o meshdscf.o
 vmatsp.o: meshphi.o precision.o
 vmb.o: m_fixed.o parallel.o precision.o sys_local.o
+wavefunctions.o: parallel.o precision.o
 write_inp_wannier.o: alloc.o atmfuncs.o atomlist.o m_ntm.o m_orderbands.o
 write_inp_wannier.o: mneighb.o parallel.o parallelsubs.o precision.o
 write_inp_wannier.o: siesta2wannier90.o siesta_geom.o
@@ -968,6 +999,7 @@
 m_ksv.o: ksv.o
 m_ksvinit.o: ksvinit.o
 m_ldau.o: ldau.o
+m_linpack.o: linpack.o
 m_local_dos.o: local_DOS.o
 m_matel_registry.o: matel_registry.o
 m_minvec.o: minvec.o
@@ -1012,6 +1044,8 @@
 m_vmat.o: vmat.o
 m_zm_broyden_optim.o: zm_broyden_optim.o
 m_zm_fire_optim.o: zm_fire_optim.o
+matdiagon.o: m_matdiag.o
+matswinversion.o: m_matswinvers.o
 mtprng.o: m_uuid.o
 neighbour.o: mneighb.o
 on_core.o: onmod.o

=== modified file 'Util/Denchar/Src/Makefile'
--- Util/Denchar/Src/Makefile	2016-08-30 14:02:33 +0000
+++ Util/Denchar/Src/Makefile	2016-11-22 08:16:45 +0000
@@ -175,7 +175,7 @@
 bsc_cellxc.o: alloc.o bsc_xcmod.o cellxc_mod.o mesh.o moremeshsubs.o parallel.o
 bsc_cellxc.o: parallelsubs.o precision.o sys.o
 bsc_xcmod.o: parallel.o precision.o sys.o
-cart2frac.o: sys.o
+cart2frac.o: precision.o sys.o
 cdiag.o: alloc.o m_diagon_opt.o memoryinfo.o parallel.o precision.o sys.o
 cell_broyden_optim.o: m_broyddj_nocomm.o parallel.o precision.o sys.o units.o
 cell_broyden_optim.o: zmatrix.o
@@ -187,6 +187,8 @@
 cgvc_zmatrix.o: alloc.o conjgr.o m_mpi_utils.o parallel.o precision.o sys.o
 cgvc_zmatrix.o: units.o zmatrix.o
 cgwf.o: egandd.o onmod.o parallel.o precision.o sys.o
+changebasis.o: alloc.o m_matdiag.o parallel.o parallelsubs.o precision.o
+changebasis.o: sparse_matrices.o sys.o wavefunctions.o
 chemical.o: parallel.o precision.o sys.o
 chempot.o: alloc.o mneighb.o parallel.o parallelsubs.o precision.o sys.o
 chkdim.o: sys.o
@@ -238,17 +240,17 @@
 denmat.o: alloc.o globalise.o m_mpi_utils.o onmod.o precision.o
 denmatlomem.o: alloc.o globalise.o onmod.o precision.o
 densematrix.o: alloc.o precision.o
-detover.o: alloc.o parallel.o parallelsubs.o precision.o
+detover.o: alloc.o linpack.o parallel.o parallelsubs.o precision.o
 dfscf.o: alloc.o atm_types.o atmfuncs.o atomlist.o listsc.o m_spin.o mesh.o
 dfscf.o: meshdscf.o meshphi.o parallel.o parallelsubs.o precision.o sys.o
 dhscf.o: alloc.o atmfuncs.o bsc_xcmod.o cellxc_mod.o delk.o dfscf.o
 dhscf.o: doping_uniform.o files.o forhar.o iogrid_netcdf.o m_charge_add.o
-dhscf.o: m_efield.o m_hartree_add.o m_iorho.o m_mesh_node.o m_ncdf_io.o
-dhscf.o: m_ncdf_siesta.o m_partial_charges.o m_rhog.o m_spin.o
-dhscf.o: m_ts_global_vars.o m_ts_hartree.o m_ts_options.o m_ts_tdir.o
-dhscf.o: m_ts_voltage.o mesh.o meshdscf.o meshsubs.o moremeshsubs.o parallel.o
-dhscf.o: parsing.o precision.o rhofft.o rhoofd.o siesta_options.o sys.o units.o
-dhscf.o: vmat.o
+dhscf.o: m_efield.o m_hartree_add.o m_iorho.o m_iotddft.o m_mesh_node.o
+dhscf.o: m_ncdf_io.o m_ncdf_siesta.o m_partial_charges.o m_poisson_fdf.o
+dhscf.o: m_poisson_switch.o m_rhog.o m_spin.o m_ts_global_vars.o m_ts_hartree.o
+dhscf.o: m_ts_options.o m_ts_tdir.o m_ts_voltage.o mesh.o meshdscf.o meshsubs.o
+dhscf.o: moremeshsubs.o parallel.o parsing.o precision.o rhofft.o rhoofd.o
+dhscf.o: siesta_options.o sys.o units.o vmat.o
 diag2g.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
 diag2k.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
 diag3g.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
@@ -276,7 +278,7 @@
 domain_decom.o: sparse_matrices.o sys.o
 doping_uniform.o: alloc.o m_ntm.o mesh.o parallel.o precision.o sys.o
 dynamics.o: alloc.o atomlist.o files.o ioxv.o m_mpi_utils.o parallel.o
-dynamics.o: precision.o sys.o units.o
+dynamics.o: precision.o siesta_options.o sys.o units.o
 egandd.o: alloc.o denmat.o ener3.o globalise.o gradient.o on_subs.o onmod.o
 egandd.o: onmod.o precision.o sys.o
 eggbox.o: parallel.o precision.o
@@ -285,6 +287,10 @@
 elpa2.o: elpa1.o
 ener3.o: alloc.o globalise.o m_mpi_utils.o onmod.o precision.o
 ener3lomem.o: alloc.o globalise.o m_mpi_utils.o onmod.o precision.o
+evolg.o: alloc.o m_eo.o m_matswinvers.o m_steps.o parallel.o parallelsubs.o
+evolg.o: precision.o siesta_options.o sparse_matrices.o sys.o wavefunctions.o
+evolk.o: m_diagon_opt.o m_matswinvers.o parallel.o precision.o siesta_options.o
+evolk.o: sys.o wavefunctions.o
 extrapolateSpData2D.o: class_Data2D.o class_OrbitalDistribution.o
 extrapolateSpData2D.o: class_SpData2D.o class_Sparsity.o
 extrapolon.o: parallel.o precision.o sys.o
@@ -325,9 +331,9 @@
 initatom.o: atmparams.o atom.o atom_options.o basis_io.o basis_specs.o
 initatom.o: basis_types.o chemical.o electrostatic.o ldau_specs.o m_spin.o
 initatom.o: old_atmfuncs.o precision.o pseudopotential.o
-initparallel.o: alloc.o class_OrbitalDistribution.o domain_decom.o parallel.o
-initparallel.o: parallelsubs.o precision.o siesta_options.o sparse_matrices.o
-initparallel.o: spatial.o sys.o
+initparallel.o: alloc.o class_OrbitalDistribution.o domain_decom.o
+initparallel.o: m_poisson_fdf.o parallel.o parallelsubs.o precision.o
+initparallel.o: siesta_options.o sparse_matrices.o spatial.o sys.o
 io.o: sys.o
 iocg.o: files.o precision.o sys.o
 iodm.o: alloc.o files.o parallel.o parallelsubs.o precision.o sys.o
@@ -386,7 +392,9 @@
 m_efield.o: siesta_geom.o sys.o units.o
 m_energies.o: precision.o
 m_eo.o: precision.o
+m_evolve.o: alloc.o m_memory.o parallel.o parallelsubs.o precision.o sys.o
 m_exp_coord.o: files.o m_os.o parallel.o precision.o units.o
+m_fastconv_dims.o: precision.o
 m_filter.o: bessph.o precision.o radfft.o sorting.o sys.o
 m_fire.o: parallel.o precision.o
 m_fire_mixing.o: alloc.o m_fire_para.o m_mpi_utils.o parallel.o precision.o
@@ -408,10 +416,14 @@
 m_hs_matrix.o: alloc.o cellsubs.o geom_helper.o precision.o sys.o
 m_hsx.o: atm_types.o atmfuncs.o atomlist.o files.o parallel.o parallelsubs.o
 m_hsx.o: precision.o siesta_geom.o sys.o
+m_initwf.o: alloc.o densematrix.o fermid.o m_diagon_opt.o m_eo.o m_memory.o
+m_initwf.o: parallel.o parallelsubs.o precision.o sparse_matrices.o sys.o
+m_initwf.o: wavefunctions.o
 m_integrate.o: precision.o
 m_io.o: sys.o
 m_io_s.o: class_OrbitalDistribution.o class_SpData1D.o class_SpData2D.o
 m_io_s.o: class_Sparsity.o parallel.o precision.o
+m_iodipol.o: files.o m_dipol.o m_steps.o
 m_iodm.o: class_OrbitalDistribution.o class_SpData2D.o class_Sparsity.o
 m_iodm.o: m_io_s.o parallel.o
 m_iodm_old.o: alloc.o files.o m_matio.o m_mpi_utils.o parallel.o parallelsubs.o
@@ -419,9 +431,13 @@
 m_iorho.o: alloc.o parallel.o parallelsubs.o precision.o sys.o
 m_iostruct.o: alloc.o files.o m_mpi_utils.o parallel.o precision.o
 m_iostruct.o: siesta_geom.o sys.o units.o
+m_iotddft.o: files.o m_dipol.o m_steps.o parallel.o siesta_options.o
+m_iotddft.o: wavefunctions.o
 m_kinetic.o: precision.o
 m_mat_invert.o: intrinsic_missing.o m_pivot_array.o precision.o
+m_matdiag.o: precision.o
 m_matio.o: alloc.o parallel.o
+m_matswinvers.o: precision.o
 m_memory.o: memoryinfo.o parallel.o sys.o
 m_mesh_node.o: intrinsic_missing.o parallel.o precision.o
 m_mixing.o: class_Data1D.o class_Fstack_Data1D.o parallel.o precision.o
@@ -469,6 +485,15 @@
 m_planewavematrix.o: atomlist.o dhscf.o m_planewavematrixvar.o precision.o
 m_planewavematrix.o: siesta_geom.o sparse_matrices.o
 m_planewavematrixvar.o: precision.o
+m_plot_1d_arr3d.o: precision.o
+m_poisson_fdf.o: m_char.o precision.o
+m_poisson_obc_comm.o: atm_types.o atmfuncs.o m_mpi_stop.o parallel.o
+m_poisson_obc_comm.o: precision.o siesta_geom.o sys.o
+m_poisson_obc_fastconv.o: cellsubs.o m_mpi_stop.o m_poisson_obc_fastconv_r2c.o
+m_poisson_obc_fastconv.o: parallel.o precision.o sys.o
+m_poisson_obc_fastconv_r2c.o: cellsubs.o m_fastconv_dims.o m_plot_1d_arr3d.o
+m_poisson_obc_fastconv_r2c.o: m_poisson_obc_comm.o precision.o sys.o
+m_poisson_switch.o: m_poisson_obc_fastconv.o precision.o sys.o
 m_pulay.o: alloc.o atomlist.o m_mpi_utils.o m_spin.o m_svd.o parallel.o
 m_pulay.o: precision.o siesta_options.o sparse_matrices.o sys.o
 m_redist_spmatrix.o: alloc.o class_Distribution.o
@@ -484,7 +509,8 @@
 m_sparsity_handling.o: class_Data1D.o class_Data2D.o
 m_sparsity_handling.o: class_OrbitalDistribution.o class_SpData1D.o
 m_sparsity_handling.o: class_SpData2D.o class_Sparsity.o geom_helper.o
-m_sparsity_handling.o: m_interpolate.o m_region.o precision.o
+m_sparsity_handling.o: intrinsic_missing.o m_interpolate.o m_region.o
+m_sparsity_handling.o: precision.o
 m_spin.o: alloc.o parallel.o precision.o sys.o units.o
 m_stress.o: precision.o
 m_supercell.o: atom_graph.o class_OrbitalDistribution.o class_SpData2D.o
@@ -596,9 +622,10 @@
 m_ts_options.o: m_ts_hartree.o m_ts_method.o m_ts_mumps_init.o m_ts_tdir.o
 m_ts_options.o: m_ts_weight.o parallel.o precision.o siesta_options.o units.o
 m_ts_pivot.o: alloc.o class_OrbitalDistribution.o class_Sparsity.o
-m_ts_pivot.o: create_Sparsity_SC.o create_Sparsity_Union.o intrinsic_missing.o
-m_ts_pivot.o: m_char.o m_pivot.o m_region.o m_sparsity_handling.o m_ts_debug.o
-m_ts_pivot.o: m_ts_electype.o parallel.o precision.o
+m_ts_pivot.o: create_Sparsity_SC.o create_Sparsity_Union.o geom_helper.o
+m_ts_pivot.o: intrinsic_missing.o m_char.o m_pivot.o m_region.o
+m_ts_pivot.o: m_sparsity_handling.o m_ts_debug.o m_ts_electype.o parallel.o
+m_ts_pivot.o: precision.o
 m_ts_rgn2trimat.o: alloc.o class_OrbitalDistribution.o class_Sparsity.o
 m_ts_rgn2trimat.o: create_Sparsity_Union.o geom_helper.o m_region.o
 m_ts_rgn2trimat.o: m_ts_electype.o m_ts_method.o m_ts_tri_common.o parallel.o
@@ -741,8 +768,8 @@
 rdiag.o: alloc.o m_diagon_opt.o memoryinfo.o parallel.o precision.o sys.o
 rdiag_elpa.o: alloc.o elpa1.o elpa2.o parallel.o precision.o sys.o
 rdiag_mrrr.o: alloc.o parallel.o precision.o sys.o
-read_options.o: files.o m_charge_add.o m_hartree_add.o m_mixing_scf.o m_spin.o
-read_options.o: m_target_stress.o memoryinfo.o parallel.o precision.o
+read_options.o: files.o m_charge_add.o m_cite.o m_hartree_add.o m_mixing_scf.o
+read_options.o: m_spin.o m_target_stress.o memoryinfo.o parallel.o precision.o
 read_options.o: siesta_cml.o siesta_options.o sys.o units.o
 read_xc_info.o: parallel.o precision.o sys.o
 readsp.o: parallel.o precision.o sys.o
@@ -809,30 +836,32 @@
 siesta_end.o: meshdscf.o meshphi.o moremeshsubs.o parallel.o siesta_cml.o
 siesta_end.o: siesta_cmlsubs.o siesta_master.o siesta_options.o
 siesta_end.o: sparse_matrices.o timestamp.o
-siesta_forces.o: atomlist.o class_Fstack_Data1D.o class_SpData2D.o compute_dm.o
-siesta_forces.o: compute_energies.o compute_max_diff.o files.o
-siesta_forces.o: final_H_f_stress.o flook_siesta.o kpoint_grid.o
-siesta_forces.o: m_convergence.o m_energies.o m_forces.o m_iodm_old.o
-siesta_forces.o: m_mixing.o m_mixing_scf.o m_mpi_utils.o m_ncdf_siesta.o
-siesta_forces.o: m_pexsi.o m_pexsi_driver.o m_rhog.o m_spin.o m_steps.o
-siesta_forces.o: m_stress.o m_transiesta.o m_ts_charge.o m_ts_electype.o
-siesta_forces.o: m_ts_global_vars.o m_ts_method.o m_ts_options.o mixer.o
-siesta_forces.o: parallel.o post_scf_work.o precision.o save_density_matrix.o
-siesta_forces.o: scfconvergence_test.o setup_H0.o setup_hamiltonian.o
-siesta_forces.o: siesta_cml.o siesta_dicts.o siesta_geom.o siesta_master.o
-siesta_forces.o: siesta_options.o sparse_matrices.o state_analysis.o
-siesta_forces.o: state_init.o sys.o units.o write_subs.o
+siesta_forces.o: alloc.o atomlist.o class_Fstack_Data1D.o class_SpData2D.o
+siesta_forces.o: compute_dm.o compute_energies.o compute_max_diff.o
+siesta_forces.o: densematrix.o files.o final_H_f_stress.o flook_siesta.o
+siesta_forces.o: kpoint_grid.o m_convergence.o m_energies.o m_eo.o m_evolve.o
+siesta_forces.o: m_forces.o m_gamma.o m_initwf.o m_iodm.o m_iodm_old.o
+siesta_forces.o: m_iotddft.o m_mixing.o m_mixing_scf.o m_mpi_utils.o
+siesta_forces.o: m_ncdf_siesta.o m_pexsi.o m_pexsi_driver.o m_rhog.o m_spin.o
+siesta_forces.o: m_steps.o m_stress.o m_transiesta.o m_ts_charge.o
+siesta_forces.o: m_ts_electype.o m_ts_global_vars.o m_ts_method.o
+siesta_forces.o: m_ts_options.o mixer.o overfsm.o parallel.o post_scf_work.o
+siesta_forces.o: precision.o save_density_matrix.o scfconvergence_test.o
+siesta_forces.o: setup_H0.o setup_hamiltonian.o siesta_cml.o siesta_dicts.o
+siesta_forces.o: siesta_geom.o siesta_master.o siesta_options.o
+siesta_forces.o: sparse_matrices.o state_analysis.o state_init.o sys.o units.o
+siesta_forces.o: wavefunctions.o write_subs.o
 siesta_geom.o: precision.o
 siesta_init.o: alloc.o atomlist.o bands.o bsc_xcmod.o
 siesta_init.o: class_Fstack_Pair_Geometry_SpData2D.o densematrix.o files.o
 siesta_init.o: flook_siesta.o ioxv.o kpoint_grid.o kpoint_pdos.o ksvinit.o
 siesta_init.o: m_cite.o m_energies.o m_eo.o m_fixed.o m_forces.o m_gamma.o
-siesta_init.o: m_iostruct.o m_mpi_utils.o m_new_dm.o m_rmaxh.o m_spin.o
-siesta_init.o: m_steps.o m_supercell.o m_timer.o m_wallclock.o metaforce.o
-siesta_init.o: molecularmechanics.o object_debug.o parallel.o parallelsubs.o
-siesta_init.o: projected_DOS.o siesta_cmlsubs.o siesta_dicts.o siesta_geom.o
-siesta_init.o: siesta_options.o sparse_matrices.o struct_init.o sys.o timer.o
-siesta_init.o: timestamp.o ts_init.o writewave.o zmatrix.o
+siesta_init.o: m_iostruct.o m_mpi_utils.o m_new_dm.o m_poisson_fdf.o m_rmaxh.o
+siesta_init.o: m_spin.o m_steps.o m_supercell.o m_timer.o m_wallclock.o
+siesta_init.o: metaforce.o molecularmechanics.o object_debug.o parallel.o
+siesta_init.o: parallelsubs.o projected_DOS.o siesta_cmlsubs.o siesta_dicts.o
+siesta_init.o: siesta_geom.o siesta_options.o sparse_matrices.o struct_init.o
+siesta_init.o: sys.o timer.o timestamp.o ts_init.o writewave.o zmatrix.o
 siesta_master.o: iopipes.o iosockets.o precision.o sys.o
 siesta_move.o: atomlist.o broyden_optim.o cell_broyden_optim.o
 siesta_move.o: cell_fire_optim.o dynamics.o fire_optim.o flook_siesta.o ioxv.o
@@ -864,18 +893,17 @@
 state_init.o: proximity_check.o siesta_cml.o siesta_geom.o siesta_options.o
 state_init.o: sparse_matrices.o sys.o units.o write_subs.o zmatrix.o
 struct_init.o: alloc.o atmfuncs.o atomlist.o files.o ioxv.o m_exp_coord.o
-struct_init.o: m_iostruct.o m_mpi_utils.o m_steps.o parallel.o periodic_table.o
-struct_init.o: siesta_cml.o siesta_geom.o siesta_master.o siesta_options.o
-struct_init.o: units.o zmatrix.o
+struct_init.o: m_iostruct.o m_mpi_utils.o m_poisson_fdf.o m_poisson_obc_comm.o
+struct_init.o: m_steps.o parallel.o periodic_table.o siesta_cml.o siesta_geom.o
+struct_init.o: siesta_master.o siesta_options.o units.o zmatrix.o
 sys.o: parallel.o
 timer.o: extrae_eventllist.o extrae_module.o m_timer.o parallel.o sys.o
 timer.o: timer_tree.o
 timer_tree.o: m_walltime.o
 transition_rate.o: alloc.o fermid.o parallel.o parallelsubs.o precision.o sys.o
-ts_init.o: m_cite.o m_fixed.o m_os.o m_ts_cctype.o m_ts_charge.o
-ts_init.o: m_ts_electrode.o m_ts_electype.o m_ts_gf.o m_ts_global_vars.o
-ts_init.o: m_ts_kpoints.o m_ts_method.o m_ts_options.o parallel.o
-ts_init.o: siesta_options.o
+ts_init.o: m_fixed.o m_os.o m_ts_cctype.o m_ts_charge.o m_ts_electrode.o
+ts_init.o: m_ts_electype.o m_ts_gf.o m_ts_global_vars.o m_ts_kpoints.o
+ts_init.o: m_ts_method.o m_ts_options.o parallel.o siesta_options.o
 ts_show_regions.o: m_region.o m_ts_electype.o m_ts_method.o parallel.o
 ts_show_regions.o: precision.o units.o
 typecell.o: precision.o
@@ -887,6 +915,7 @@
 vmatsp.o: alloc.o atm_types.o atmfuncs.o atomlist.o listsc.o mesh.o meshdscf.o
 vmatsp.o: meshphi.o precision.o
 vmb.o: m_fixed.o parallel.o precision.o sys.o
+wavefunctions.o: parallel.o precision.o
 write_inp_wannier.o: alloc.o atmfuncs.o atomlist.o m_ntm.o m_orderbands.o
 write_inp_wannier.o: mneighb.o parallel.o parallelsubs.o precision.o
 write_inp_wannier.o: siesta2wannier90.o siesta_geom.o
@@ -911,7 +940,9 @@
 atompla.o: precision.o sys.o
 colinear.o: precision.o
 denchar.o: basis_io.o listsc.o parallel.o precision.o
+length.o: precision.o
 local_reinit.o: files.o
+matvect.o: precision.o
 neighb.o: precision.o
 planed.o: precision.o
 ranger.o: precision.o
@@ -1020,6 +1051,7 @@
 m_ksv.o: ksv.o
 m_ksvinit.o: ksvinit.o
 m_ldau.o: ldau.o
+m_linpack.o: linpack.o
 m_local_dos.o: local_DOS.o
 m_matel_registry.o: matel_registry.o
 m_minvec.o: minvec.o
@@ -1064,6 +1096,8 @@
 m_vmat.o: vmat.o
 m_zm_broyden_optim.o: zm_broyden_optim.o
 m_zm_fire_optim.o: zm_fire_optim.o
+matdiagon.o: m_matdiag.o
+matswinversion.o: m_matswinvers.o
 mtprng.o: m_uuid.o
 neighbour.o: mneighb.o
 on_core.o: onmod.o

=== modified file 'Util/Gen-basis/Makefile'
--- Util/Gen-basis/Makefile	2016-08-30 14:02:33 +0000
+++ Util/Gen-basis/Makefile	2016-11-22 08:16:45 +0000
@@ -171,7 +171,7 @@
 bsc_cellxc.o: alloc.o bsc_xcmod.o cellxc_mod.o mesh.o moremeshsubs.o parallel.o
 bsc_cellxc.o: parallelsubs.o precision.o sys.o
 bsc_xcmod.o: parallel.o precision.o sys.o
-cart2frac.o: sys.o
+cart2frac.o: precision.o sys.o
 cdiag.o: alloc.o m_diagon_opt.o memoryinfo.o parallel.o precision.o sys.o
 cell_broyden_optim.o: m_broyddj_nocomm.o parallel.o precision.o sys.o units.o
 cell_broyden_optim.o: zmatrix.o
@@ -183,6 +183,8 @@
 cgvc_zmatrix.o: alloc.o conjgr.o m_mpi_utils.o parallel.o precision.o sys.o
 cgvc_zmatrix.o: units.o zmatrix.o
 cgwf.o: egandd.o onmod.o parallel.o precision.o sys.o
+changebasis.o: alloc.o m_matdiag.o parallel.o parallelsubs.o precision.o
+changebasis.o: sparse_matrices.o sys.o wavefunctions.o
 chemical.o: parallel.o precision.o sys.o
 chempot.o: alloc.o mneighb.o parallel.o parallelsubs.o precision.o sys.o
 chkdim.o: sys.o
@@ -234,17 +236,17 @@
 denmat.o: alloc.o globalise.o m_mpi_utils.o onmod.o precision.o
 denmatlomem.o: alloc.o globalise.o onmod.o precision.o
 densematrix.o: alloc.o precision.o
-detover.o: alloc.o parallel.o parallelsubs.o precision.o
+detover.o: alloc.o linpack.o parallel.o parallelsubs.o precision.o
 dfscf.o: alloc.o atm_types.o atmfuncs.o atomlist.o listsc.o m_spin.o mesh.o
 dfscf.o: meshdscf.o meshphi.o parallel.o parallelsubs.o precision.o sys.o
 dhscf.o: alloc.o atmfuncs.o bsc_xcmod.o cellxc_mod.o delk.o dfscf.o
 dhscf.o: doping_uniform.o files.o forhar.o iogrid_netcdf.o m_charge_add.o
-dhscf.o: m_efield.o m_hartree_add.o m_iorho.o m_mesh_node.o m_ncdf_io.o
-dhscf.o: m_ncdf_siesta.o m_partial_charges.o m_rhog.o m_spin.o
-dhscf.o: m_ts_global_vars.o m_ts_hartree.o m_ts_options.o m_ts_tdir.o
-dhscf.o: m_ts_voltage.o mesh.o meshdscf.o meshsubs.o moremeshsubs.o parallel.o
-dhscf.o: parsing.o precision.o rhofft.o rhoofd.o siesta_options.o sys.o units.o
-dhscf.o: vmat.o
+dhscf.o: m_efield.o m_hartree_add.o m_iorho.o m_iotddft.o m_mesh_node.o
+dhscf.o: m_ncdf_io.o m_ncdf_siesta.o m_partial_charges.o m_poisson_fdf.o
+dhscf.o: m_poisson_switch.o m_rhog.o m_spin.o m_ts_global_vars.o m_ts_hartree.o
+dhscf.o: m_ts_options.o m_ts_tdir.o m_ts_voltage.o mesh.o meshdscf.o meshsubs.o
+dhscf.o: moremeshsubs.o parallel.o parsing.o precision.o rhofft.o rhoofd.o
+dhscf.o: siesta_options.o sys.o units.o vmat.o
 diag2g.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
 diag2k.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
 diag3g.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
@@ -272,7 +274,7 @@
 domain_decom.o: sparse_matrices.o sys.o
 doping_uniform.o: alloc.o m_ntm.o mesh.o parallel.o precision.o sys.o
 dynamics.o: alloc.o atomlist.o files.o ioxv.o m_mpi_utils.o parallel.o
-dynamics.o: precision.o sys.o units.o
+dynamics.o: precision.o siesta_options.o sys.o units.o
 egandd.o: alloc.o denmat.o ener3.o globalise.o gradient.o on_subs.o onmod.o
 egandd.o: onmod.o precision.o sys.o
 eggbox.o: parallel.o precision.o
@@ -281,6 +283,10 @@
 elpa2.o: elpa1.o
 ener3.o: alloc.o globalise.o m_mpi_utils.o onmod.o precision.o
 ener3lomem.o: alloc.o globalise.o m_mpi_utils.o onmod.o precision.o
+evolg.o: alloc.o m_eo.o m_matswinvers.o m_steps.o parallel.o parallelsubs.o
+evolg.o: precision.o siesta_options.o sparse_matrices.o sys.o wavefunctions.o
+evolk.o: m_diagon_opt.o m_matswinvers.o parallel.o precision.o siesta_options.o
+evolk.o: sys.o wavefunctions.o
 extrapolateSpData2D.o: class_Data2D.o class_OrbitalDistribution.o
 extrapolateSpData2D.o: class_SpData2D.o class_Sparsity.o
 extrapolon.o: parallel.o precision.o sys.o
@@ -321,9 +327,9 @@
 initatom.o: atmparams.o atom.o atom_options.o basis_io.o basis_specs.o
 initatom.o: basis_types.o chemical.o electrostatic.o ldau_specs.o m_spin.o
 initatom.o: old_atmfuncs.o precision.o pseudopotential.o
-initparallel.o: alloc.o class_OrbitalDistribution.o domain_decom.o parallel.o
-initparallel.o: parallelsubs.o precision.o siesta_options.o sparse_matrices.o
-initparallel.o: spatial.o sys.o
+initparallel.o: alloc.o class_OrbitalDistribution.o domain_decom.o
+initparallel.o: m_poisson_fdf.o parallel.o parallelsubs.o precision.o
+initparallel.o: siesta_options.o sparse_matrices.o spatial.o sys.o
 io.o: sys.o
 iocg.o: files.o precision.o sys.o
 iodm.o: alloc.o files.o parallel.o parallelsubs.o precision.o sys.o
@@ -382,7 +388,9 @@
 m_efield.o: siesta_geom.o sys.o units.o
 m_energies.o: precision.o
 m_eo.o: precision.o
+m_evolve.o: alloc.o m_memory.o parallel.o parallelsubs.o precision.o sys.o
 m_exp_coord.o: files.o m_os.o parallel.o precision.o units.o
+m_fastconv_dims.o: precision.o
 m_filter.o: bessph.o precision.o radfft.o sorting.o sys.o
 m_fire.o: parallel.o precision.o
 m_fire_mixing.o: alloc.o m_fire_para.o m_mpi_utils.o parallel.o precision.o
@@ -404,10 +412,14 @@
 m_hs_matrix.o: alloc.o cellsubs.o geom_helper.o precision.o sys.o
 m_hsx.o: atm_types.o atmfuncs.o atomlist.o files.o parallel.o parallelsubs.o
 m_hsx.o: precision.o siesta_geom.o sys.o
+m_initwf.o: alloc.o densematrix.o fermid.o m_diagon_opt.o m_eo.o m_memory.o
+m_initwf.o: parallel.o parallelsubs.o precision.o sparse_matrices.o sys.o
+m_initwf.o: wavefunctions.o
 m_integrate.o: precision.o
 m_io.o: sys.o
 m_io_s.o: class_OrbitalDistribution.o class_SpData1D.o class_SpData2D.o
 m_io_s.o: class_Sparsity.o parallel.o precision.o
+m_iodipol.o: files.o m_dipol.o m_steps.o
 m_iodm.o: class_OrbitalDistribution.o class_SpData2D.o class_Sparsity.o
 m_iodm.o: m_io_s.o parallel.o
 m_iodm_old.o: alloc.o files.o m_matio.o m_mpi_utils.o parallel.o parallelsubs.o
@@ -415,9 +427,13 @@
 m_iorho.o: alloc.o parallel.o parallelsubs.o precision.o sys.o
 m_iostruct.o: alloc.o files.o m_mpi_utils.o parallel.o precision.o
 m_iostruct.o: siesta_geom.o sys.o units.o
+m_iotddft.o: files.o m_dipol.o m_steps.o parallel.o siesta_options.o
+m_iotddft.o: wavefunctions.o
 m_kinetic.o: precision.o
 m_mat_invert.o: intrinsic_missing.o m_pivot_array.o precision.o
+m_matdiag.o: precision.o
 m_matio.o: alloc.o parallel.o
+m_matswinvers.o: precision.o
 m_memory.o: memoryinfo.o parallel.o sys.o
 m_mesh_node.o: intrinsic_missing.o parallel.o precision.o
 m_mixing.o: class_Data1D.o class_Fstack_Data1D.o parallel.o precision.o
@@ -465,6 +481,15 @@
 m_planewavematrix.o: atomlist.o dhscf.o m_planewavematrixvar.o precision.o
 m_planewavematrix.o: siesta_geom.o sparse_matrices.o
 m_planewavematrixvar.o: precision.o
+m_plot_1d_arr3d.o: precision.o
+m_poisson_fdf.o: m_char.o precision.o
+m_poisson_obc_comm.o: atm_types.o atmfuncs.o m_mpi_stop.o parallel.o
+m_poisson_obc_comm.o: precision.o siesta_geom.o sys.o
+m_poisson_obc_fastconv.o: cellsubs.o m_mpi_stop.o m_poisson_obc_fastconv_r2c.o
+m_poisson_obc_fastconv.o: parallel.o precision.o sys.o
+m_poisson_obc_fastconv_r2c.o: cellsubs.o m_fastconv_dims.o m_plot_1d_arr3d.o
+m_poisson_obc_fastconv_r2c.o: m_poisson_obc_comm.o precision.o sys.o
+m_poisson_switch.o: m_poisson_obc_fastconv.o precision.o sys.o
 m_pulay.o: alloc.o atomlist.o m_mpi_utils.o m_spin.o m_svd.o parallel.o
 m_pulay.o: precision.o siesta_options.o sparse_matrices.o sys.o
 m_redist_spmatrix.o: alloc.o class_Distribution.o
@@ -480,7 +505,8 @@
 m_sparsity_handling.o: class_Data1D.o class_Data2D.o
 m_sparsity_handling.o: class_OrbitalDistribution.o class_SpData1D.o
 m_sparsity_handling.o: class_SpData2D.o class_Sparsity.o geom_helper.o
-m_sparsity_handling.o: m_interpolate.o m_region.o precision.o
+m_sparsity_handling.o: intrinsic_missing.o m_interpolate.o m_region.o
+m_sparsity_handling.o: precision.o
 m_spin.o: alloc.o parallel.o precision.o sys.o units.o
 m_stress.o: precision.o
 m_supercell.o: atom_graph.o class_OrbitalDistribution.o class_SpData2D.o
@@ -592,9 +618,10 @@
 m_ts_options.o: m_ts_hartree.o m_ts_method.o m_ts_mumps_init.o m_ts_tdir.o
 m_ts_options.o: m_ts_weight.o parallel.o precision.o siesta_options.o units.o
 m_ts_pivot.o: alloc.o class_OrbitalDistribution.o class_Sparsity.o
-m_ts_pivot.o: create_Sparsity_SC.o create_Sparsity_Union.o intrinsic_missing.o
-m_ts_pivot.o: m_char.o m_pivot.o m_region.o m_sparsity_handling.o m_ts_debug.o
-m_ts_pivot.o: m_ts_electype.o parallel.o precision.o
+m_ts_pivot.o: create_Sparsity_SC.o create_Sparsity_Union.o geom_helper.o
+m_ts_pivot.o: intrinsic_missing.o m_char.o m_pivot.o m_region.o
+m_ts_pivot.o: m_sparsity_handling.o m_ts_debug.o m_ts_electype.o parallel.o
+m_ts_pivot.o: precision.o
 m_ts_rgn2trimat.o: alloc.o class_OrbitalDistribution.o class_Sparsity.o
 m_ts_rgn2trimat.o: create_Sparsity_Union.o geom_helper.o m_region.o
 m_ts_rgn2trimat.o: m_ts_electype.o m_ts_method.o m_ts_tri_common.o parallel.o
@@ -737,8 +764,8 @@
 rdiag.o: alloc.o m_diagon_opt.o memoryinfo.o parallel.o precision.o sys.o
 rdiag_elpa.o: alloc.o elpa1.o elpa2.o parallel.o precision.o sys.o
 rdiag_mrrr.o: alloc.o parallel.o precision.o sys.o
-read_options.o: files.o m_charge_add.o m_hartree_add.o m_mixing_scf.o m_spin.o
-read_options.o: m_target_stress.o memoryinfo.o parallel.o precision.o
+read_options.o: files.o m_charge_add.o m_cite.o m_hartree_add.o m_mixing_scf.o
+read_options.o: m_spin.o m_target_stress.o memoryinfo.o parallel.o precision.o
 read_options.o: siesta_cml.o siesta_options.o sys.o units.o
 read_xc_info.o: parallel.o precision.o sys.o
 readsp.o: parallel.o precision.o sys.o
@@ -805,30 +832,32 @@
 siesta_end.o: meshdscf.o meshphi.o moremeshsubs.o parallel.o siesta_cml.o
 siesta_end.o: siesta_cmlsubs.o siesta_master.o siesta_options.o
 siesta_end.o: sparse_matrices.o timestamp.o
-siesta_forces.o: atomlist.o class_Fstack_Data1D.o class_SpData2D.o compute_dm.o
-siesta_forces.o: compute_energies.o compute_max_diff.o files.o
-siesta_forces.o: final_H_f_stress.o flook_siesta.o kpoint_grid.o
-siesta_forces.o: m_convergence.o m_energies.o m_forces.o m_iodm_old.o
-siesta_forces.o: m_mixing.o m_mixing_scf.o m_mpi_utils.o m_ncdf_siesta.o
-siesta_forces.o: m_pexsi.o m_pexsi_driver.o m_rhog.o m_spin.o m_steps.o
-siesta_forces.o: m_stress.o m_transiesta.o m_ts_charge.o m_ts_electype.o
-siesta_forces.o: m_ts_global_vars.o m_ts_method.o m_ts_options.o mixer.o
-siesta_forces.o: parallel.o post_scf_work.o precision.o save_density_matrix.o
-siesta_forces.o: scfconvergence_test.o setup_H0.o setup_hamiltonian.o
-siesta_forces.o: siesta_cml.o siesta_dicts.o siesta_geom.o siesta_master.o
-siesta_forces.o: siesta_options.o sparse_matrices.o state_analysis.o
-siesta_forces.o: state_init.o sys.o units.o write_subs.o
+siesta_forces.o: alloc.o atomlist.o class_Fstack_Data1D.o class_SpData2D.o
+siesta_forces.o: compute_dm.o compute_energies.o compute_max_diff.o
+siesta_forces.o: densematrix.o files.o final_H_f_stress.o flook_siesta.o
+siesta_forces.o: kpoint_grid.o m_convergence.o m_energies.o m_eo.o m_evolve.o
+siesta_forces.o: m_forces.o m_gamma.o m_initwf.o m_iodm.o m_iodm_old.o
+siesta_forces.o: m_iotddft.o m_mixing.o m_mixing_scf.o m_mpi_utils.o
+siesta_forces.o: m_ncdf_siesta.o m_pexsi.o m_pexsi_driver.o m_rhog.o m_spin.o
+siesta_forces.o: m_steps.o m_stress.o m_transiesta.o m_ts_charge.o
+siesta_forces.o: m_ts_electype.o m_ts_global_vars.o m_ts_method.o
+siesta_forces.o: m_ts_options.o mixer.o overfsm.o parallel.o post_scf_work.o
+siesta_forces.o: precision.o save_density_matrix.o scfconvergence_test.o
+siesta_forces.o: setup_H0.o setup_hamiltonian.o siesta_cml.o siesta_dicts.o
+siesta_forces.o: siesta_geom.o siesta_master.o siesta_options.o
+siesta_forces.o: sparse_matrices.o state_analysis.o state_init.o sys.o units.o
+siesta_forces.o: wavefunctions.o write_subs.o
 siesta_geom.o: precision.o
 siesta_init.o: alloc.o atomlist.o bands.o bsc_xcmod.o
 siesta_init.o: class_Fstack_Pair_Geometry_SpData2D.o densematrix.o files.o
 siesta_init.o: flook_siesta.o ioxv.o kpoint_grid.o kpoint_pdos.o ksvinit.o
 siesta_init.o: m_cite.o m_energies.o m_eo.o m_fixed.o m_forces.o m_gamma.o
-siesta_init.o: m_iostruct.o m_mpi_utils.o m_new_dm.o m_rmaxh.o m_spin.o
-siesta_init.o: m_steps.o m_supercell.o m_timer.o m_wallclock.o metaforce.o
-siesta_init.o: molecularmechanics.o object_debug.o parallel.o parallelsubs.o
-siesta_init.o: projected_DOS.o siesta_cmlsubs.o siesta_dicts.o siesta_geom.o
-siesta_init.o: siesta_options.o sparse_matrices.o struct_init.o sys.o timer.o
-siesta_init.o: timestamp.o ts_init.o writewave.o zmatrix.o
+siesta_init.o: m_iostruct.o m_mpi_utils.o m_new_dm.o m_poisson_fdf.o m_rmaxh.o
+siesta_init.o: m_spin.o m_steps.o m_supercell.o m_timer.o m_wallclock.o
+siesta_init.o: metaforce.o molecularmechanics.o object_debug.o parallel.o
+siesta_init.o: parallelsubs.o projected_DOS.o siesta_cmlsubs.o siesta_dicts.o
+siesta_init.o: siesta_geom.o siesta_options.o sparse_matrices.o struct_init.o
+siesta_init.o: sys.o timer.o timestamp.o ts_init.o writewave.o zmatrix.o
 siesta_master.o: iopipes.o iosockets.o precision.o sys.o
 siesta_move.o: atomlist.o broyden_optim.o cell_broyden_optim.o
 siesta_move.o: cell_fire_optim.o dynamics.o fire_optim.o flook_siesta.o ioxv.o
@@ -860,18 +889,17 @@
 state_init.o: proximity_check.o siesta_cml.o siesta_geom.o siesta_options.o
 state_init.o: sparse_matrices.o sys.o units.o write_subs.o zmatrix.o
 struct_init.o: alloc.o atmfuncs.o atomlist.o files.o ioxv.o m_exp_coord.o
-struct_init.o: m_iostruct.o m_mpi_utils.o m_steps.o parallel.o periodic_table.o
-struct_init.o: siesta_cml.o siesta_geom.o siesta_master.o siesta_options.o
-struct_init.o: units.o zmatrix.o
+struct_init.o: m_iostruct.o m_mpi_utils.o m_poisson_fdf.o m_poisson_obc_comm.o
+struct_init.o: m_steps.o parallel.o periodic_table.o siesta_cml.o siesta_geom.o
+struct_init.o: siesta_master.o siesta_options.o units.o zmatrix.o
 sys.o: parallel.o
 timer.o: extrae_eventllist.o extrae_module.o m_timer.o parallel.o sys.o
 timer.o: timer_tree.o
 timer_tree.o: m_walltime.o
 transition_rate.o: alloc.o fermid.o parallel.o parallelsubs.o precision.o sys.o
-ts_init.o: m_cite.o m_fixed.o m_os.o m_ts_cctype.o m_ts_charge.o
-ts_init.o: m_ts_electrode.o m_ts_electype.o m_ts_gf.o m_ts_global_vars.o
-ts_init.o: m_ts_kpoints.o m_ts_method.o m_ts_options.o parallel.o
-ts_init.o: siesta_options.o
+ts_init.o: m_fixed.o m_os.o m_ts_cctype.o m_ts_charge.o m_ts_electrode.o
+ts_init.o: m_ts_electype.o m_ts_gf.o m_ts_global_vars.o m_ts_kpoints.o
+ts_init.o: m_ts_method.o m_ts_options.o parallel.o siesta_options.o
 ts_show_regions.o: m_region.o m_ts_electype.o m_ts_method.o parallel.o
 ts_show_regions.o: precision.o units.o
 typecell.o: precision.o
@@ -883,6 +911,7 @@
 vmatsp.o: alloc.o atm_types.o atmfuncs.o atomlist.o listsc.o mesh.o meshdscf.o
 vmatsp.o: meshphi.o precision.o
 vmb.o: m_fixed.o parallel.o precision.o sys.o
+wavefunctions.o: parallel.o precision.o
 write_inp_wannier.o: alloc.o atmfuncs.o atomlist.o m_ntm.o m_orderbands.o
 write_inp_wannier.o: mneighb.o parallel.o parallelsubs.o precision.o
 write_inp_wannier.o: siesta2wannier90.o siesta_geom.o
@@ -1004,6 +1033,7 @@
 m_ksv.o: ksv.o
 m_ksvinit.o: ksvinit.o
 m_ldau.o: ldau.o
+m_linpack.o: linpack.o
 m_local_dos.o: local_DOS.o
 m_matel_registry.o: matel_registry.o
 m_minvec.o: minvec.o
@@ -1048,6 +1078,8 @@
 m_vmat.o: vmat.o
 m_zm_broyden_optim.o: zm_broyden_optim.o
 m_zm_fire_optim.o: zm_fire_optim.o
+matdiagon.o: m_matdiag.o
+matswinversion.o: m_matswinvers.o
 mtprng.o: m_uuid.o
 neighbour.o: mneighb.o
 on_core.o: onmod.o

=== modified file 'Util/Grimme/Makefile'
--- Util/Grimme/Makefile	2016-08-30 14:02:33 +0000
+++ Util/Grimme/Makefile	2016-11-22 08:16:45 +0000
@@ -109,7 +109,7 @@
 bsc_cellxc.o: alloc.o bsc_xcmod.o cellxc_mod.o mesh.o moremeshsubs.o parallel.o
 bsc_cellxc.o: parallelsubs.o precision.o sys.o
 bsc_xcmod.o: parallel.o precision.o sys.o
-cart2frac.o: sys.o
+cart2frac.o: precision.o sys.o
 cdiag.o: alloc.o m_diagon_opt.o memoryinfo.o parallel.o precision.o sys.o
 cell_broyden_optim.o: m_broyddj_nocomm.o parallel.o precision.o sys.o units.o
 cell_broyden_optim.o: zmatrix.o
@@ -121,6 +121,8 @@
 cgvc_zmatrix.o: alloc.o conjgr.o m_mpi_utils.o parallel.o precision.o sys.o
 cgvc_zmatrix.o: units.o zmatrix.o
 cgwf.o: egandd.o onmod.o parallel.o precision.o sys.o
+changebasis.o: alloc.o m_matdiag.o parallel.o parallelsubs.o precision.o
+changebasis.o: sparse_matrices.o sys.o wavefunctions.o
 chemical.o: parallel.o precision.o sys.o
 chempot.o: alloc.o mneighb.o parallel.o parallelsubs.o precision.o sys.o
 chkdim.o: sys.o
@@ -172,17 +174,17 @@
 denmat.o: alloc.o globalise.o m_mpi_utils.o onmod.o precision.o
 denmatlomem.o: alloc.o globalise.o onmod.o precision.o
 densematrix.o: alloc.o precision.o
-detover.o: alloc.o parallel.o parallelsubs.o precision.o
+detover.o: alloc.o linpack.o parallel.o parallelsubs.o precision.o
 dfscf.o: alloc.o atm_types.o atmfuncs.o atomlist.o listsc.o m_spin.o mesh.o
 dfscf.o: meshdscf.o meshphi.o parallel.o parallelsubs.o precision.o sys.o
 dhscf.o: alloc.o atmfuncs.o bsc_xcmod.o cellxc_mod.o delk.o dfscf.o
 dhscf.o: doping_uniform.o files.o forhar.o iogrid_netcdf.o m_charge_add.o
-dhscf.o: m_efield.o m_hartree_add.o m_iorho.o m_mesh_node.o m_ncdf_io.o
-dhscf.o: m_ncdf_siesta.o m_partial_charges.o m_rhog.o m_spin.o
-dhscf.o: m_ts_global_vars.o m_ts_hartree.o m_ts_options.o m_ts_tdir.o
-dhscf.o: m_ts_voltage.o mesh.o meshdscf.o meshsubs.o moremeshsubs.o parallel.o
-dhscf.o: parsing.o precision.o rhofft.o rhoofd.o siesta_options.o sys.o units.o
-dhscf.o: vmat.o
+dhscf.o: m_efield.o m_hartree_add.o m_iorho.o m_iotddft.o m_mesh_node.o
+dhscf.o: m_ncdf_io.o m_ncdf_siesta.o m_partial_charges.o m_poisson_fdf.o
+dhscf.o: m_poisson_switch.o m_rhog.o m_spin.o m_ts_global_vars.o m_ts_hartree.o
+dhscf.o: m_ts_options.o m_ts_tdir.o m_ts_voltage.o mesh.o meshdscf.o meshsubs.o
+dhscf.o: moremeshsubs.o parallel.o parsing.o precision.o rhofft.o rhoofd.o
+dhscf.o: siesta_options.o sys.o units.o vmat.o
 diag2g.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
 diag2k.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
 diag3g.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
@@ -210,7 +212,7 @@
 domain_decom.o: sparse_matrices.o sys.o
 doping_uniform.o: alloc.o m_ntm.o mesh.o parallel.o precision.o sys.o
 dynamics.o: alloc.o atomlist.o files.o ioxv.o m_mpi_utils.o parallel.o
-dynamics.o: precision.o sys.o units.o
+dynamics.o: precision.o siesta_options.o sys.o units.o
 egandd.o: alloc.o denmat.o ener3.o globalise.o gradient.o on_subs.o onmod.o
 egandd.o: onmod.o precision.o sys.o
 eggbox.o: parallel.o precision.o
@@ -219,6 +221,10 @@
 elpa2.o: elpa1.o
 ener3.o: alloc.o globalise.o m_mpi_utils.o onmod.o precision.o
 ener3lomem.o: alloc.o globalise.o m_mpi_utils.o onmod.o precision.o
+evolg.o: alloc.o m_eo.o m_matswinvers.o m_steps.o parallel.o parallelsubs.o
+evolg.o: precision.o siesta_options.o sparse_matrices.o sys.o wavefunctions.o
+evolk.o: m_diagon_opt.o m_matswinvers.o parallel.o precision.o siesta_options.o
+evolk.o: sys.o wavefunctions.o
 extrapolateSpData2D.o: class_Data2D.o class_OrbitalDistribution.o
 extrapolateSpData2D.o: class_SpData2D.o class_Sparsity.o
 extrapolon.o: parallel.o precision.o sys.o
@@ -259,9 +265,9 @@
 initatom.o: atmparams.o atom.o atom_options.o basis_io.o basis_specs.o
 initatom.o: basis_types.o chemical.o electrostatic.o ldau_specs.o m_spin.o
 initatom.o: old_atmfuncs.o precision.o pseudopotential.o
-initparallel.o: alloc.o class_OrbitalDistribution.o domain_decom.o parallel.o
-initparallel.o: parallelsubs.o precision.o siesta_options.o sparse_matrices.o
-initparallel.o: spatial.o sys.o
+initparallel.o: alloc.o class_OrbitalDistribution.o domain_decom.o
+initparallel.o: m_poisson_fdf.o parallel.o parallelsubs.o precision.o
+initparallel.o: siesta_options.o sparse_matrices.o spatial.o sys.o
 io.o: sys.o
 iocg.o: files.o precision.o sys.o
 iodm.o: alloc.o files.o parallel.o parallelsubs.o precision.o sys.o
@@ -320,7 +326,9 @@
 m_efield.o: siesta_geom.o sys.o units.o
 m_energies.o: precision.o
 m_eo.o: precision.o
+m_evolve.o: alloc.o m_memory.o parallel.o parallelsubs.o precision.o sys.o
 m_exp_coord.o: files.o m_os.o parallel.o precision.o units.o
+m_fastconv_dims.o: precision.o
 m_filter.o: bessph.o precision.o radfft.o sorting.o sys.o
 m_fire.o: parallel.o precision.o
 m_fire_mixing.o: alloc.o m_fire_para.o m_mpi_utils.o parallel.o precision.o
@@ -342,10 +350,14 @@
 m_hs_matrix.o: alloc.o cellsubs.o geom_helper.o precision.o sys.o
 m_hsx.o: atm_types.o atmfuncs.o atomlist.o files.o parallel.o parallelsubs.o
 m_hsx.o: precision.o siesta_geom.o sys.o
+m_initwf.o: alloc.o densematrix.o fermid.o m_diagon_opt.o m_eo.o m_memory.o
+m_initwf.o: parallel.o parallelsubs.o precision.o sparse_matrices.o sys.o
+m_initwf.o: wavefunctions.o
 m_integrate.o: precision.o
 m_io.o: sys.o
 m_io_s.o: class_OrbitalDistribution.o class_SpData1D.o class_SpData2D.o
 m_io_s.o: class_Sparsity.o parallel.o precision.o
+m_iodipol.o: files.o m_dipol.o m_steps.o
 m_iodm.o: class_OrbitalDistribution.o class_SpData2D.o class_Sparsity.o
 m_iodm.o: m_io_s.o parallel.o
 m_iodm_old.o: alloc.o files.o m_matio.o m_mpi_utils.o parallel.o parallelsubs.o
@@ -353,9 +365,13 @@
 m_iorho.o: alloc.o parallel.o parallelsubs.o precision.o sys.o
 m_iostruct.o: alloc.o files.o m_mpi_utils.o parallel.o precision.o
 m_iostruct.o: siesta_geom.o sys.o units.o
+m_iotddft.o: files.o m_dipol.o m_steps.o parallel.o siesta_options.o
+m_iotddft.o: wavefunctions.o
 m_kinetic.o: precision.o
 m_mat_invert.o: intrinsic_missing.o m_pivot_array.o precision.o
+m_matdiag.o: precision.o
 m_matio.o: alloc.o parallel.o
+m_matswinvers.o: precision.o
 m_memory.o: memoryinfo.o parallel.o sys.o
 m_mesh_node.o: intrinsic_missing.o parallel.o precision.o
 m_mixing.o: class_Data1D.o class_Fstack_Data1D.o parallel.o precision.o
@@ -403,6 +419,15 @@
 m_planewavematrix.o: atomlist.o dhscf.o m_planewavematrixvar.o precision.o
 m_planewavematrix.o: siesta_geom.o sparse_matrices.o
 m_planewavematrixvar.o: precision.o
+m_plot_1d_arr3d.o: precision.o
+m_poisson_fdf.o: m_char.o precision.o
+m_poisson_obc_comm.o: atm_types.o atmfuncs.o m_mpi_stop.o parallel.o
+m_poisson_obc_comm.o: precision.o siesta_geom.o sys.o
+m_poisson_obc_fastconv.o: cellsubs.o m_mpi_stop.o m_poisson_obc_fastconv_r2c.o
+m_poisson_obc_fastconv.o: parallel.o precision.o sys.o
+m_poisson_obc_fastconv_r2c.o: cellsubs.o m_fastconv_dims.o m_plot_1d_arr3d.o
+m_poisson_obc_fastconv_r2c.o: m_poisson_obc_comm.o precision.o sys.o
+m_poisson_switch.o: m_poisson_obc_fastconv.o precision.o sys.o
 m_pulay.o: alloc.o atomlist.o m_mpi_utils.o m_spin.o m_svd.o parallel.o
 m_pulay.o: precision.o siesta_options.o sparse_matrices.o sys.o
 m_redist_spmatrix.o: alloc.o class_Distribution.o
@@ -418,7 +443,8 @@
 m_sparsity_handling.o: class_Data1D.o class_Data2D.o
 m_sparsity_handling.o: class_OrbitalDistribution.o class_SpData1D.o
 m_sparsity_handling.o: class_SpData2D.o class_Sparsity.o geom_helper.o
-m_sparsity_handling.o: m_interpolate.o m_region.o precision.o
+m_sparsity_handling.o: intrinsic_missing.o m_interpolate.o m_region.o
+m_sparsity_handling.o: precision.o
 m_spin.o: alloc.o parallel.o precision.o sys.o units.o
 m_stress.o: precision.o
 m_supercell.o: atom_graph.o class_OrbitalDistribution.o class_SpData2D.o
@@ -530,9 +556,10 @@
 m_ts_options.o: m_ts_hartree.o m_ts_method.o m_ts_mumps_init.o m_ts_tdir.o
 m_ts_options.o: m_ts_weight.o parallel.o precision.o siesta_options.o units.o
 m_ts_pivot.o: alloc.o class_OrbitalDistribution.o class_Sparsity.o
-m_ts_pivot.o: create_Sparsity_SC.o create_Sparsity_Union.o intrinsic_missing.o
-m_ts_pivot.o: m_char.o m_pivot.o m_region.o m_sparsity_handling.o m_ts_debug.o
-m_ts_pivot.o: m_ts_electype.o parallel.o precision.o
+m_ts_pivot.o: create_Sparsity_SC.o create_Sparsity_Union.o geom_helper.o
+m_ts_pivot.o: intrinsic_missing.o m_char.o m_pivot.o m_region.o
+m_ts_pivot.o: m_sparsity_handling.o m_ts_debug.o m_ts_electype.o parallel.o
+m_ts_pivot.o: precision.o
 m_ts_rgn2trimat.o: alloc.o class_OrbitalDistribution.o class_Sparsity.o
 m_ts_rgn2trimat.o: create_Sparsity_Union.o geom_helper.o m_region.o
 m_ts_rgn2trimat.o: m_ts_electype.o m_ts_method.o m_ts_tri_common.o parallel.o
@@ -675,8 +702,8 @@
 rdiag.o: alloc.o m_diagon_opt.o memoryinfo.o parallel.o precision.o sys.o
 rdiag_elpa.o: alloc.o elpa1.o elpa2.o parallel.o precision.o sys.o
 rdiag_mrrr.o: alloc.o parallel.o precision.o sys.o
-read_options.o: files.o m_charge_add.o m_hartree_add.o m_mixing_scf.o m_spin.o
-read_options.o: m_target_stress.o memoryinfo.o parallel.o precision.o
+read_options.o: files.o m_charge_add.o m_cite.o m_hartree_add.o m_mixing_scf.o
+read_options.o: m_spin.o m_target_stress.o memoryinfo.o parallel.o precision.o
 read_options.o: siesta_cml.o siesta_options.o sys.o units.o
 read_xc_info.o: parallel.o precision.o sys.o
 readsp.o: parallel.o precision.o sys.o
@@ -743,30 +770,32 @@
 siesta_end.o: meshdscf.o meshphi.o moremeshsubs.o parallel.o siesta_cml.o
 siesta_end.o: siesta_cmlsubs.o siesta_master.o siesta_options.o
 siesta_end.o: sparse_matrices.o timestamp.o
-siesta_forces.o: atomlist.o class_Fstack_Data1D.o class_SpData2D.o compute_dm.o
-siesta_forces.o: compute_energies.o compute_max_diff.o files.o
-siesta_forces.o: final_H_f_stress.o flook_siesta.o kpoint_grid.o
-siesta_forces.o: m_convergence.o m_energies.o m_forces.o m_iodm_old.o
-siesta_forces.o: m_mixing.o m_mixing_scf.o m_mpi_utils.o m_ncdf_siesta.o
-siesta_forces.o: m_pexsi.o m_pexsi_driver.o m_rhog.o m_spin.o m_steps.o
-siesta_forces.o: m_stress.o m_transiesta.o m_ts_charge.o m_ts_electype.o
-siesta_forces.o: m_ts_global_vars.o m_ts_method.o m_ts_options.o mixer.o
-siesta_forces.o: parallel.o post_scf_work.o precision.o save_density_matrix.o
-siesta_forces.o: scfconvergence_test.o setup_H0.o setup_hamiltonian.o
-siesta_forces.o: siesta_cml.o siesta_dicts.o siesta_geom.o siesta_master.o
-siesta_forces.o: siesta_options.o sparse_matrices.o state_analysis.o
-siesta_forces.o: state_init.o sys.o units.o write_subs.o
+siesta_forces.o: alloc.o atomlist.o class_Fstack_Data1D.o class_SpData2D.o
+siesta_forces.o: compute_dm.o compute_energies.o compute_max_diff.o
+siesta_forces.o: densematrix.o files.o final_H_f_stress.o flook_siesta.o
+siesta_forces.o: kpoint_grid.o m_convergence.o m_energies.o m_eo.o m_evolve.o
+siesta_forces.o: m_forces.o m_gamma.o m_initwf.o m_iodm.o m_iodm_old.o
+siesta_forces.o: m_iotddft.o m_mixing.o m_mixing_scf.o m_mpi_utils.o
+siesta_forces.o: m_ncdf_siesta.o m_pexsi.o m_pexsi_driver.o m_rhog.o m_spin.o
+siesta_forces.o: m_steps.o m_stress.o m_transiesta.o m_ts_charge.o
+siesta_forces.o: m_ts_electype.o m_ts_global_vars.o m_ts_method.o
+siesta_forces.o: m_ts_options.o mixer.o overfsm.o parallel.o post_scf_work.o
+siesta_forces.o: precision.o save_density_matrix.o scfconvergence_test.o
+siesta_forces.o: setup_H0.o setup_hamiltonian.o siesta_cml.o siesta_dicts.o
+siesta_forces.o: siesta_geom.o siesta_master.o siesta_options.o
+siesta_forces.o: sparse_matrices.o state_analysis.o state_init.o sys.o units.o
+siesta_forces.o: wavefunctions.o write_subs.o
 siesta_geom.o: precision.o
 siesta_init.o: alloc.o atomlist.o bands.o bsc_xcmod.o
 siesta_init.o: class_Fstack_Pair_Geometry_SpData2D.o densematrix.o files.o
 siesta_init.o: flook_siesta.o ioxv.o kpoint_grid.o kpoint_pdos.o ksvinit.o
 siesta_init.o: m_cite.o m_energies.o m_eo.o m_fixed.o m_forces.o m_gamma.o
-siesta_init.o: m_iostruct.o m_mpi_utils.o m_new_dm.o m_rmaxh.o m_spin.o
-siesta_init.o: m_steps.o m_supercell.o m_timer.o m_wallclock.o metaforce.o
-siesta_init.o: molecularmechanics.o object_debug.o parallel.o parallelsubs.o
-siesta_init.o: projected_DOS.o siesta_cmlsubs.o siesta_dicts.o siesta_geom.o
-siesta_init.o: siesta_options.o sparse_matrices.o struct_init.o sys.o timer.o
-siesta_init.o: timestamp.o ts_init.o writewave.o zmatrix.o
+siesta_init.o: m_iostruct.o m_mpi_utils.o m_new_dm.o m_poisson_fdf.o m_rmaxh.o
+siesta_init.o: m_spin.o m_steps.o m_supercell.o m_timer.o m_wallclock.o
+siesta_init.o: metaforce.o molecularmechanics.o object_debug.o parallel.o
+siesta_init.o: parallelsubs.o projected_DOS.o siesta_cmlsubs.o siesta_dicts.o
+siesta_init.o: siesta_geom.o siesta_options.o sparse_matrices.o struct_init.o
+siesta_init.o: sys.o timer.o timestamp.o ts_init.o writewave.o zmatrix.o
 siesta_master.o: iopipes.o iosockets.o precision.o sys.o
 siesta_move.o: atomlist.o broyden_optim.o cell_broyden_optim.o
 siesta_move.o: cell_fire_optim.o dynamics.o fire_optim.o flook_siesta.o ioxv.o
@@ -798,18 +827,17 @@
 state_init.o: proximity_check.o siesta_cml.o siesta_geom.o siesta_options.o
 state_init.o: sparse_matrices.o sys.o units.o write_subs.o zmatrix.o
 struct_init.o: alloc.o atmfuncs.o atomlist.o files.o ioxv.o m_exp_coord.o
-struct_init.o: m_iostruct.o m_mpi_utils.o m_steps.o parallel.o periodic_table.o
-struct_init.o: siesta_cml.o siesta_geom.o siesta_master.o siesta_options.o
-struct_init.o: units.o zmatrix.o
+struct_init.o: m_iostruct.o m_mpi_utils.o m_poisson_fdf.o m_poisson_obc_comm.o
+struct_init.o: m_steps.o parallel.o periodic_table.o siesta_cml.o siesta_geom.o
+struct_init.o: siesta_master.o siesta_options.o units.o zmatrix.o
 sys.o: parallel.o
 timer.o: extrae_eventllist.o extrae_module.o m_timer.o parallel.o sys.o
 timer.o: timer_tree.o
 timer_tree.o: m_walltime.o
 transition_rate.o: alloc.o fermid.o parallel.o parallelsubs.o precision.o sys.o
-ts_init.o: m_cite.o m_fixed.o m_os.o m_ts_cctype.o m_ts_charge.o
-ts_init.o: m_ts_electrode.o m_ts_electype.o m_ts_gf.o m_ts_global_vars.o
-ts_init.o: m_ts_kpoints.o m_ts_method.o m_ts_options.o parallel.o
-ts_init.o: siesta_options.o
+ts_init.o: m_fixed.o m_os.o m_ts_cctype.o m_ts_charge.o m_ts_electrode.o
+ts_init.o: m_ts_electype.o m_ts_gf.o m_ts_global_vars.o m_ts_kpoints.o
+ts_init.o: m_ts_method.o m_ts_options.o parallel.o siesta_options.o
 ts_show_regions.o: m_region.o m_ts_electype.o m_ts_method.o parallel.o
 ts_show_regions.o: precision.o units.o
 typecell.o: precision.o
@@ -821,6 +849,7 @@
 vmatsp.o: alloc.o atm_types.o atmfuncs.o atomlist.o listsc.o mesh.o meshdscf.o
 vmatsp.o: meshphi.o precision.o
 vmb.o: m_fixed.o parallel.o precision.o sys.o
+wavefunctions.o: parallel.o precision.o
 write_inp_wannier.o: alloc.o atmfuncs.o atomlist.o m_ntm.o m_orderbands.o
 write_inp_wannier.o: mneighb.o parallel.o parallelsubs.o precision.o
 write_inp_wannier.o: siesta2wannier90.o siesta_geom.o
@@ -938,6 +967,7 @@
 m_ksv.o: ksv.o
 m_ksvinit.o: ksvinit.o
 m_ldau.o: ldau.o
+m_linpack.o: linpack.o
 m_local_dos.o: local_DOS.o
 m_matel_registry.o: matel_registry.o
 m_minvec.o: minvec.o
@@ -982,6 +1012,8 @@
 m_vmat.o: vmat.o
 m_zm_broyden_optim.o: zm_broyden_optim.o
 m_zm_fire_optim.o: zm_fire_optim.o
+matdiagon.o: m_matdiag.o
+matswinversion.o: m_matswinvers.o
 mtprng.o: m_uuid.o
 neighbour.o: mneighb.o
 on_core.o: onmod.o

=== modified file 'Util/Helpers/Makefile'
--- Util/Helpers/Makefile	2016-08-30 14:02:33 +0000
+++ Util/Helpers/Makefile	2016-11-22 08:16:45 +0000
@@ -112,7 +112,7 @@
 bsc_cellxc.o: alloc.o bsc_xcmod.o cellxc_mod.o mesh.o moremeshsubs.o parallel.o
 bsc_cellxc.o: parallelsubs.o precision.o sys.o
 bsc_xcmod.o: parallel.o precision.o sys.o
-cart2frac.o: sys.o
+cart2frac.o: precision.o sys.o
 cdiag.o: alloc.o m_diagon_opt.o memoryinfo.o parallel.o precision.o sys.o
 cell_broyden_optim.o: m_broyddj_nocomm.o parallel.o precision.o sys.o units.o
 cell_broyden_optim.o: zmatrix.o
@@ -124,6 +124,8 @@
 cgvc_zmatrix.o: alloc.o conjgr.o m_mpi_utils.o parallel.o precision.o sys.o
 cgvc_zmatrix.o: units.o zmatrix.o
 cgwf.o: egandd.o onmod.o parallel.o precision.o sys.o
+changebasis.o: alloc.o m_matdiag.o parallel.o parallelsubs.o precision.o
+changebasis.o: sparse_matrices.o sys.o wavefunctions.o
 chemical.o: parallel.o precision.o sys.o
 chempot.o: alloc.o mneighb.o parallel.o parallelsubs.o precision.o sys.o
 chkdim.o: sys.o
@@ -175,17 +177,17 @@
 denmat.o: alloc.o globalise.o m_mpi_utils.o onmod.o precision.o
 denmatlomem.o: alloc.o globalise.o onmod.o precision.o
 densematrix.o: alloc.o precision.o
-detover.o: alloc.o parallel.o parallelsubs.o precision.o
+detover.o: alloc.o linpack.o parallel.o parallelsubs.o precision.o
 dfscf.o: alloc.o atm_types.o atmfuncs.o atomlist.o listsc.o m_spin.o mesh.o
 dfscf.o: meshdscf.o meshphi.o parallel.o parallelsubs.o precision.o sys.o
 dhscf.o: alloc.o atmfuncs.o bsc_xcmod.o cellxc_mod.o delk.o dfscf.o
 dhscf.o: doping_uniform.o files.o forhar.o iogrid_netcdf.o m_charge_add.o
-dhscf.o: m_efield.o m_hartree_add.o m_iorho.o m_mesh_node.o m_ncdf_io.o
-dhscf.o: m_ncdf_siesta.o m_partial_charges.o m_rhog.o m_spin.o
-dhscf.o: m_ts_global_vars.o m_ts_hartree.o m_ts_options.o m_ts_tdir.o
-dhscf.o: m_ts_voltage.o mesh.o meshdscf.o meshsubs.o moremeshsubs.o parallel.o
-dhscf.o: parsing.o precision.o rhofft.o rhoofd.o siesta_options.o sys.o units.o
-dhscf.o: vmat.o
+dhscf.o: m_efield.o m_hartree_add.o m_iorho.o m_iotddft.o m_mesh_node.o
+dhscf.o: m_ncdf_io.o m_ncdf_siesta.o m_partial_charges.o m_poisson_fdf.o
+dhscf.o: m_poisson_switch.o m_rhog.o m_spin.o m_ts_global_vars.o m_ts_hartree.o
+dhscf.o: m_ts_options.o m_ts_tdir.o m_ts_voltage.o mesh.o meshdscf.o meshsubs.o
+dhscf.o: moremeshsubs.o parallel.o parsing.o precision.o rhofft.o rhoofd.o
+dhscf.o: siesta_options.o sys.o units.o vmat.o
 diag2g.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
 diag2k.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
 diag3g.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
@@ -213,7 +215,7 @@
 domain_decom.o: sparse_matrices.o sys.o
 doping_uniform.o: alloc.o m_ntm.o mesh.o parallel.o precision.o sys.o
 dynamics.o: alloc.o atomlist.o files.o ioxv.o m_mpi_utils.o parallel.o
-dynamics.o: precision.o sys.o units.o
+dynamics.o: precision.o siesta_options.o sys.o units.o
 egandd.o: alloc.o denmat.o ener3.o globalise.o gradient.o on_subs.o onmod.o
 egandd.o: onmod.o precision.o sys.o
 eggbox.o: parallel.o precision.o
@@ -222,6 +224,10 @@
 elpa2.o: elpa1.o
 ener3.o: alloc.o globalise.o m_mpi_utils.o onmod.o precision.o
 ener3lomem.o: alloc.o globalise.o m_mpi_utils.o onmod.o precision.o
+evolg.o: alloc.o m_eo.o m_matswinvers.o m_steps.o parallel.o parallelsubs.o
+evolg.o: precision.o siesta_options.o sparse_matrices.o sys.o wavefunctions.o
+evolk.o: m_diagon_opt.o m_matswinvers.o parallel.o precision.o siesta_options.o
+evolk.o: sys.o wavefunctions.o
 extrapolateSpData2D.o: class_Data2D.o class_OrbitalDistribution.o
 extrapolateSpData2D.o: class_SpData2D.o class_Sparsity.o
 extrapolon.o: parallel.o precision.o sys.o
@@ -262,9 +268,9 @@
 initatom.o: atmparams.o atom.o atom_options.o basis_io.o basis_specs.o
 initatom.o: basis_types.o chemical.o electrostatic.o ldau_specs.o m_spin.o
 initatom.o: old_atmfuncs.o precision.o pseudopotential.o
-initparallel.o: alloc.o class_OrbitalDistribution.o domain_decom.o parallel.o
-initparallel.o: parallelsubs.o precision.o siesta_options.o sparse_matrices.o
-initparallel.o: spatial.o sys.o
+initparallel.o: alloc.o class_OrbitalDistribution.o domain_decom.o
+initparallel.o: m_poisson_fdf.o parallel.o parallelsubs.o precision.o
+initparallel.o: siesta_options.o sparse_matrices.o spatial.o sys.o
 io.o: sys.o
 iocg.o: files.o precision.o sys.o
 iodm.o: alloc.o files.o parallel.o parallelsubs.o precision.o sys.o
@@ -323,7 +329,9 @@
 m_efield.o: siesta_geom.o sys.o units.o
 m_energies.o: precision.o
 m_eo.o: precision.o
+m_evolve.o: alloc.o m_memory.o parallel.o parallelsubs.o precision.o sys.o
 m_exp_coord.o: files.o m_os.o parallel.o precision.o units.o
+m_fastconv_dims.o: precision.o
 m_filter.o: bessph.o precision.o radfft.o sorting.o sys.o
 m_fire.o: parallel.o precision.o
 m_fire_mixing.o: alloc.o m_fire_para.o m_mpi_utils.o parallel.o precision.o
@@ -345,10 +353,14 @@
 m_hs_matrix.o: alloc.o cellsubs.o geom_helper.o precision.o sys.o
 m_hsx.o: atm_types.o atmfuncs.o atomlist.o files.o parallel.o parallelsubs.o
 m_hsx.o: precision.o siesta_geom.o sys.o
+m_initwf.o: alloc.o densematrix.o fermid.o m_diagon_opt.o m_eo.o m_memory.o
+m_initwf.o: parallel.o parallelsubs.o precision.o sparse_matrices.o sys.o
+m_initwf.o: wavefunctions.o
 m_integrate.o: precision.o
 m_io.o: sys.o
 m_io_s.o: class_OrbitalDistribution.o class_SpData1D.o class_SpData2D.o
 m_io_s.o: class_Sparsity.o parallel.o precision.o
+m_iodipol.o: files.o m_dipol.o m_steps.o
 m_iodm.o: class_OrbitalDistribution.o class_SpData2D.o class_Sparsity.o
 m_iodm.o: m_io_s.o parallel.o
 m_iodm_old.o: alloc.o files.o m_matio.o m_mpi_utils.o parallel.o parallelsubs.o
@@ -356,9 +368,13 @@
 m_iorho.o: alloc.o parallel.o parallelsubs.o precision.o sys.o
 m_iostruct.o: alloc.o files.o m_mpi_utils.o parallel.o precision.o
 m_iostruct.o: siesta_geom.o sys.o units.o
+m_iotddft.o: files.o m_dipol.o m_steps.o parallel.o siesta_options.o
+m_iotddft.o: wavefunctions.o
 m_kinetic.o: precision.o
 m_mat_invert.o: intrinsic_missing.o m_pivot_array.o precision.o
+m_matdiag.o: precision.o
 m_matio.o: alloc.o parallel.o
+m_matswinvers.o: precision.o
 m_memory.o: memoryinfo.o parallel.o sys.o
 m_mesh_node.o: intrinsic_missing.o parallel.o precision.o
 m_mixing.o: class_Data1D.o class_Fstack_Data1D.o parallel.o precision.o
@@ -406,6 +422,15 @@
 m_planewavematrix.o: atomlist.o dhscf.o m_planewavematrixvar.o precision.o
 m_planewavematrix.o: siesta_geom.o sparse_matrices.o
 m_planewavematrixvar.o: precision.o
+m_plot_1d_arr3d.o: precision.o
+m_poisson_fdf.o: m_char.o precision.o
+m_poisson_obc_comm.o: atm_types.o atmfuncs.o m_mpi_stop.o parallel.o
+m_poisson_obc_comm.o: precision.o siesta_geom.o sys.o
+m_poisson_obc_fastconv.o: cellsubs.o m_mpi_stop.o m_poisson_obc_fastconv_r2c.o
+m_poisson_obc_fastconv.o: parallel.o precision.o sys.o
+m_poisson_obc_fastconv_r2c.o: cellsubs.o m_fastconv_dims.o m_plot_1d_arr3d.o
+m_poisson_obc_fastconv_r2c.o: m_poisson_obc_comm.o precision.o sys.o
+m_poisson_switch.o: m_poisson_obc_fastconv.o precision.o sys.o
 m_pulay.o: alloc.o atomlist.o m_mpi_utils.o m_spin.o m_svd.o parallel.o
 m_pulay.o: precision.o siesta_options.o sparse_matrices.o sys.o
 m_redist_spmatrix.o: alloc.o class_Distribution.o
@@ -421,7 +446,8 @@
 m_sparsity_handling.o: class_Data1D.o class_Data2D.o
 m_sparsity_handling.o: class_OrbitalDistribution.o class_SpData1D.o
 m_sparsity_handling.o: class_SpData2D.o class_Sparsity.o geom_helper.o
-m_sparsity_handling.o: m_interpolate.o m_region.o precision.o
+m_sparsity_handling.o: intrinsic_missing.o m_interpolate.o m_region.o
+m_sparsity_handling.o: precision.o
 m_spin.o: alloc.o parallel.o precision.o sys.o units.o
 m_stress.o: precision.o
 m_supercell.o: atom_graph.o class_OrbitalDistribution.o class_SpData2D.o
@@ -533,9 +559,10 @@
 m_ts_options.o: m_ts_hartree.o m_ts_method.o m_ts_mumps_init.o m_ts_tdir.o
 m_ts_options.o: m_ts_weight.o parallel.o precision.o siesta_options.o units.o
 m_ts_pivot.o: alloc.o class_OrbitalDistribution.o class_Sparsity.o
-m_ts_pivot.o: create_Sparsity_SC.o create_Sparsity_Union.o intrinsic_missing.o
-m_ts_pivot.o: m_char.o m_pivot.o m_region.o m_sparsity_handling.o m_ts_debug.o
-m_ts_pivot.o: m_ts_electype.o parallel.o precision.o
+m_ts_pivot.o: create_Sparsity_SC.o create_Sparsity_Union.o geom_helper.o
+m_ts_pivot.o: intrinsic_missing.o m_char.o m_pivot.o m_region.o
+m_ts_pivot.o: m_sparsity_handling.o m_ts_debug.o m_ts_electype.o parallel.o
+m_ts_pivot.o: precision.o
 m_ts_rgn2trimat.o: alloc.o class_OrbitalDistribution.o class_Sparsity.o
 m_ts_rgn2trimat.o: create_Sparsity_Union.o geom_helper.o m_region.o
 m_ts_rgn2trimat.o: m_ts_electype.o m_ts_method.o m_ts_tri_common.o parallel.o
@@ -678,8 +705,8 @@
 rdiag.o: alloc.o m_diagon_opt.o memoryinfo.o parallel.o precision.o sys.o
 rdiag_elpa.o: alloc.o elpa1.o elpa2.o parallel.o precision.o sys.o
 rdiag_mrrr.o: alloc.o parallel.o precision.o sys.o
-read_options.o: files.o m_charge_add.o m_hartree_add.o m_mixing_scf.o m_spin.o
-read_options.o: m_target_stress.o memoryinfo.o parallel.o precision.o
+read_options.o: files.o m_charge_add.o m_cite.o m_hartree_add.o m_mixing_scf.o
+read_options.o: m_spin.o m_target_stress.o memoryinfo.o parallel.o precision.o
 read_options.o: siesta_cml.o siesta_options.o sys.o units.o
 read_xc_info.o: parallel.o precision.o sys.o
 readsp.o: parallel.o precision.o sys.o
@@ -746,30 +773,32 @@
 siesta_end.o: meshdscf.o meshphi.o moremeshsubs.o parallel.o siesta_cml.o
 siesta_end.o: siesta_cmlsubs.o siesta_master.o siesta_options.o
 siesta_end.o: sparse_matrices.o timestamp.o
-siesta_forces.o: atomlist.o class_Fstack_Data1D.o class_SpData2D.o compute_dm.o
-siesta_forces.o: compute_energies.o compute_max_diff.o files.o
-siesta_forces.o: final_H_f_stress.o flook_siesta.o kpoint_grid.o
-siesta_forces.o: m_convergence.o m_energies.o m_forces.o m_iodm_old.o
-siesta_forces.o: m_mixing.o m_mixing_scf.o m_mpi_utils.o m_ncdf_siesta.o
-siesta_forces.o: m_pexsi.o m_pexsi_driver.o m_rhog.o m_spin.o m_steps.o
-siesta_forces.o: m_stress.o m_transiesta.o m_ts_charge.o m_ts_electype.o
-siesta_forces.o: m_ts_global_vars.o m_ts_method.o m_ts_options.o mixer.o
-siesta_forces.o: parallel.o post_scf_work.o precision.o save_density_matrix.o
-siesta_forces.o: scfconvergence_test.o setup_H0.o setup_hamiltonian.o
-siesta_forces.o: siesta_cml.o siesta_dicts.o siesta_geom.o siesta_master.o
-siesta_forces.o: siesta_options.o sparse_matrices.o state_analysis.o
-siesta_forces.o: state_init.o sys.o units.o write_subs.o
+siesta_forces.o: alloc.o atomlist.o class_Fstack_Data1D.o class_SpData2D.o
+siesta_forces.o: compute_dm.o compute_energies.o compute_max_diff.o
+siesta_forces.o: densematrix.o files.o final_H_f_stress.o flook_siesta.o
+siesta_forces.o: kpoint_grid.o m_convergence.o m_energies.o m_eo.o m_evolve.o
+siesta_forces.o: m_forces.o m_gamma.o m_initwf.o m_iodm.o m_iodm_old.o
+siesta_forces.o: m_iotddft.o m_mixing.o m_mixing_scf.o m_mpi_utils.o
+siesta_forces.o: m_ncdf_siesta.o m_pexsi.o m_pexsi_driver.o m_rhog.o m_spin.o
+siesta_forces.o: m_steps.o m_stress.o m_transiesta.o m_ts_charge.o
+siesta_forces.o: m_ts_electype.o m_ts_global_vars.o m_ts_method.o
+siesta_forces.o: m_ts_options.o mixer.o overfsm.o parallel.o post_scf_work.o
+siesta_forces.o: precision.o save_density_matrix.o scfconvergence_test.o
+siesta_forces.o: setup_H0.o setup_hamiltonian.o siesta_cml.o siesta_dicts.o
+siesta_forces.o: siesta_geom.o siesta_master.o siesta_options.o
+siesta_forces.o: sparse_matrices.o state_analysis.o state_init.o sys.o units.o
+siesta_forces.o: wavefunctions.o write_subs.o
 siesta_geom.o: precision.o
 siesta_init.o: alloc.o atomlist.o bands.o bsc_xcmod.o
 siesta_init.o: class_Fstack_Pair_Geometry_SpData2D.o densematrix.o files.o
 siesta_init.o: flook_siesta.o ioxv.o kpoint_grid.o kpoint_pdos.o ksvinit.o
 siesta_init.o: m_cite.o m_energies.o m_eo.o m_fixed.o m_forces.o m_gamma.o
-siesta_init.o: m_iostruct.o m_mpi_utils.o m_new_dm.o m_rmaxh.o m_spin.o
-siesta_init.o: m_steps.o m_supercell.o m_timer.o m_wallclock.o metaforce.o
-siesta_init.o: molecularmechanics.o object_debug.o parallel.o parallelsubs.o
-siesta_init.o: projected_DOS.o siesta_cmlsubs.o siesta_dicts.o siesta_geom.o
-siesta_init.o: siesta_options.o sparse_matrices.o struct_init.o sys.o timer.o
-siesta_init.o: timestamp.o ts_init.o writewave.o zmatrix.o
+siesta_init.o: m_iostruct.o m_mpi_utils.o m_new_dm.o m_poisson_fdf.o m_rmaxh.o
+siesta_init.o: m_spin.o m_steps.o m_supercell.o m_timer.o m_wallclock.o
+siesta_init.o: metaforce.o molecularmechanics.o object_debug.o parallel.o
+siesta_init.o: parallelsubs.o projected_DOS.o siesta_cmlsubs.o siesta_dicts.o
+siesta_init.o: siesta_geom.o siesta_options.o sparse_matrices.o struct_init.o
+siesta_init.o: sys.o timer.o timestamp.o ts_init.o writewave.o zmatrix.o
 siesta_master.o: iopipes.o iosockets.o precision.o sys.o
 siesta_move.o: atomlist.o broyden_optim.o cell_broyden_optim.o
 siesta_move.o: cell_fire_optim.o dynamics.o fire_optim.o flook_siesta.o ioxv.o
@@ -801,18 +830,17 @@
 state_init.o: proximity_check.o siesta_cml.o siesta_geom.o siesta_options.o
 state_init.o: sparse_matrices.o sys.o units.o write_subs.o zmatrix.o
 struct_init.o: alloc.o atmfuncs.o atomlist.o files.o ioxv.o m_exp_coord.o
-struct_init.o: m_iostruct.o m_mpi_utils.o m_steps.o parallel.o periodic_table.o
-struct_init.o: siesta_cml.o siesta_geom.o siesta_master.o siesta_options.o
-struct_init.o: units.o zmatrix.o
+struct_init.o: m_iostruct.o m_mpi_utils.o m_poisson_fdf.o m_poisson_obc_comm.o
+struct_init.o: m_steps.o parallel.o periodic_table.o siesta_cml.o siesta_geom.o
+struct_init.o: siesta_master.o siesta_options.o units.o zmatrix.o
 sys.o: parallel.o
 timer.o: extrae_eventllist.o extrae_module.o m_timer.o parallel.o sys.o
 timer.o: timer_tree.o
 timer_tree.o: m_walltime.o
 transition_rate.o: alloc.o fermid.o parallel.o parallelsubs.o precision.o sys.o
-ts_init.o: m_cite.o m_fixed.o m_os.o m_ts_cctype.o m_ts_charge.o
-ts_init.o: m_ts_electrode.o m_ts_electype.o m_ts_gf.o m_ts_global_vars.o
-ts_init.o: m_ts_kpoints.o m_ts_method.o m_ts_options.o parallel.o
-ts_init.o: siesta_options.o
+ts_init.o: m_fixed.o m_os.o m_ts_cctype.o m_ts_charge.o m_ts_electrode.o
+ts_init.o: m_ts_electype.o m_ts_gf.o m_ts_global_vars.o m_ts_kpoints.o
+ts_init.o: m_ts_method.o m_ts_options.o parallel.o siesta_options.o
 ts_show_regions.o: m_region.o m_ts_electype.o m_ts_method.o parallel.o
 ts_show_regions.o: precision.o units.o
 typecell.o: precision.o
@@ -824,6 +852,7 @@
 vmatsp.o: alloc.o atm_types.o atmfuncs.o atomlist.o listsc.o mesh.o meshdscf.o
 vmatsp.o: meshphi.o precision.o
 vmb.o: m_fixed.o parallel.o precision.o sys.o
+wavefunctions.o: parallel.o precision.o
 write_inp_wannier.o: alloc.o atmfuncs.o atomlist.o m_ntm.o m_orderbands.o
 write_inp_wannier.o: mneighb.o parallel.o parallelsubs.o precision.o
 write_inp_wannier.o: siesta2wannier90.o siesta_geom.o
@@ -941,6 +970,7 @@
 m_ksv.o: ksv.o
 m_ksvinit.o: ksvinit.o
 m_ldau.o: ldau.o
+m_linpack.o: linpack.o
 m_local_dos.o: local_DOS.o
 m_matel_registry.o: matel_registry.o
 m_minvec.o: minvec.o
@@ -985,6 +1015,8 @@
 m_vmat.o: vmat.o
 m_zm_broyden_optim.o: zm_broyden_optim.o
 m_zm_fire_optim.o: zm_fire_optim.o
+matdiagon.o: m_matdiag.o
+matswinversion.o: m_matswinvers.o
 mtprng.o: m_uuid.o
 neighbour.o: mneighb.o
 on_core.o: onmod.o

=== modified file 'Util/STM/ol-stm/Src/Makefile'
--- Util/STM/ol-stm/Src/Makefile	2016-08-30 14:02:33 +0000
+++ Util/STM/ol-stm/Src/Makefile	2016-11-22 08:16:45 +0000
@@ -185,7 +185,7 @@
 bsc_cellxc.o: alloc.o bsc_xcmod.o cellxc_mod.o mesh.o moremeshsubs.o parallel.o
 bsc_cellxc.o: parallelsubs.o precision.o sys.o
 bsc_xcmod.o: parallel.o precision.o sys.o
-cart2frac.o: sys.o
+cart2frac.o: precision.o sys.o
 cdiag.o: alloc.o m_diagon_opt.o memoryinfo.o parallel.o precision.o sys.o
 cell_broyden_optim.o: m_broyddj_nocomm.o parallel.o precision.o sys.o units.o
 cell_broyden_optim.o: zmatrix.o
@@ -197,6 +197,8 @@
 cgvc_zmatrix.o: alloc.o conjgr.o m_mpi_utils.o parallel.o precision.o sys.o
 cgvc_zmatrix.o: units.o zmatrix.o
 cgwf.o: egandd.o onmod.o parallel.o precision.o sys.o
+changebasis.o: alloc.o m_matdiag.o parallel.o parallelsubs.o precision.o
+changebasis.o: sparse_matrices.o sys.o wavefunctions.o
 chemical.o: parallel.o precision.o sys.o
 chempot.o: alloc.o mneighb.o parallel.o parallelsubs.o precision.o sys.o
 chkdim.o: sys.o
@@ -248,17 +250,17 @@
 denmat.o: alloc.o globalise.o m_mpi_utils.o onmod.o precision.o
 denmatlomem.o: alloc.o globalise.o onmod.o precision.o
 densematrix.o: alloc.o precision.o
-detover.o: alloc.o parallel.o parallelsubs.o precision.o
+detover.o: alloc.o linpack.o parallel.o parallelsubs.o precision.o
 dfscf.o: alloc.o atm_types.o atmfuncs.o atomlist.o listsc.o m_spin.o mesh.o
 dfscf.o: meshdscf.o meshphi.o parallel.o parallelsubs.o precision.o sys.o
 dhscf.o: alloc.o atmfuncs.o bsc_xcmod.o cellxc_mod.o delk.o dfscf.o
 dhscf.o: doping_uniform.o files.o forhar.o iogrid_netcdf.o m_charge_add.o
-dhscf.o: m_efield.o m_hartree_add.o m_iorho.o m_mesh_node.o m_ncdf_io.o
-dhscf.o: m_ncdf_siesta.o m_partial_charges.o m_rhog.o m_spin.o
-dhscf.o: m_ts_global_vars.o m_ts_hartree.o m_ts_options.o m_ts_tdir.o
-dhscf.o: m_ts_voltage.o mesh.o meshdscf.o meshsubs.o moremeshsubs.o parallel.o
-dhscf.o: parsing.o precision.o rhofft.o rhoofd.o siesta_options.o sys.o units.o
-dhscf.o: vmat.o
+dhscf.o: m_efield.o m_hartree_add.o m_iorho.o m_iotddft.o m_mesh_node.o
+dhscf.o: m_ncdf_io.o m_ncdf_siesta.o m_partial_charges.o m_poisson_fdf.o
+dhscf.o: m_poisson_switch.o m_rhog.o m_spin.o m_ts_global_vars.o m_ts_hartree.o
+dhscf.o: m_ts_options.o m_ts_tdir.o m_ts_voltage.o mesh.o meshdscf.o meshsubs.o
+dhscf.o: moremeshsubs.o parallel.o parsing.o precision.o rhofft.o rhoofd.o
+dhscf.o: siesta_options.o sys.o units.o vmat.o
 diag2g.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
 diag2k.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
 diag3g.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
@@ -286,7 +288,7 @@
 domain_decom.o: sparse_matrices.o sys.o
 doping_uniform.o: alloc.o m_ntm.o mesh.o parallel.o precision.o sys.o
 dynamics.o: alloc.o atomlist.o files.o ioxv.o m_mpi_utils.o parallel.o
-dynamics.o: precision.o sys.o units.o
+dynamics.o: precision.o siesta_options.o sys.o units.o
 egandd.o: alloc.o denmat.o ener3.o globalise.o gradient.o on_subs.o onmod.o
 egandd.o: onmod.o precision.o sys.o
 eggbox.o: parallel.o precision.o
@@ -295,6 +297,10 @@
 elpa2.o: elpa1.o
 ener3.o: alloc.o globalise.o m_mpi_utils.o onmod.o precision.o
 ener3lomem.o: alloc.o globalise.o m_mpi_utils.o onmod.o precision.o
+evolg.o: alloc.o m_eo.o m_matswinvers.o m_steps.o parallel.o parallelsubs.o
+evolg.o: precision.o siesta_options.o sparse_matrices.o sys.o wavefunctions.o
+evolk.o: m_diagon_opt.o m_matswinvers.o parallel.o precision.o siesta_options.o
+evolk.o: sys.o wavefunctions.o
 extrapolateSpData2D.o: class_Data2D.o class_OrbitalDistribution.o
 extrapolateSpData2D.o: class_SpData2D.o class_Sparsity.o
 extrapolon.o: parallel.o precision.o sys.o
@@ -335,9 +341,9 @@
 initatom.o: atmparams.o atom.o atom_options.o basis_io.o basis_specs.o
 initatom.o: basis_types.o chemical.o electrostatic.o ldau_specs.o m_spin.o
 initatom.o: old_atmfuncs.o precision.o pseudopotential.o
-initparallel.o: alloc.o class_OrbitalDistribution.o domain_decom.o parallel.o
-initparallel.o: parallelsubs.o precision.o siesta_options.o sparse_matrices.o
-initparallel.o: spatial.o sys.o
+initparallel.o: alloc.o class_OrbitalDistribution.o domain_decom.o
+initparallel.o: m_poisson_fdf.o parallel.o parallelsubs.o precision.o
+initparallel.o: siesta_options.o sparse_matrices.o spatial.o sys.o
 io.o: sys.o
 iocg.o: files.o precision.o sys.o
 iodm.o: alloc.o files.o parallel.o parallelsubs.o precision.o sys.o
@@ -396,7 +402,9 @@
 m_efield.o: siesta_geom.o sys.o units.o
 m_energies.o: precision.o
 m_eo.o: precision.o
+m_evolve.o: alloc.o m_memory.o parallel.o parallelsubs.o precision.o sys.o
 m_exp_coord.o: files.o m_os.o parallel.o precision.o units.o
+m_fastconv_dims.o: precision.o
 m_filter.o: bessph.o precision.o radfft.o sorting.o sys.o
 m_fire.o: parallel.o precision.o
 m_fire_mixing.o: alloc.o m_fire_para.o m_mpi_utils.o parallel.o precision.o
@@ -418,10 +426,14 @@
 m_hs_matrix.o: alloc.o cellsubs.o geom_helper.o precision.o sys.o
 m_hsx.o: atm_types.o atmfuncs.o atomlist.o files.o parallel.o parallelsubs.o
 m_hsx.o: precision.o siesta_geom.o sys.o
+m_initwf.o: alloc.o densematrix.o fermid.o m_diagon_opt.o m_eo.o m_memory.o
+m_initwf.o: parallel.o parallelsubs.o precision.o sparse_matrices.o sys.o
+m_initwf.o: wavefunctions.o
 m_integrate.o: precision.o
 m_io.o: sys.o
 m_io_s.o: class_OrbitalDistribution.o class_SpData1D.o class_SpData2D.o
 m_io_s.o: class_Sparsity.o parallel.o precision.o
+m_iodipol.o: files.o m_dipol.o m_steps.o
 m_iodm.o: class_OrbitalDistribution.o class_SpData2D.o class_Sparsity.o
 m_iodm.o: m_io_s.o parallel.o
 m_iodm_old.o: alloc.o files.o m_matio.o m_mpi_utils.o parallel.o parallelsubs.o
@@ -429,9 +441,13 @@
 m_iorho.o: alloc.o parallel.o parallelsubs.o precision.o sys.o
 m_iostruct.o: alloc.o files.o m_mpi_utils.o parallel.o precision.o
 m_iostruct.o: siesta_geom.o sys.o units.o
+m_iotddft.o: files.o m_dipol.o m_steps.o parallel.o siesta_options.o
+m_iotddft.o: wavefunctions.o
 m_kinetic.o: precision.o
 m_mat_invert.o: intrinsic_missing.o m_pivot_array.o precision.o
+m_matdiag.o: precision.o
 m_matio.o: alloc.o parallel.o
+m_matswinvers.o: precision.o
 m_memory.o: memoryinfo.o parallel.o sys.o
 m_mesh_node.o: intrinsic_missing.o parallel.o precision.o
 m_mixing.o: class_Data1D.o class_Fstack_Data1D.o parallel.o precision.o
@@ -479,6 +495,15 @@
 m_planewavematrix.o: atomlist.o dhscf.o m_planewavematrixvar.o precision.o
 m_planewavematrix.o: siesta_geom.o sparse_matrices.o
 m_planewavematrixvar.o: precision.o
+m_plot_1d_arr3d.o: precision.o
+m_poisson_fdf.o: m_char.o precision.o
+m_poisson_obc_comm.o: atm_types.o atmfuncs.o m_mpi_stop.o parallel.o
+m_poisson_obc_comm.o: precision.o siesta_geom.o sys.o
+m_poisson_obc_fastconv.o: cellsubs.o m_mpi_stop.o m_poisson_obc_fastconv_r2c.o
+m_poisson_obc_fastconv.o: parallel.o precision.o sys.o
+m_poisson_obc_fastconv_r2c.o: cellsubs.o m_fastconv_dims.o m_plot_1d_arr3d.o
+m_poisson_obc_fastconv_r2c.o: m_poisson_obc_comm.o precision.o sys.o
+m_poisson_switch.o: m_poisson_obc_fastconv.o precision.o sys.o
 m_pulay.o: alloc.o atomlist.o m_mpi_utils.o m_spin.o m_svd.o parallel.o
 m_pulay.o: precision.o siesta_options.o sparse_matrices.o sys.o
 m_redist_spmatrix.o: alloc.o class_Distribution.o
@@ -494,7 +519,8 @@
 m_sparsity_handling.o: class_Data1D.o class_Data2D.o
 m_sparsity_handling.o: class_OrbitalDistribution.o class_SpData1D.o
 m_sparsity_handling.o: class_SpData2D.o class_Sparsity.o geom_helper.o
-m_sparsity_handling.o: m_interpolate.o m_region.o precision.o
+m_sparsity_handling.o: intrinsic_missing.o m_interpolate.o m_region.o
+m_sparsity_handling.o: precision.o
 m_spin.o: alloc.o parallel.o precision.o sys.o units.o
 m_stress.o: precision.o
 m_supercell.o: atom_graph.o class_OrbitalDistribution.o class_SpData2D.o
@@ -606,9 +632,10 @@
 m_ts_options.o: m_ts_hartree.o m_ts_method.o m_ts_mumps_init.o m_ts_tdir.o
 m_ts_options.o: m_ts_weight.o parallel.o precision.o siesta_options.o units.o
 m_ts_pivot.o: alloc.o class_OrbitalDistribution.o class_Sparsity.o
-m_ts_pivot.o: create_Sparsity_SC.o create_Sparsity_Union.o intrinsic_missing.o
-m_ts_pivot.o: m_char.o m_pivot.o m_region.o m_sparsity_handling.o m_ts_debug.o
-m_ts_pivot.o: m_ts_electype.o parallel.o precision.o
+m_ts_pivot.o: create_Sparsity_SC.o create_Sparsity_Union.o geom_helper.o
+m_ts_pivot.o: intrinsic_missing.o m_char.o m_pivot.o m_region.o
+m_ts_pivot.o: m_sparsity_handling.o m_ts_debug.o m_ts_electype.o parallel.o
+m_ts_pivot.o: precision.o
 m_ts_rgn2trimat.o: alloc.o class_OrbitalDistribution.o class_Sparsity.o
 m_ts_rgn2trimat.o: create_Sparsity_Union.o geom_helper.o m_region.o
 m_ts_rgn2trimat.o: m_ts_electype.o m_ts_method.o m_ts_tri_common.o parallel.o
@@ -751,8 +778,8 @@
 rdiag.o: alloc.o m_diagon_opt.o memoryinfo.o parallel.o precision.o sys.o
 rdiag_elpa.o: alloc.o elpa1.o elpa2.o parallel.o precision.o sys.o
 rdiag_mrrr.o: alloc.o parallel.o precision.o sys.o
-read_options.o: files.o m_charge_add.o m_hartree_add.o m_mixing_scf.o m_spin.o
-read_options.o: m_target_stress.o memoryinfo.o parallel.o precision.o
+read_options.o: files.o m_charge_add.o m_cite.o m_hartree_add.o m_mixing_scf.o
+read_options.o: m_spin.o m_target_stress.o memoryinfo.o parallel.o precision.o
 read_options.o: siesta_cml.o siesta_options.o sys.o units.o
 read_xc_info.o: parallel.o precision.o sys.o
 readsp.o: parallel.o precision.o sys.o
@@ -819,30 +846,32 @@
 siesta_end.o: meshdscf.o meshphi.o moremeshsubs.o parallel.o siesta_cml.o
 siesta_end.o: siesta_cmlsubs.o siesta_master.o siesta_options.o
 siesta_end.o: sparse_matrices.o timestamp.o
-siesta_forces.o: atomlist.o class_Fstack_Data1D.o class_SpData2D.o compute_dm.o
-siesta_forces.o: compute_energies.o compute_max_diff.o files.o
-siesta_forces.o: final_H_f_stress.o flook_siesta.o kpoint_grid.o
-siesta_forces.o: m_convergence.o m_energies.o m_forces.o m_iodm_old.o
-siesta_forces.o: m_mixing.o m_mixing_scf.o m_mpi_utils.o m_ncdf_siesta.o
-siesta_forces.o: m_pexsi.o m_pexsi_driver.o m_rhog.o m_spin.o m_steps.o
-siesta_forces.o: m_stress.o m_transiesta.o m_ts_charge.o m_ts_electype.o
-siesta_forces.o: m_ts_global_vars.o m_ts_method.o m_ts_options.o mixer.o
-siesta_forces.o: parallel.o post_scf_work.o precision.o save_density_matrix.o
-siesta_forces.o: scfconvergence_test.o setup_H0.o setup_hamiltonian.o
-siesta_forces.o: siesta_cml.o siesta_dicts.o siesta_geom.o siesta_master.o
-siesta_forces.o: siesta_options.o sparse_matrices.o state_analysis.o
-siesta_forces.o: state_init.o sys.o units.o write_subs.o
+siesta_forces.o: alloc.o atomlist.o class_Fstack_Data1D.o class_SpData2D.o
+siesta_forces.o: compute_dm.o compute_energies.o compute_max_diff.o
+siesta_forces.o: densematrix.o files.o final_H_f_stress.o flook_siesta.o
+siesta_forces.o: kpoint_grid.o m_convergence.o m_energies.o m_eo.o m_evolve.o
+siesta_forces.o: m_forces.o m_gamma.o m_initwf.o m_iodm.o m_iodm_old.o
+siesta_forces.o: m_iotddft.o m_mixing.o m_mixing_scf.o m_mpi_utils.o
+siesta_forces.o: m_ncdf_siesta.o m_pexsi.o m_pexsi_driver.o m_rhog.o m_spin.o
+siesta_forces.o: m_steps.o m_stress.o m_transiesta.o m_ts_charge.o
+siesta_forces.o: m_ts_electype.o m_ts_global_vars.o m_ts_method.o
+siesta_forces.o: m_ts_options.o mixer.o overfsm.o parallel.o post_scf_work.o
+siesta_forces.o: precision.o save_density_matrix.o scfconvergence_test.o
+siesta_forces.o: setup_H0.o setup_hamiltonian.o siesta_cml.o siesta_dicts.o
+siesta_forces.o: siesta_geom.o siesta_master.o siesta_options.o
+siesta_forces.o: sparse_matrices.o state_analysis.o state_init.o sys.o units.o
+siesta_forces.o: wavefunctions.o write_subs.o
 siesta_geom.o: precision.o
 siesta_init.o: alloc.o atomlist.o bands.o bsc_xcmod.o
 siesta_init.o: class_Fstack_Pair_Geometry_SpData2D.o densematrix.o files.o
 siesta_init.o: flook_siesta.o ioxv.o kpoint_grid.o kpoint_pdos.o ksvinit.o
 siesta_init.o: m_cite.o m_energies.o m_eo.o m_fixed.o m_forces.o m_gamma.o
-siesta_init.o: m_iostruct.o m_mpi_utils.o m_new_dm.o m_rmaxh.o m_spin.o
-siesta_init.o: m_steps.o m_supercell.o m_timer.o m_wallclock.o metaforce.o
-siesta_init.o: molecularmechanics.o object_debug.o parallel.o parallelsubs.o
-siesta_init.o: projected_DOS.o siesta_cmlsubs.o siesta_dicts.o siesta_geom.o
-siesta_init.o: siesta_options.o sparse_matrices.o struct_init.o sys.o timer.o
-siesta_init.o: timestamp.o ts_init.o writewave.o zmatrix.o
+siesta_init.o: m_iostruct.o m_mpi_utils.o m_new_dm.o m_poisson_fdf.o m_rmaxh.o
+siesta_init.o: m_spin.o m_steps.o m_supercell.o m_timer.o m_wallclock.o
+siesta_init.o: metaforce.o molecularmechanics.o object_debug.o parallel.o
+siesta_init.o: parallelsubs.o projected_DOS.o siesta_cmlsubs.o siesta_dicts.o
+siesta_init.o: siesta_geom.o siesta_options.o sparse_matrices.o struct_init.o
+siesta_init.o: sys.o timer.o timestamp.o ts_init.o writewave.o zmatrix.o
 siesta_master.o: iopipes.o iosockets.o precision.o sys.o
 siesta_move.o: atomlist.o broyden_optim.o cell_broyden_optim.o
 siesta_move.o: cell_fire_optim.o dynamics.o fire_optim.o flook_siesta.o ioxv.o
@@ -874,18 +903,17 @@
 state_init.o: proximity_check.o siesta_cml.o siesta_geom.o siesta_options.o
 state_init.o: sparse_matrices.o sys.o units.o write_subs.o zmatrix.o
 struct_init.o: alloc.o atmfuncs.o atomlist.o files.o ioxv.o m_exp_coord.o
-struct_init.o: m_iostruct.o m_mpi_utils.o m_steps.o parallel.o periodic_table.o
-struct_init.o: siesta_cml.o siesta_geom.o siesta_master.o siesta_options.o
-struct_init.o: units.o zmatrix.o
+struct_init.o: m_iostruct.o m_mpi_utils.o m_poisson_fdf.o m_poisson_obc_comm.o
+struct_init.o: m_steps.o parallel.o periodic_table.o siesta_cml.o siesta_geom.o
+struct_init.o: siesta_master.o siesta_options.o units.o zmatrix.o
 sys.o: parallel.o
 timer.o: extrae_eventllist.o extrae_module.o m_timer.o parallel.o sys.o
 timer.o: timer_tree.o
 timer_tree.o: m_walltime.o
 transition_rate.o: alloc.o fermid.o parallel.o parallelsubs.o precision.o sys.o
-ts_init.o: m_cite.o m_fixed.o m_os.o m_ts_cctype.o m_ts_charge.o
-ts_init.o: m_ts_electrode.o m_ts_electype.o m_ts_gf.o m_ts_global_vars.o
-ts_init.o: m_ts_kpoints.o m_ts_method.o m_ts_options.o parallel.o
-ts_init.o: siesta_options.o
+ts_init.o: m_fixed.o m_os.o m_ts_cctype.o m_ts_charge.o m_ts_electrode.o
+ts_init.o: m_ts_electype.o m_ts_gf.o m_ts_global_vars.o m_ts_kpoints.o
+ts_init.o: m_ts_method.o m_ts_options.o parallel.o siesta_options.o
 ts_show_regions.o: m_region.o m_ts_electype.o m_ts_method.o parallel.o
 ts_show_regions.o: precision.o units.o
 typecell.o: precision.o
@@ -897,6 +925,7 @@
 vmatsp.o: alloc.o atm_types.o atmfuncs.o atomlist.o listsc.o mesh.o meshdscf.o
 vmatsp.o: meshphi.o precision.o
 vmb.o: m_fixed.o parallel.o precision.o sys.o
+wavefunctions.o: parallel.o precision.o
 write_inp_wannier.o: alloc.o atmfuncs.o atomlist.o m_ntm.o m_orderbands.o
 write_inp_wannier.o: mneighb.o parallel.o parallelsubs.o precision.o
 write_inp_wannier.o: siesta2wannier90.o siesta_geom.o
@@ -1020,6 +1049,7 @@
 m_ksv.o: ksv.o
 m_ksvinit.o: ksvinit.o
 m_ldau.o: ldau.o
+m_linpack.o: linpack.o
 m_local_dos.o: local_DOS.o
 m_matel_registry.o: matel_registry.o
 m_minvec.o: minvec.o
@@ -1064,6 +1094,8 @@
 m_vmat.o: vmat.o
 m_zm_broyden_optim.o: zm_broyden_optim.o
 m_zm_fire_optim.o: zm_fire_optim.o
+matdiagon.o: m_matdiag.o
+matswinversion.o: m_matswinvers.o
 mtprng.o: m_uuid.o
 neighbour.o: mneighb.o
 on_core.o: onmod.o

=== modified file 'Util/SpPivot/Makefile'
--- Util/SpPivot/Makefile	2016-08-30 14:02:33 +0000
+++ Util/SpPivot/Makefile	2016-11-22 08:16:45 +0000
@@ -116,7 +116,7 @@
 bsc_cellxc.o: alloc.o bsc_xcmod.o cellxc_mod.o mesh.o moremeshsubs.o parallel.o
 bsc_cellxc.o: parallelsubs.o precision.o sys.o
 bsc_xcmod.o: parallel.o precision.o sys.o
-cart2frac.o: sys.o
+cart2frac.o: precision.o sys.o
 cdiag.o: alloc.o m_diagon_opt.o memoryinfo.o parallel.o precision.o sys.o
 cell_broyden_optim.o: m_broyddj_nocomm.o parallel.o precision.o sys.o units.o
 cell_broyden_optim.o: zmatrix.o
@@ -128,6 +128,8 @@
 cgvc_zmatrix.o: alloc.o conjgr.o m_mpi_utils.o parallel.o precision.o sys.o
 cgvc_zmatrix.o: units.o zmatrix.o
 cgwf.o: egandd.o onmod.o parallel.o precision.o sys.o
+changebasis.o: alloc.o m_matdiag.o parallel.o parallelsubs.o precision.o
+changebasis.o: sparse_matrices.o sys.o wavefunctions.o
 chemical.o: parallel.o precision.o sys.o
 chempot.o: alloc.o mneighb.o parallel.o parallelsubs.o precision.o sys.o
 chkdim.o: sys.o
@@ -179,17 +181,17 @@
 denmat.o: alloc.o globalise.o m_mpi_utils.o onmod.o precision.o
 denmatlomem.o: alloc.o globalise.o onmod.o precision.o
 densematrix.o: alloc.o precision.o
-detover.o: alloc.o parallel.o parallelsubs.o precision.o
+detover.o: alloc.o linpack.o parallel.o parallelsubs.o precision.o
 dfscf.o: alloc.o atm_types.o atmfuncs.o atomlist.o listsc.o m_spin.o mesh.o
 dfscf.o: meshdscf.o meshphi.o parallel.o parallelsubs.o precision.o sys.o
 dhscf.o: alloc.o atmfuncs.o bsc_xcmod.o cellxc_mod.o delk.o dfscf.o
 dhscf.o: doping_uniform.o files.o forhar.o iogrid_netcdf.o m_charge_add.o
-dhscf.o: m_efield.o m_hartree_add.o m_iorho.o m_mesh_node.o m_ncdf_io.o
-dhscf.o: m_ncdf_siesta.o m_partial_charges.o m_rhog.o m_spin.o
-dhscf.o: m_ts_global_vars.o m_ts_hartree.o m_ts_options.o m_ts_tdir.o
-dhscf.o: m_ts_voltage.o mesh.o meshdscf.o meshsubs.o moremeshsubs.o parallel.o
-dhscf.o: parsing.o precision.o rhofft.o rhoofd.o siesta_options.o sys.o units.o
-dhscf.o: vmat.o
+dhscf.o: m_efield.o m_hartree_add.o m_iorho.o m_iotddft.o m_mesh_node.o
+dhscf.o: m_ncdf_io.o m_ncdf_siesta.o m_partial_charges.o m_poisson_fdf.o
+dhscf.o: m_poisson_switch.o m_rhog.o m_spin.o m_ts_global_vars.o m_ts_hartree.o
+dhscf.o: m_ts_options.o m_ts_tdir.o m_ts_voltage.o mesh.o meshdscf.o meshsubs.o
+dhscf.o: moremeshsubs.o parallel.o parsing.o precision.o rhofft.o rhoofd.o
+dhscf.o: siesta_options.o sys.o units.o vmat.o
 diag2g.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
 diag2k.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
 diag3g.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
@@ -217,7 +219,7 @@
 domain_decom.o: sparse_matrices.o sys.o
 doping_uniform.o: alloc.o m_ntm.o mesh.o parallel.o precision.o sys.o
 dynamics.o: alloc.o atomlist.o files.o ioxv.o m_mpi_utils.o parallel.o
-dynamics.o: precision.o sys.o units.o
+dynamics.o: precision.o siesta_options.o sys.o units.o
 egandd.o: alloc.o denmat.o ener3.o globalise.o gradient.o on_subs.o onmod.o
 egandd.o: onmod.o precision.o sys.o
 eggbox.o: parallel.o precision.o
@@ -226,6 +228,10 @@
 elpa2.o: elpa1.o
 ener3.o: alloc.o globalise.o m_mpi_utils.o onmod.o precision.o
 ener3lomem.o: alloc.o globalise.o m_mpi_utils.o onmod.o precision.o
+evolg.o: alloc.o m_eo.o m_matswinvers.o m_steps.o parallel.o parallelsubs.o
+evolg.o: precision.o siesta_options.o sparse_matrices.o sys.o wavefunctions.o
+evolk.o: m_diagon_opt.o m_matswinvers.o parallel.o precision.o siesta_options.o
+evolk.o: sys.o wavefunctions.o
 extrapolateSpData2D.o: class_Data2D.o class_OrbitalDistribution.o
 extrapolateSpData2D.o: class_SpData2D.o class_Sparsity.o
 extrapolon.o: parallel.o precision.o sys.o
@@ -266,9 +272,9 @@
 initatom.o: atmparams.o atom.o atom_options.o basis_io.o basis_specs.o
 initatom.o: basis_types.o chemical.o electrostatic.o ldau_specs.o m_spin.o
 initatom.o: old_atmfuncs.o precision.o pseudopotential.o
-initparallel.o: alloc.o class_OrbitalDistribution.o domain_decom.o parallel.o
-initparallel.o: parallelsubs.o precision.o siesta_options.o sparse_matrices.o
-initparallel.o: spatial.o sys.o
+initparallel.o: alloc.o class_OrbitalDistribution.o domain_decom.o
+initparallel.o: m_poisson_fdf.o parallel.o parallelsubs.o precision.o
+initparallel.o: siesta_options.o sparse_matrices.o spatial.o sys.o
 io.o: sys.o
 iocg.o: files.o precision.o sys.o
 iodm.o: alloc.o files.o parallel.o parallelsubs.o precision.o sys.o
@@ -327,7 +333,9 @@
 m_efield.o: siesta_geom.o sys.o units.o
 m_energies.o: precision.o
 m_eo.o: precision.o
+m_evolve.o: alloc.o m_memory.o parallel.o parallelsubs.o precision.o sys.o
 m_exp_coord.o: files.o m_os.o parallel.o precision.o units.o
+m_fastconv_dims.o: precision.o
 m_filter.o: bessph.o precision.o radfft.o sorting.o sys.o
 m_fire.o: parallel.o precision.o
 m_fire_mixing.o: alloc.o m_fire_para.o m_mpi_utils.o parallel.o precision.o
@@ -349,10 +357,14 @@
 m_hs_matrix.o: alloc.o cellsubs.o geom_helper.o precision.o sys.o
 m_hsx.o: atm_types.o atmfuncs.o atomlist.o files.o parallel.o parallelsubs.o
 m_hsx.o: precision.o siesta_geom.o sys.o
+m_initwf.o: alloc.o densematrix.o fermid.o m_diagon_opt.o m_eo.o m_memory.o
+m_initwf.o: parallel.o parallelsubs.o precision.o sparse_matrices.o sys.o
+m_initwf.o: wavefunctions.o
 m_integrate.o: precision.o
 m_io.o: sys.o
 m_io_s.o: class_OrbitalDistribution.o class_SpData1D.o class_SpData2D.o
 m_io_s.o: class_Sparsity.o parallel.o precision.o
+m_iodipol.o: files.o m_dipol.o m_steps.o
 m_iodm.o: class_OrbitalDistribution.o class_SpData2D.o class_Sparsity.o
 m_iodm.o: m_io_s.o parallel.o
 m_iodm_old.o: alloc.o files.o m_matio.o m_mpi_utils.o parallel.o parallelsubs.o
@@ -360,9 +372,13 @@
 m_iorho.o: alloc.o parallel.o parallelsubs.o precision.o sys.o
 m_iostruct.o: alloc.o files.o m_mpi_utils.o parallel.o precision.o
 m_iostruct.o: siesta_geom.o sys.o units.o
+m_iotddft.o: files.o m_dipol.o m_steps.o parallel.o siesta_options.o
+m_iotddft.o: wavefunctions.o
 m_kinetic.o: precision.o
 m_mat_invert.o: intrinsic_missing.o m_pivot_array.o precision.o
+m_matdiag.o: precision.o
 m_matio.o: alloc.o parallel.o
+m_matswinvers.o: precision.o
 m_memory.o: memoryinfo.o parallel.o sys.o
 m_mesh_node.o: intrinsic_missing.o parallel.o precision.o
 m_mixing.o: class_Data1D.o class_Fstack_Data1D.o parallel.o precision.o
@@ -410,6 +426,15 @@
 m_planewavematrix.o: atomlist.o dhscf.o m_planewavematrixvar.o precision.o
 m_planewavematrix.o: siesta_geom.o sparse_matrices.o
 m_planewavematrixvar.o: precision.o
+m_plot_1d_arr3d.o: precision.o
+m_poisson_fdf.o: m_char.o precision.o
+m_poisson_obc_comm.o: atm_types.o atmfuncs.o m_mpi_stop.o parallel.o
+m_poisson_obc_comm.o: precision.o siesta_geom.o sys.o
+m_poisson_obc_fastconv.o: cellsubs.o m_mpi_stop.o m_poisson_obc_fastconv_r2c.o
+m_poisson_obc_fastconv.o: parallel.o precision.o sys.o
+m_poisson_obc_fastconv_r2c.o: cellsubs.o m_fastconv_dims.o m_plot_1d_arr3d.o
+m_poisson_obc_fastconv_r2c.o: m_poisson_obc_comm.o precision.o sys.o
+m_poisson_switch.o: m_poisson_obc_fastconv.o precision.o sys.o
 m_pulay.o: alloc.o atomlist.o m_mpi_utils.o m_spin.o m_svd.o parallel.o
 m_pulay.o: precision.o siesta_options.o sparse_matrices.o sys.o
 m_redist_spmatrix.o: alloc.o class_Distribution.o
@@ -425,7 +450,8 @@
 m_sparsity_handling.o: class_Data1D.o class_Data2D.o
 m_sparsity_handling.o: class_OrbitalDistribution.o class_SpData1D.o
 m_sparsity_handling.o: class_SpData2D.o class_Sparsity.o geom_helper.o
-m_sparsity_handling.o: m_interpolate.o m_region.o precision.o
+m_sparsity_handling.o: intrinsic_missing.o m_interpolate.o m_region.o
+m_sparsity_handling.o: precision.o
 m_spin.o: alloc.o parallel.o precision.o sys.o units.o
 m_stress.o: precision.o
 m_supercell.o: atom_graph.o class_OrbitalDistribution.o class_SpData2D.o
@@ -537,9 +563,10 @@
 m_ts_options.o: m_ts_hartree.o m_ts_method.o m_ts_mumps_init.o m_ts_tdir.o
 m_ts_options.o: m_ts_weight.o parallel.o precision.o siesta_options.o units.o
 m_ts_pivot.o: alloc.o class_OrbitalDistribution.o class_Sparsity.o
-m_ts_pivot.o: create_Sparsity_SC.o create_Sparsity_Union.o intrinsic_missing.o
-m_ts_pivot.o: m_char.o m_pivot.o m_region.o m_sparsity_handling.o m_ts_debug.o
-m_ts_pivot.o: m_ts_electype.o parallel.o precision.o
+m_ts_pivot.o: create_Sparsity_SC.o create_Sparsity_Union.o geom_helper.o
+m_ts_pivot.o: intrinsic_missing.o m_char.o m_pivot.o m_region.o
+m_ts_pivot.o: m_sparsity_handling.o m_ts_debug.o m_ts_electype.o parallel.o
+m_ts_pivot.o: precision.o
 m_ts_rgn2trimat.o: alloc.o class_OrbitalDistribution.o class_Sparsity.o
 m_ts_rgn2trimat.o: create_Sparsity_Union.o geom_helper.o m_region.o
 m_ts_rgn2trimat.o: m_ts_electype.o m_ts_method.o m_ts_tri_common.o parallel.o
@@ -682,8 +709,8 @@
 rdiag.o: alloc.o m_diagon_opt.o memoryinfo.o parallel.o precision.o sys.o
 rdiag_elpa.o: alloc.o elpa1.o elpa2.o parallel.o precision.o sys.o
 rdiag_mrrr.o: alloc.o parallel.o precision.o sys.o
-read_options.o: files.o m_charge_add.o m_hartree_add.o m_mixing_scf.o m_spin.o
-read_options.o: m_target_stress.o memoryinfo.o parallel.o precision.o
+read_options.o: files.o m_charge_add.o m_cite.o m_hartree_add.o m_mixing_scf.o
+read_options.o: m_spin.o m_target_stress.o memoryinfo.o parallel.o precision.o
 read_options.o: siesta_cml.o siesta_options.o sys.o units.o
 read_xc_info.o: parallel.o precision.o sys.o
 readsp.o: parallel.o precision.o sys.o
@@ -750,30 +777,32 @@
 siesta_end.o: meshdscf.o meshphi.o moremeshsubs.o parallel.o siesta_cml.o
 siesta_end.o: siesta_cmlsubs.o siesta_master.o siesta_options.o
 siesta_end.o: sparse_matrices.o timestamp.o
-siesta_forces.o: atomlist.o class_Fstack_Data1D.o class_SpData2D.o compute_dm.o
-siesta_forces.o: compute_energies.o compute_max_diff.o files.o
-siesta_forces.o: final_H_f_stress.o flook_siesta.o kpoint_grid.o
-siesta_forces.o: m_convergence.o m_energies.o m_forces.o m_iodm_old.o
-siesta_forces.o: m_mixing.o m_mixing_scf.o m_mpi_utils.o m_ncdf_siesta.o
-siesta_forces.o: m_pexsi.o m_pexsi_driver.o m_rhog.o m_spin.o m_steps.o
-siesta_forces.o: m_stress.o m_transiesta.o m_ts_charge.o m_ts_electype.o
-siesta_forces.o: m_ts_global_vars.o m_ts_method.o m_ts_options.o mixer.o
-siesta_forces.o: parallel.o post_scf_work.o precision.o save_density_matrix.o
-siesta_forces.o: scfconvergence_test.o setup_H0.o setup_hamiltonian.o
-siesta_forces.o: siesta_cml.o siesta_dicts.o siesta_geom.o siesta_master.o
-siesta_forces.o: siesta_options.o sparse_matrices.o state_analysis.o
-siesta_forces.o: state_init.o sys.o units.o write_subs.o
+siesta_forces.o: alloc.o atomlist.o class_Fstack_Data1D.o class_SpData2D.o
+siesta_forces.o: compute_dm.o compute_energies.o compute_max_diff.o
+siesta_forces.o: densematrix.o files.o final_H_f_stress.o flook_siesta.o
+siesta_forces.o: kpoint_grid.o m_convergence.o m_energies.o m_eo.o m_evolve.o
+siesta_forces.o: m_forces.o m_gamma.o m_initwf.o m_iodm.o m_iodm_old.o
+siesta_forces.o: m_iotddft.o m_mixing.o m_mixing_scf.o m_mpi_utils.o
+siesta_forces.o: m_ncdf_siesta.o m_pexsi.o m_pexsi_driver.o m_rhog.o m_spin.o
+siesta_forces.o: m_steps.o m_stress.o m_transiesta.o m_ts_charge.o
+siesta_forces.o: m_ts_electype.o m_ts_global_vars.o m_ts_method.o
+siesta_forces.o: m_ts_options.o mixer.o overfsm.o parallel.o post_scf_work.o
+siesta_forces.o: precision.o save_density_matrix.o scfconvergence_test.o
+siesta_forces.o: setup_H0.o setup_hamiltonian.o siesta_cml.o siesta_dicts.o
+siesta_forces.o: siesta_geom.o siesta_master.o siesta_options.o
+siesta_forces.o: sparse_matrices.o state_analysis.o state_init.o sys.o units.o
+siesta_forces.o: wavefunctions.o write_subs.o
 siesta_geom.o: precision.o
 siesta_init.o: alloc.o atomlist.o bands.o bsc_xcmod.o
 siesta_init.o: class_Fstack_Pair_Geometry_SpData2D.o densematrix.o files.o
 siesta_init.o: flook_siesta.o ioxv.o kpoint_grid.o kpoint_pdos.o ksvinit.o
 siesta_init.o: m_cite.o m_energies.o m_eo.o m_fixed.o m_forces.o m_gamma.o
-siesta_init.o: m_iostruct.o m_mpi_utils.o m_new_dm.o m_rmaxh.o m_spin.o
-siesta_init.o: m_steps.o m_supercell.o m_timer.o m_wallclock.o metaforce.o
-siesta_init.o: molecularmechanics.o object_debug.o parallel.o parallelsubs.o
-siesta_init.o: projected_DOS.o siesta_cmlsubs.o siesta_dicts.o siesta_geom.o
-siesta_init.o: siesta_options.o sparse_matrices.o struct_init.o sys.o timer.o
-siesta_init.o: timestamp.o ts_init.o writewave.o zmatrix.o
+siesta_init.o: m_iostruct.o m_mpi_utils.o m_new_dm.o m_poisson_fdf.o m_rmaxh.o
+siesta_init.o: m_spin.o m_steps.o m_supercell.o m_timer.o m_wallclock.o
+siesta_init.o: metaforce.o molecularmechanics.o object_debug.o parallel.o
+siesta_init.o: parallelsubs.o projected_DOS.o siesta_cmlsubs.o siesta_dicts.o
+siesta_init.o: siesta_geom.o siesta_options.o sparse_matrices.o struct_init.o
+siesta_init.o: sys.o timer.o timestamp.o ts_init.o writewave.o zmatrix.o
 siesta_master.o: iopipes.o iosockets.o precision.o sys.o
 siesta_move.o: atomlist.o broyden_optim.o cell_broyden_optim.o
 siesta_move.o: cell_fire_optim.o dynamics.o fire_optim.o flook_siesta.o ioxv.o
@@ -805,18 +834,17 @@
 state_init.o: proximity_check.o siesta_cml.o siesta_geom.o siesta_options.o
 state_init.o: sparse_matrices.o sys.o units.o write_subs.o zmatrix.o
 struct_init.o: alloc.o atmfuncs.o atomlist.o files.o ioxv.o m_exp_coord.o
-struct_init.o: m_iostruct.o m_mpi_utils.o m_steps.o parallel.o periodic_table.o
-struct_init.o: siesta_cml.o siesta_geom.o siesta_master.o siesta_options.o
-struct_init.o: units.o zmatrix.o
+struct_init.o: m_iostruct.o m_mpi_utils.o m_poisson_fdf.o m_poisson_obc_comm.o
+struct_init.o: m_steps.o parallel.o periodic_table.o siesta_cml.o siesta_geom.o
+struct_init.o: siesta_master.o siesta_options.o units.o zmatrix.o
 sys.o: parallel.o
 timer.o: extrae_eventllist.o extrae_module.o m_timer.o parallel.o sys.o
 timer.o: timer_tree.o
 timer_tree.o: m_walltime.o
 transition_rate.o: alloc.o fermid.o parallel.o parallelsubs.o precision.o sys.o
-ts_init.o: m_cite.o m_fixed.o m_os.o m_ts_cctype.o m_ts_charge.o
-ts_init.o: m_ts_electrode.o m_ts_electype.o m_ts_gf.o m_ts_global_vars.o
-ts_init.o: m_ts_kpoints.o m_ts_method.o m_ts_options.o parallel.o
-ts_init.o: siesta_options.o
+ts_init.o: m_fixed.o m_os.o m_ts_cctype.o m_ts_charge.o m_ts_electrode.o
+ts_init.o: m_ts_electype.o m_ts_gf.o m_ts_global_vars.o m_ts_kpoints.o
+ts_init.o: m_ts_method.o m_ts_options.o parallel.o siesta_options.o
 ts_show_regions.o: m_region.o m_ts_electype.o m_ts_method.o parallel.o
 ts_show_regions.o: precision.o units.o
 typecell.o: precision.o
@@ -828,6 +856,7 @@
 vmatsp.o: alloc.o atm_types.o atmfuncs.o atomlist.o listsc.o mesh.o meshdscf.o
 vmatsp.o: meshphi.o precision.o
 vmb.o: m_fixed.o parallel.o precision.o sys.o
+wavefunctions.o: parallel.o precision.o
 write_inp_wannier.o: alloc.o atmfuncs.o atomlist.o m_ntm.o m_orderbands.o
 write_inp_wannier.o: mneighb.o parallel.o parallelsubs.o precision.o
 write_inp_wannier.o: siesta2wannier90.o siesta_geom.o
@@ -948,6 +977,7 @@
 m_ksv.o: ksv.o
 m_ksvinit.o: ksvinit.o
 m_ldau.o: ldau.o
+m_linpack.o: linpack.o
 m_local_dos.o: local_DOS.o
 m_matel_registry.o: matel_registry.o
 m_minvec.o: minvec.o
@@ -992,6 +1022,8 @@
 m_vmat.o: vmat.o
 m_zm_broyden_optim.o: zm_broyden_optim.o
 m_zm_fire_optim.o: zm_fire_optim.o
+matdiagon.o: m_matdiag.o
+matswinversion.o: m_matswinvers.o
 mtprng.o: m_uuid.o
 neighbour.o: mneighb.o
 on_core.o: onmod.o

=== modified file 'Util/TS/TBtrans/Makefile'
--- Util/TS/TBtrans/Makefile	2016-10-10 13:00:00 +0000
+++ Util/TS/TBtrans/Makefile	2016-11-22 08:16:45 +0000
@@ -364,6 +364,8 @@
 cgvc_zmatrix.o: alloc.o conjgr.o m_mpi_utils.o parallel.o precision.o sys.o
 cgvc_zmatrix.o: units.o zmatrix.o
 cgwf.o: egandd.o onmod.o parallel.o precision.o sys.o
+changebasis.o: alloc.o m_matdiag.o parallel.o parallelsubs.o precision.o
+changebasis.o: sparse_matrices.o sys.o wavefunctions.o
 chemical.o: parallel.o precision.o sys.o
 chempot.o: alloc.o mneighb.o parallel.o parallelsubs.o precision.o sys.o
 chkdim.o: sys.o
@@ -420,12 +422,12 @@
 dfscf.o: meshdscf.o meshphi.o parallel.o parallelsubs.o precision.o sys.o
 dhscf.o: alloc.o atmfuncs.o bsc_xcmod.o cellxc_mod.o delk.o dfscf.o
 dhscf.o: doping_uniform.o files.o forhar.o iogrid_netcdf.o m_charge_add.o
-dhscf.o: m_efield.o m_hartree_add.o m_iorho.o m_mesh_node.o m_ncdf_io.o
-dhscf.o: m_ncdf_siesta.o m_partial_charges.o m_rhog.o m_spin.o
-dhscf.o: m_ts_global_vars.o m_ts_hartree.o m_ts_options.o m_ts_tdir.o
-dhscf.o: m_ts_voltage.o mesh.o meshdscf.o meshsubs.o moremeshsubs.o parallel.o
-dhscf.o: parsing.o precision.o rhofft.o rhoofd.o siesta_options.o sys.o units.o
-dhscf.o: vmat.o
+dhscf.o: m_efield.o m_hartree_add.o m_iorho.o m_iotddft.o m_mesh_node.o
+dhscf.o: m_ncdf_io.o m_ncdf_siesta.o m_partial_charges.o m_poisson_fdf.o
+dhscf.o: m_poisson_switch.o m_rhog.o m_spin.o m_ts_global_vars.o m_ts_hartree.o
+dhscf.o: m_ts_options.o m_ts_tdir.o m_ts_voltage.o mesh.o meshdscf.o meshsubs.o
+dhscf.o: moremeshsubs.o parallel.o parsing.o precision.o rhofft.o rhoofd.o
+dhscf.o: siesta_options.o sys.o units.o vmat.o
 diag2g.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
 diag2k.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
 diag3g.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
@@ -453,7 +455,7 @@
 domain_decom.o: sparse_matrices.o sys.o
 doping_uniform.o: alloc.o m_ntm.o mesh.o parallel.o precision.o sys.o
 dynamics.o: alloc.o atomlist.o files.o ioxv.o m_mpi_utils.o parallel.o
-dynamics.o: precision.o sys.o units.o
+dynamics.o: precision.o siesta_options.o sys.o units.o
 egandd.o: alloc.o denmat.o ener3.o globalise.o gradient.o on_subs.o onmod.o
 egandd.o: onmod.o precision.o sys.o
 eggbox.o: parallel.o precision.o
@@ -462,6 +464,10 @@
 elpa2.o: elpa1.o
 ener3.o: alloc.o globalise.o m_mpi_utils.o onmod.o precision.o
 ener3lomem.o: alloc.o globalise.o m_mpi_utils.o onmod.o precision.o
+evolg.o: alloc.o m_eo.o m_matswinvers.o m_steps.o parallel.o parallelsubs.o
+evolg.o: precision.o siesta_options.o sparse_matrices.o sys.o wavefunctions.o
+evolk.o: m_diagon_opt.o m_matswinvers.o parallel.o precision.o siesta_options.o
+evolk.o: sys.o wavefunctions.o
 extrapolateSpData2D.o: class_Data2D.o class_OrbitalDistribution.o
 extrapolateSpData2D.o: class_SpData2D.o class_Sparsity.o
 extrapolon.o: parallel.o precision.o sys.o
@@ -502,9 +508,9 @@
 initatom.o: atmparams.o atom.o atom_options.o basis_io.o basis_specs.o
 initatom.o: basis_types.o chemical.o electrostatic.o ldau_specs.o m_spin.o
 initatom.o: old_atmfuncs.o precision.o pseudopotential.o
-initparallel.o: alloc.o class_OrbitalDistribution.o domain_decom.o parallel.o
-initparallel.o: parallelsubs.o precision.o siesta_options.o sparse_matrices.o
-initparallel.o: spatial.o sys.o
+initparallel.o: alloc.o class_OrbitalDistribution.o domain_decom.o
+initparallel.o: m_poisson_fdf.o parallel.o parallelsubs.o precision.o
+initparallel.o: siesta_options.o sparse_matrices.o spatial.o sys.o
 io.o: sys.o
 iocg.o: files.o precision.o sys.o
 iodm.o: alloc.o files.o parallel.o parallelsubs.o precision.o sys.o
@@ -563,7 +569,9 @@
 m_efield.o: siesta_geom.o sys.o units.o
 m_energies.o: precision.o
 m_eo.o: precision.o
+m_evolve.o: alloc.o m_memory.o parallel.o parallelsubs.o precision.o sys.o
 m_exp_coord.o: files.o m_os.o parallel.o precision.o units.o
+m_fastconv_dims.o: precision.o
 m_filter.o: bessph.o precision.o radfft.o sorting.o sys.o
 m_fire.o: parallel.o precision.o
 m_fire_mixing.o: alloc.o m_fire_para.o m_mpi_utils.o parallel.o precision.o
@@ -585,10 +593,14 @@
 m_hs_matrix.o: alloc.o cellsubs.o geom_helper.o precision.o sys.o
 m_hsx.o: atm_types.o atmfuncs.o atomlist.o files.o parallel.o parallelsubs.o
 m_hsx.o: precision.o siesta_geom.o sys.o
+m_initwf.o: alloc.o densematrix.o fermid.o m_diagon_opt.o m_eo.o m_memory.o
+m_initwf.o: parallel.o parallelsubs.o precision.o sparse_matrices.o sys.o
+m_initwf.o: wavefunctions.o
 m_integrate.o: precision.o
 m_io.o: sys.o
 m_io_s.o: class_OrbitalDistribution.o class_SpData1D.o class_SpData2D.o
 m_io_s.o: class_Sparsity.o parallel.o precision.o
+m_iodipol.o: files.o m_dipol.o m_steps.o
 m_iodm.o: class_OrbitalDistribution.o class_SpData2D.o class_Sparsity.o
 m_iodm.o: m_io_s.o parallel.o
 m_iodm_old.o: alloc.o files.o m_matio.o m_mpi_utils.o parallel.o parallelsubs.o
@@ -596,9 +608,13 @@
 m_iorho.o: alloc.o parallel.o parallelsubs.o precision.o sys.o
 m_iostruct.o: alloc.o files.o m_mpi_utils.o parallel.o precision.o
 m_iostruct.o: siesta_geom.o sys.o units.o
+m_iotddft.o: files.o m_dipol.o m_steps.o parallel.o siesta_options.o
+m_iotddft.o: wavefunctions.o
 m_kinetic.o: precision.o
 m_mat_invert.o: intrinsic_missing.o m_pivot_array.o precision.o
+m_matdiag.o: precision.o
 m_matio.o: alloc.o parallel.o
+m_matswinvers.o: precision.o
 m_memory.o: memoryinfo.o parallel.o sys.o
 m_mesh_node.o: intrinsic_missing.o parallel.o precision.o
 m_mixing.o: class_Data1D.o class_Fstack_Data1D.o parallel.o precision.o
@@ -646,6 +662,15 @@
 m_planewavematrix.o: atomlist.o dhscf.o m_planewavematrixvar.o precision.o
 m_planewavematrix.o: siesta_geom.o sparse_matrices.o
 m_planewavematrixvar.o: precision.o
+m_plot_1d_arr3d.o: precision.o
+m_poisson_fdf.o: m_char.o precision.o
+m_poisson_obc_comm.o: atm_types.o atmfuncs.o m_mpi_stop.o parallel.o
+m_poisson_obc_comm.o: precision.o siesta_geom.o sys.o
+m_poisson_obc_fastconv.o: cellsubs.o m_mpi_stop.o m_poisson_obc_fastconv_r2c.o
+m_poisson_obc_fastconv.o: parallel.o precision.o sys.o
+m_poisson_obc_fastconv_r2c.o: cellsubs.o m_fastconv_dims.o m_plot_1d_arr3d.o
+m_poisson_obc_fastconv_r2c.o: m_poisson_obc_comm.o precision.o sys.o
+m_poisson_switch.o: m_poisson_obc_fastconv.o precision.o sys.o
 m_pulay.o: alloc.o atomlist.o m_mpi_utils.o m_spin.o m_svd.o parallel.o
 m_pulay.o: precision.o siesta_options.o sparse_matrices.o sys.o
 m_redist_spmatrix.o: alloc.o class_Distribution.o
@@ -988,30 +1013,32 @@
 siesta_end.o: meshdscf.o meshphi.o moremeshsubs.o parallel.o siesta_cml.o
 siesta_end.o: siesta_cmlsubs.o siesta_master.o siesta_options.o
 siesta_end.o: sparse_matrices.o timestamp.o
-siesta_forces.o: atomlist.o class_Fstack_Data1D.o class_SpData2D.o compute_dm.o
-siesta_forces.o: compute_energies.o compute_max_diff.o files.o
-siesta_forces.o: final_H_f_stress.o flook_siesta.o kpoint_grid.o
-siesta_forces.o: m_convergence.o m_energies.o m_forces.o m_iodm_old.o
-siesta_forces.o: m_mixing.o m_mixing_scf.o m_mpi_utils.o m_ncdf_siesta.o
-siesta_forces.o: m_pexsi.o m_pexsi_driver.o m_rhog.o m_spin.o m_steps.o
-siesta_forces.o: m_stress.o m_transiesta.o m_ts_charge.o m_ts_electype.o
-siesta_forces.o: m_ts_global_vars.o m_ts_method.o m_ts_options.o mixer.o
-siesta_forces.o: parallel.o post_scf_work.o precision.o save_density_matrix.o
-siesta_forces.o: scfconvergence_test.o setup_H0.o setup_hamiltonian.o
-siesta_forces.o: siesta_cml.o siesta_dicts.o siesta_geom.o siesta_master.o
-siesta_forces.o: siesta_options.o sparse_matrices.o state_analysis.o
-siesta_forces.o: state_init.o sys.o units.o write_subs.o
+siesta_forces.o: alloc.o atomlist.o class_Fstack_Data1D.o class_SpData2D.o
+siesta_forces.o: compute_dm.o compute_energies.o compute_max_diff.o
+siesta_forces.o: densematrix.o files.o final_H_f_stress.o flook_siesta.o
+siesta_forces.o: kpoint_grid.o m_convergence.o m_energies.o m_eo.o m_evolve.o
+siesta_forces.o: m_forces.o m_gamma.o m_initwf.o m_iodm.o m_iodm_old.o
+siesta_forces.o: m_iotddft.o m_mixing.o m_mixing_scf.o m_mpi_utils.o
+siesta_forces.o: m_ncdf_siesta.o m_pexsi.o m_pexsi_driver.o m_rhog.o m_spin.o
+siesta_forces.o: m_steps.o m_stress.o m_transiesta.o m_ts_charge.o
+siesta_forces.o: m_ts_electype.o m_ts_global_vars.o m_ts_method.o
+siesta_forces.o: m_ts_options.o mixer.o overfsm.o parallel.o post_scf_work.o
+siesta_forces.o: precision.o save_density_matrix.o scfconvergence_test.o
+siesta_forces.o: setup_H0.o setup_hamiltonian.o siesta_cml.o siesta_dicts.o
+siesta_forces.o: siesta_geom.o siesta_master.o siesta_options.o
+siesta_forces.o: sparse_matrices.o state_analysis.o state_init.o sys.o units.o
+siesta_forces.o: wavefunctions.o write_subs.o
 siesta_geom.o: precision.o
 siesta_init.o: alloc.o atomlist.o bands.o bsc_xcmod.o
 siesta_init.o: class_Fstack_Pair_Geometry_SpData2D.o densematrix.o files.o
 siesta_init.o: flook_siesta.o ioxv.o kpoint_grid.o kpoint_pdos.o ksvinit.o
 siesta_init.o: m_cite.o m_energies.o m_eo.o m_fixed.o m_forces.o m_gamma.o
-siesta_init.o: m_iostruct.o m_mpi_utils.o m_new_dm.o m_rmaxh.o m_spin.o
-siesta_init.o: m_steps.o m_supercell.o m_timer.o m_wallclock.o metaforce.o
-siesta_init.o: molecularmechanics.o object_debug.o parallel.o parallelsubs.o
-siesta_init.o: projected_DOS.o siesta_cmlsubs.o siesta_dicts.o siesta_geom.o
-siesta_init.o: siesta_options.o sparse_matrices.o struct_init.o sys.o timer.o
-siesta_init.o: timestamp.o ts_init.o writewave.o zmatrix.o
+siesta_init.o: m_iostruct.o m_mpi_utils.o m_new_dm.o m_poisson_fdf.o m_rmaxh.o
+siesta_init.o: m_spin.o m_steps.o m_supercell.o m_timer.o m_wallclock.o
+siesta_init.o: metaforce.o molecularmechanics.o object_debug.o parallel.o
+siesta_init.o: parallelsubs.o projected_DOS.o siesta_cmlsubs.o siesta_dicts.o
+siesta_init.o: siesta_geom.o siesta_options.o sparse_matrices.o struct_init.o
+siesta_init.o: sys.o timer.o timestamp.o ts_init.o writewave.o zmatrix.o
 siesta_master.o: iopipes.o iosockets.o precision.o sys.o
 siesta_move.o: atomlist.o broyden_optim.o cell_broyden_optim.o
 siesta_move.o: cell_fire_optim.o dynamics.o fire_optim.o flook_siesta.o ioxv.o
@@ -1043,9 +1070,9 @@
 state_init.o: proximity_check.o siesta_cml.o siesta_geom.o siesta_options.o
 state_init.o: sparse_matrices.o sys.o units.o write_subs.o zmatrix.o
 struct_init.o: alloc.o atmfuncs.o atomlist.o files.o ioxv.o m_exp_coord.o
-struct_init.o: m_iostruct.o m_mpi_utils.o m_steps.o parallel.o periodic_table.o
-struct_init.o: siesta_cml.o siesta_geom.o siesta_master.o siesta_options.o
-struct_init.o: units.o zmatrix.o
+struct_init.o: m_iostruct.o m_mpi_utils.o m_poisson_fdf.o m_poisson_obc_comm.o
+struct_init.o: m_steps.o parallel.o periodic_table.o siesta_cml.o siesta_geom.o
+struct_init.o: siesta_master.o siesta_options.o units.o zmatrix.o
 sys.o: parallel.o
 timer.o: extrae_eventllist.o extrae_module.o m_timer.o parallel.o sys.o
 timer.o: timer_tree.o
@@ -1065,6 +1092,7 @@
 vmatsp.o: alloc.o atm_types.o atmfuncs.o atomlist.o listsc.o mesh.o meshdscf.o
 vmatsp.o: meshphi.o precision.o
 vmb.o: m_fixed.o parallel.o precision.o sys.o
+wavefunctions.o: parallel.o precision.o
 write_inp_wannier.o: alloc.o atmfuncs.o atomlist.o m_ntm.o m_orderbands.o
 write_inp_wannier.o: mneighb.o parallel.o parallelsubs.o precision.o
 write_inp_wannier.o: siesta2wannier90.o siesta_geom.o
@@ -1310,6 +1338,8 @@
 m_vmat.o: vmat.o
 m_zm_broyden_optim.o: zm_broyden_optim.o
 m_zm_fire_optim.o: zm_fire_optim.o
+matdiagon.o: m_matdiag.o
+matswinversion.o: m_matswinvers.o
 mtprng.o: m_uuid.o
 neighbour.o: mneighb.o
 on_core.o: onmod.o

=== modified file 'Util/TS/ts2ts/Makefile'
--- Util/TS/ts2ts/Makefile	2016-08-30 14:02:33 +0000
+++ Util/TS/ts2ts/Makefile	2016-11-22 08:16:45 +0000
@@ -111,7 +111,7 @@
 bsc_cellxc.o: alloc.o bsc_xcmod.o cellxc_mod.o mesh.o moremeshsubs.o parallel.o
 bsc_cellxc.o: parallelsubs.o precision.o sys.o
 bsc_xcmod.o: parallel.o precision.o sys.o
-cart2frac.o: sys.o
+cart2frac.o: precision.o sys.o
 cdiag.o: alloc.o m_diagon_opt.o memoryinfo.o parallel.o precision.o sys.o
 cell_broyden_optim.o: m_broyddj_nocomm.o parallel.o precision.o sys.o units.o
 cell_broyden_optim.o: zmatrix.o
@@ -123,6 +123,8 @@
 cgvc_zmatrix.o: alloc.o conjgr.o m_mpi_utils.o parallel.o precision.o sys.o
 cgvc_zmatrix.o: units.o zmatrix.o
 cgwf.o: egandd.o onmod.o parallel.o precision.o sys.o
+changebasis.o: alloc.o m_matdiag.o parallel.o parallelsubs.o precision.o
+changebasis.o: sparse_matrices.o sys.o wavefunctions.o
 chemical.o: parallel.o precision.o sys.o
 chempot.o: alloc.o mneighb.o parallel.o parallelsubs.o precision.o sys.o
 chkdim.o: sys.o
@@ -179,17 +181,17 @@
 denmat.o: alloc.o globalise.o m_mpi_utils.o onmod.o precision.o
 denmatlomem.o: alloc.o globalise.o onmod.o precision.o
 densematrix.o: alloc.o precision.o
-detover.o: alloc.o parallel.o parallelsubs.o precision.o
+detover.o: alloc.o linpack.o parallel.o parallelsubs.o precision.o
 dfscf.o: alloc.o atm_types.o atmfuncs.o atomlist.o listsc.o m_spin.o mesh.o
 dfscf.o: meshdscf.o meshphi.o parallel.o parallelsubs.o precision.o sys.o
 dhscf.o: alloc.o atmfuncs.o bsc_xcmod.o cellxc_mod.o delk.o dfscf.o
 dhscf.o: doping_uniform.o files.o forhar.o iogrid_netcdf.o m_charge_add.o
-dhscf.o: m_efield.o m_hartree_add.o m_iorho.o m_mesh_node.o m_ncdf_io.o
-dhscf.o: m_ncdf_siesta.o m_partial_charges.o m_rhog.o m_spin.o
-dhscf.o: m_ts_global_vars.o m_ts_hartree.o m_ts_options.o m_ts_tdir.o
-dhscf.o: m_ts_voltage.o mesh.o meshdscf.o meshsubs.o moremeshsubs.o parallel.o
-dhscf.o: parsing.o precision.o rhofft.o rhoofd.o siesta_options.o sys.o units.o
-dhscf.o: vmat.o
+dhscf.o: m_efield.o m_hartree_add.o m_iorho.o m_iotddft.o m_mesh_node.o
+dhscf.o: m_ncdf_io.o m_ncdf_siesta.o m_partial_charges.o m_poisson_fdf.o
+dhscf.o: m_poisson_switch.o m_rhog.o m_spin.o m_ts_global_vars.o m_ts_hartree.o
+dhscf.o: m_ts_options.o m_ts_tdir.o m_ts_voltage.o mesh.o meshdscf.o meshsubs.o
+dhscf.o: moremeshsubs.o parallel.o parsing.o precision.o rhofft.o rhoofd.o
+dhscf.o: siesta_options.o sys.o units.o vmat.o
 diag2g.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
 diag2k.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
 diag3g.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
@@ -217,7 +219,7 @@
 domain_decom.o: sparse_matrices.o sys.o
 doping_uniform.o: alloc.o m_ntm.o mesh.o parallel.o precision.o sys.o
 dynamics.o: alloc.o atomlist.o files.o ioxv.o m_mpi_utils.o parallel.o
-dynamics.o: precision.o sys.o units.o
+dynamics.o: precision.o siesta_options.o sys.o units.o
 egandd.o: alloc.o denmat.o ener3.o globalise.o gradient.o on_subs.o onmod.o
 egandd.o: onmod.o precision.o sys.o
 eggbox.o: parallel.o precision.o
@@ -226,6 +228,10 @@
 elpa2.o: elpa1.o
 ener3.o: alloc.o globalise.o m_mpi_utils.o onmod.o precision.o
 ener3lomem.o: alloc.o globalise.o m_mpi_utils.o onmod.o precision.o
+evolg.o: alloc.o m_eo.o m_matswinvers.o m_steps.o parallel.o parallelsubs.o
+evolg.o: precision.o siesta_options.o sparse_matrices.o sys.o wavefunctions.o
+evolk.o: m_diagon_opt.o m_matswinvers.o parallel.o precision.o siesta_options.o
+evolk.o: sys.o wavefunctions.o
 extrapolateSpData2D.o: class_Data2D.o class_OrbitalDistribution.o
 extrapolateSpData2D.o: class_SpData2D.o class_Sparsity.o
 extrapolon.o: parallel.o precision.o sys.o
@@ -266,9 +272,9 @@
 initatom.o: atmparams.o atom.o atom_options.o basis_io.o basis_specs.o
 initatom.o: basis_types.o chemical.o electrostatic.o ldau_specs.o m_spin.o
 initatom.o: old_atmfuncs.o precision.o pseudopotential.o
-initparallel.o: alloc.o class_OrbitalDistribution.o domain_decom.o parallel.o
-initparallel.o: parallelsubs.o precision.o siesta_options.o sparse_matrices.o
-initparallel.o: spatial.o sys.o
+initparallel.o: alloc.o class_OrbitalDistribution.o domain_decom.o
+initparallel.o: m_poisson_fdf.o parallel.o parallelsubs.o precision.o
+initparallel.o: siesta_options.o sparse_matrices.o spatial.o sys.o
 io.o: sys.o
 iocg.o: files.o precision.o sys.o
 iodm.o: alloc.o files.o parallel.o parallelsubs.o precision.o sys.o
@@ -327,7 +333,9 @@
 m_efield.o: siesta_geom.o sys.o units.o
 m_energies.o: precision.o
 m_eo.o: precision.o
+m_evolve.o: alloc.o m_memory.o parallel.o parallelsubs.o precision.o sys.o
 m_exp_coord.o: files.o m_os.o parallel.o precision.o units.o
+m_fastconv_dims.o: precision.o
 m_filter.o: bessph.o precision.o radfft.o sorting.o sys.o
 m_fire.o: parallel.o precision.o
 m_fire_mixing.o: alloc.o m_fire_para.o m_mpi_utils.o parallel.o precision.o
@@ -349,10 +357,14 @@
 m_hs_matrix.o: alloc.o cellsubs.o geom_helper.o precision.o sys.o
 m_hsx.o: atm_types.o atmfuncs.o atomlist.o files.o parallel.o parallelsubs.o
 m_hsx.o: precision.o siesta_geom.o sys.o
+m_initwf.o: alloc.o densematrix.o fermid.o m_diagon_opt.o m_eo.o m_memory.o
+m_initwf.o: parallel.o parallelsubs.o precision.o sparse_matrices.o sys.o
+m_initwf.o: wavefunctions.o
 m_integrate.o: precision.o
 m_io.o: sys.o
 m_io_s.o: class_OrbitalDistribution.o class_SpData1D.o class_SpData2D.o
 m_io_s.o: class_Sparsity.o parallel.o precision.o
+m_iodipol.o: files.o m_dipol.o m_steps.o
 m_iodm.o: class_OrbitalDistribution.o class_SpData2D.o class_Sparsity.o
 m_iodm.o: m_io_s.o parallel.o
 m_iodm_old.o: alloc.o files.o m_matio.o m_mpi_utils.o parallel.o parallelsubs.o
@@ -360,9 +372,13 @@
 m_iorho.o: alloc.o parallel.o parallelsubs.o precision.o sys.o
 m_iostruct.o: alloc.o files.o m_mpi_utils.o parallel.o precision.o
 m_iostruct.o: siesta_geom.o sys.o units.o
+m_iotddft.o: files.o m_dipol.o m_steps.o parallel.o siesta_options.o
+m_iotddft.o: wavefunctions.o
 m_kinetic.o: precision.o
 m_mat_invert.o: intrinsic_missing.o m_pivot_array.o precision.o
+m_matdiag.o: precision.o
 m_matio.o: alloc.o parallel.o
+m_matswinvers.o: precision.o
 m_memory.o: memoryinfo.o parallel.o sys.o
 m_mesh_node.o: intrinsic_missing.o parallel.o precision.o
 m_mixing.o: class_Data1D.o class_Fstack_Data1D.o parallel.o precision.o
@@ -410,6 +426,15 @@
 m_planewavematrix.o: atomlist.o dhscf.o m_planewavematrixvar.o precision.o
 m_planewavematrix.o: siesta_geom.o sparse_matrices.o
 m_planewavematrixvar.o: precision.o
+m_plot_1d_arr3d.o: precision.o
+m_poisson_fdf.o: m_char.o precision.o
+m_poisson_obc_comm.o: atm_types.o atmfuncs.o m_mpi_stop.o parallel.o
+m_poisson_obc_comm.o: precision.o siesta_geom.o sys.o
+m_poisson_obc_fastconv.o: cellsubs.o m_mpi_stop.o m_poisson_obc_fastconv_r2c.o
+m_poisson_obc_fastconv.o: parallel.o precision.o sys.o
+m_poisson_obc_fastconv_r2c.o: cellsubs.o m_fastconv_dims.o m_plot_1d_arr3d.o
+m_poisson_obc_fastconv_r2c.o: m_poisson_obc_comm.o precision.o sys.o
+m_poisson_switch.o: m_poisson_obc_fastconv.o precision.o sys.o
 m_pulay.o: alloc.o atomlist.o m_mpi_utils.o m_spin.o m_svd.o parallel.o
 m_pulay.o: precision.o siesta_options.o sparse_matrices.o sys.o
 m_redist_spmatrix.o: alloc.o class_Distribution.o
@@ -425,7 +450,8 @@
 m_sparsity_handling.o: class_Data1D.o class_Data2D.o
 m_sparsity_handling.o: class_OrbitalDistribution.o class_SpData1D.o
 m_sparsity_handling.o: class_SpData2D.o class_Sparsity.o geom_helper.o
-m_sparsity_handling.o: m_interpolate.o m_region.o precision.o
+m_sparsity_handling.o: intrinsic_missing.o m_interpolate.o m_region.o
+m_sparsity_handling.o: precision.o
 m_spin.o: alloc.o parallel.o precision.o sys.o units.o
 m_stress.o: precision.o
 m_supercell.o: atom_graph.o class_OrbitalDistribution.o class_SpData2D.o
@@ -537,9 +563,10 @@
 m_ts_options.o: m_ts_hartree.o m_ts_method.o m_ts_mumps_init.o m_ts_tdir.o
 m_ts_options.o: m_ts_weight.o parallel.o precision.o siesta_options.o units.o
 m_ts_pivot.o: alloc.o class_OrbitalDistribution.o class_Sparsity.o
-m_ts_pivot.o: create_Sparsity_SC.o create_Sparsity_Union.o intrinsic_missing.o
-m_ts_pivot.o: m_char.o m_pivot.o m_region.o m_sparsity_handling.o m_ts_debug.o
-m_ts_pivot.o: m_ts_electype.o parallel.o precision.o
+m_ts_pivot.o: create_Sparsity_SC.o create_Sparsity_Union.o geom_helper.o
+m_ts_pivot.o: intrinsic_missing.o m_char.o m_pivot.o m_region.o
+m_ts_pivot.o: m_sparsity_handling.o m_ts_debug.o m_ts_electype.o parallel.o
+m_ts_pivot.o: precision.o
 m_ts_rgn2trimat.o: alloc.o class_OrbitalDistribution.o class_Sparsity.o
 m_ts_rgn2trimat.o: create_Sparsity_Union.o geom_helper.o m_region.o
 m_ts_rgn2trimat.o: m_ts_electype.o m_ts_method.o m_ts_tri_common.o parallel.o
@@ -682,8 +709,8 @@
 rdiag.o: alloc.o m_diagon_opt.o memoryinfo.o parallel.o precision.o sys.o
 rdiag_elpa.o: alloc.o elpa1.o elpa2.o parallel.o precision.o sys.o
 rdiag_mrrr.o: alloc.o parallel.o precision.o sys.o
-read_options.o: files.o m_charge_add.o m_hartree_add.o m_mixing_scf.o m_spin.o
-read_options.o: m_target_stress.o memoryinfo.o parallel.o precision.o
+read_options.o: files.o m_charge_add.o m_cite.o m_hartree_add.o m_mixing_scf.o
+read_options.o: m_spin.o m_target_stress.o memoryinfo.o parallel.o precision.o
 read_options.o: siesta_cml.o siesta_options.o sys.o units.o
 read_xc_info.o: parallel.o precision.o sys.o
 readsp.o: parallel.o precision.o sys.o
@@ -750,30 +777,32 @@
 siesta_end.o: meshdscf.o meshphi.o moremeshsubs.o parallel.o siesta_cml.o
 siesta_end.o: siesta_cmlsubs.o siesta_master.o siesta_options.o
 siesta_end.o: sparse_matrices.o timestamp.o
-siesta_forces.o: atomlist.o class_Fstack_Data1D.o class_SpData2D.o compute_dm.o
-siesta_forces.o: compute_energies.o compute_max_diff.o files.o
-siesta_forces.o: final_H_f_stress.o flook_siesta.o kpoint_grid.o
-siesta_forces.o: m_convergence.o m_energies.o m_forces.o m_iodm_old.o
-siesta_forces.o: m_mixing.o m_mixing_scf.o m_mpi_utils.o m_ncdf_siesta.o
-siesta_forces.o: m_pexsi.o m_pexsi_driver.o m_rhog.o m_spin.o m_steps.o
-siesta_forces.o: m_stress.o m_transiesta.o m_ts_charge.o m_ts_electype.o
-siesta_forces.o: m_ts_global_vars.o m_ts_method.o m_ts_options.o mixer.o
-siesta_forces.o: parallel.o post_scf_work.o precision.o save_density_matrix.o
-siesta_forces.o: scfconvergence_test.o setup_H0.o setup_hamiltonian.o
-siesta_forces.o: siesta_cml.o siesta_dicts.o siesta_geom.o siesta_master.o
-siesta_forces.o: siesta_options.o sparse_matrices.o state_analysis.o
-siesta_forces.o: state_init.o sys.o units.o write_subs.o
+siesta_forces.o: alloc.o atomlist.o class_Fstack_Data1D.o class_SpData2D.o
+siesta_forces.o: compute_dm.o compute_energies.o compute_max_diff.o
+siesta_forces.o: densematrix.o files.o final_H_f_stress.o flook_siesta.o
+siesta_forces.o: kpoint_grid.o m_convergence.o m_energies.o m_eo.o m_evolve.o
+siesta_forces.o: m_forces.o m_gamma.o m_initwf.o m_iodm.o m_iodm_old.o
+siesta_forces.o: m_iotddft.o m_mixing.o m_mixing_scf.o m_mpi_utils.o
+siesta_forces.o: m_ncdf_siesta.o m_pexsi.o m_pexsi_driver.o m_rhog.o m_spin.o
+siesta_forces.o: m_steps.o m_stress.o m_transiesta.o m_ts_charge.o
+siesta_forces.o: m_ts_electype.o m_ts_global_vars.o m_ts_method.o
+siesta_forces.o: m_ts_options.o mixer.o overfsm.o parallel.o post_scf_work.o
+siesta_forces.o: precision.o save_density_matrix.o scfconvergence_test.o
+siesta_forces.o: setup_H0.o setup_hamiltonian.o siesta_cml.o siesta_dicts.o
+siesta_forces.o: siesta_geom.o siesta_master.o siesta_options.o
+siesta_forces.o: sparse_matrices.o state_analysis.o state_init.o sys.o units.o
+siesta_forces.o: wavefunctions.o write_subs.o
 siesta_geom.o: precision.o
 siesta_init.o: alloc.o atomlist.o bands.o bsc_xcmod.o
 siesta_init.o: class_Fstack_Pair_Geometry_SpData2D.o densematrix.o files.o
 siesta_init.o: flook_siesta.o ioxv.o kpoint_grid.o kpoint_pdos.o ksvinit.o
 siesta_init.o: m_cite.o m_energies.o m_eo.o m_fixed.o m_forces.o m_gamma.o
-siesta_init.o: m_iostruct.o m_mpi_utils.o m_new_dm.o m_rmaxh.o m_spin.o
-siesta_init.o: m_steps.o m_supercell.o m_timer.o m_wallclock.o metaforce.o
-siesta_init.o: molecularmechanics.o object_debug.o parallel.o parallelsubs.o
-siesta_init.o: projected_DOS.o siesta_cmlsubs.o siesta_dicts.o siesta_geom.o
-siesta_init.o: siesta_options.o sparse_matrices.o struct_init.o sys.o timer.o
-siesta_init.o: timestamp.o ts_init.o writewave.o zmatrix.o
+siesta_init.o: m_iostruct.o m_mpi_utils.o m_new_dm.o m_poisson_fdf.o m_rmaxh.o
+siesta_init.o: m_spin.o m_steps.o m_supercell.o m_timer.o m_wallclock.o
+siesta_init.o: metaforce.o molecularmechanics.o object_debug.o parallel.o
+siesta_init.o: parallelsubs.o projected_DOS.o siesta_cmlsubs.o siesta_dicts.o
+siesta_init.o: siesta_geom.o siesta_options.o sparse_matrices.o struct_init.o
+siesta_init.o: sys.o timer.o timestamp.o ts_init.o writewave.o zmatrix.o
 siesta_master.o: iopipes.o iosockets.o precision.o sys.o
 siesta_move.o: atomlist.o broyden_optim.o cell_broyden_optim.o
 siesta_move.o: cell_fire_optim.o dynamics.o fire_optim.o flook_siesta.o ioxv.o
@@ -805,18 +834,17 @@
 state_init.o: proximity_check.o siesta_cml.o siesta_geom.o siesta_options.o
 state_init.o: sparse_matrices.o sys.o units.o write_subs.o zmatrix.o
 struct_init.o: alloc.o atmfuncs.o atomlist.o files.o ioxv.o m_exp_coord.o
-struct_init.o: m_iostruct.o m_mpi_utils.o m_steps.o parallel.o periodic_table.o
-struct_init.o: siesta_cml.o siesta_geom.o siesta_master.o siesta_options.o
-struct_init.o: units.o zmatrix.o
+struct_init.o: m_iostruct.o m_mpi_utils.o m_poisson_fdf.o m_poisson_obc_comm.o
+struct_init.o: m_steps.o parallel.o periodic_table.o siesta_cml.o siesta_geom.o
+struct_init.o: siesta_master.o siesta_options.o units.o zmatrix.o
 sys.o: parallel.o
 timer.o: extrae_eventllist.o extrae_module.o m_timer.o parallel.o sys.o
 timer.o: timer_tree.o
 timer_tree.o: m_walltime.o
 transition_rate.o: alloc.o fermid.o parallel.o parallelsubs.o precision.o sys.o
-ts_init.o: m_cite.o m_fixed.o m_os.o m_ts_cctype.o m_ts_charge.o
-ts_init.o: m_ts_electrode.o m_ts_electype.o m_ts_gf.o m_ts_global_vars.o
-ts_init.o: m_ts_kpoints.o m_ts_method.o m_ts_options.o parallel.o
-ts_init.o: siesta_options.o
+ts_init.o: m_fixed.o m_os.o m_ts_cctype.o m_ts_charge.o m_ts_electrode.o
+ts_init.o: m_ts_electype.o m_ts_gf.o m_ts_global_vars.o m_ts_kpoints.o
+ts_init.o: m_ts_method.o m_ts_options.o parallel.o siesta_options.o
 ts_show_regions.o: m_region.o m_ts_electype.o m_ts_method.o parallel.o
 ts_show_regions.o: precision.o units.o
 typecell.o: precision.o
@@ -828,6 +856,7 @@
 vmatsp.o: alloc.o atm_types.o atmfuncs.o atomlist.o listsc.o mesh.o meshdscf.o
 vmatsp.o: meshphi.o precision.o
 vmb.o: m_fixed.o parallel.o precision.o sys.o
+wavefunctions.o: parallel.o precision.o
 write_inp_wannier.o: alloc.o atmfuncs.o atomlist.o m_ntm.o m_orderbands.o
 write_inp_wannier.o: mneighb.o parallel.o parallelsubs.o precision.o
 write_inp_wannier.o: siesta2wannier90.o siesta_geom.o
@@ -945,6 +974,7 @@
 m_ksv.o: ksv.o
 m_ksvinit.o: ksvinit.o
 m_ldau.o: ldau.o
+m_linpack.o: linpack.o
 m_local_dos.o: local_DOS.o
 m_matel_registry.o: matel_registry.o
 m_minvec.o: minvec.o
@@ -989,6 +1019,8 @@
 m_vmat.o: vmat.o
 m_zm_broyden_optim.o: zm_broyden_optim.o
 m_zm_fire_optim.o: zm_fire_optim.o
+matdiagon.o: m_matdiag.o
+matswinversion.o: m_matswinvers.o
 mtprng.o: m_uuid.o
 neighbour.o: mneighb.o
 on_core.o: onmod.o

=== modified file 'Util/TS/tshs2tshs/Makefile'
--- Util/TS/tshs2tshs/Makefile	2016-08-30 14:02:33 +0000
+++ Util/TS/tshs2tshs/Makefile	2016-11-22 08:16:45 +0000
@@ -169,7 +169,7 @@
 bsc_cellxc.o: alloc.o bsc_xcmod.o cellxc_mod.o mesh.o moremeshsubs.o parallel.o
 bsc_cellxc.o: parallelsubs.o precision.o sys.o
 bsc_xcmod.o: parallel.o precision.o sys.o
-cart2frac.o: sys.o
+cart2frac.o: precision.o sys.o
 cdiag.o: alloc.o m_diagon_opt.o memoryinfo.o parallel.o precision.o sys.o
 cell_broyden_optim.o: m_broyddj_nocomm.o parallel.o precision.o sys.o units.o
 cell_broyden_optim.o: zmatrix.o
@@ -181,6 +181,8 @@
 cgvc_zmatrix.o: alloc.o conjgr.o m_mpi_utils.o parallel.o precision.o sys.o
 cgvc_zmatrix.o: units.o zmatrix.o
 cgwf.o: egandd.o onmod.o parallel.o precision.o sys.o
+changebasis.o: alloc.o m_matdiag.o parallel.o parallelsubs.o precision.o
+changebasis.o: sparse_matrices.o sys.o wavefunctions.o
 chemical.o: parallel.o precision.o sys.o
 chempot.o: alloc.o mneighb.o parallel.o parallelsubs.o precision.o sys.o
 chkdim.o: sys.o
@@ -237,17 +239,17 @@
 denmat.o: alloc.o globalise.o m_mpi_utils.o onmod.o precision.o
 denmatlomem.o: alloc.o globalise.o onmod.o precision.o
 densematrix.o: alloc.o precision.o
-detover.o: alloc.o parallel.o parallelsubs.o precision.o
+detover.o: alloc.o linpack.o parallel.o parallelsubs.o precision.o
 dfscf.o: alloc.o atm_types.o atmfuncs.o atomlist.o listsc.o m_spin.o mesh.o
 dfscf.o: meshdscf.o meshphi.o parallel.o parallelsubs.o precision.o sys.o
 dhscf.o: alloc.o atmfuncs.o bsc_xcmod.o cellxc_mod.o delk.o dfscf.o
 dhscf.o: doping_uniform.o files.o forhar.o iogrid_netcdf.o m_charge_add.o
-dhscf.o: m_efield.o m_hartree_add.o m_iorho.o m_mesh_node.o m_ncdf_io.o
-dhscf.o: m_ncdf_siesta.o m_partial_charges.o m_rhog.o m_spin.o
-dhscf.o: m_ts_global_vars.o m_ts_hartree.o m_ts_options.o m_ts_tdir.o
-dhscf.o: m_ts_voltage.o mesh.o meshdscf.o meshsubs.o moremeshsubs.o parallel.o
-dhscf.o: parsing.o precision.o rhofft.o rhoofd.o siesta_options.o sys.o units.o
-dhscf.o: vmat.o
+dhscf.o: m_efield.o m_hartree_add.o m_iorho.o m_iotddft.o m_mesh_node.o
+dhscf.o: m_ncdf_io.o m_ncdf_siesta.o m_partial_charges.o m_poisson_fdf.o
+dhscf.o: m_poisson_switch.o m_rhog.o m_spin.o m_ts_global_vars.o m_ts_hartree.o
+dhscf.o: m_ts_options.o m_ts_tdir.o m_ts_voltage.o mesh.o meshdscf.o meshsubs.o
+dhscf.o: moremeshsubs.o parallel.o parsing.o precision.o rhofft.o rhoofd.o
+dhscf.o: siesta_options.o sys.o units.o vmat.o
 diag2g.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
 diag2k.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
 diag3g.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
@@ -275,7 +277,7 @@
 domain_decom.o: sparse_matrices.o sys.o
 doping_uniform.o: alloc.o m_ntm.o mesh.o parallel.o precision.o sys.o
 dynamics.o: alloc.o atomlist.o files.o ioxv.o m_mpi_utils.o parallel.o
-dynamics.o: precision.o sys.o units.o
+dynamics.o: precision.o siesta_options.o sys.o units.o
 egandd.o: alloc.o denmat.o ener3.o globalise.o gradient.o on_subs.o onmod.o
 egandd.o: onmod.o precision.o sys.o
 eggbox.o: parallel.o precision.o
@@ -284,6 +286,10 @@
 elpa2.o: elpa1.o
 ener3.o: alloc.o globalise.o m_mpi_utils.o onmod.o precision.o
 ener3lomem.o: alloc.o globalise.o m_mpi_utils.o onmod.o precision.o
+evolg.o: alloc.o m_eo.o m_matswinvers.o m_steps.o parallel.o parallelsubs.o
+evolg.o: precision.o siesta_options.o sparse_matrices.o sys.o wavefunctions.o
+evolk.o: m_diagon_opt.o m_matswinvers.o parallel.o precision.o siesta_options.o
+evolk.o: sys.o wavefunctions.o
 extrapolateSpData2D.o: class_Data2D.o class_OrbitalDistribution.o
 extrapolateSpData2D.o: class_SpData2D.o class_Sparsity.o
 extrapolon.o: parallel.o precision.o sys.o
@@ -324,9 +330,9 @@
 initatom.o: atmparams.o atom.o atom_options.o basis_io.o basis_specs.o
 initatom.o: basis_types.o chemical.o electrostatic.o ldau_specs.o m_spin.o
 initatom.o: old_atmfuncs.o precision.o pseudopotential.o
-initparallel.o: alloc.o class_OrbitalDistribution.o domain_decom.o parallel.o
-initparallel.o: parallelsubs.o precision.o siesta_options.o sparse_matrices.o
-initparallel.o: spatial.o sys.o
+initparallel.o: alloc.o class_OrbitalDistribution.o domain_decom.o
+initparallel.o: m_poisson_fdf.o parallel.o parallelsubs.o precision.o
+initparallel.o: siesta_options.o sparse_matrices.o spatial.o sys.o
 io.o: sys.o
 iocg.o: files.o precision.o sys.o
 iodm.o: alloc.o files.o parallel.o parallelsubs.o precision.o sys.o
@@ -385,7 +391,9 @@
 m_efield.o: siesta_geom.o sys.o units.o
 m_energies.o: precision.o
 m_eo.o: precision.o
+m_evolve.o: alloc.o m_memory.o parallel.o parallelsubs.o precision.o sys.o
 m_exp_coord.o: files.o m_os.o parallel.o precision.o units.o
+m_fastconv_dims.o: precision.o
 m_filter.o: bessph.o precision.o radfft.o sorting.o sys.o
 m_fire.o: parallel.o precision.o
 m_fire_mixing.o: alloc.o m_fire_para.o m_mpi_utils.o parallel.o precision.o
@@ -407,10 +415,14 @@
 m_hs_matrix.o: alloc.o cellsubs.o geom_helper.o precision.o sys.o
 m_hsx.o: atm_types.o atmfuncs.o atomlist.o files.o parallel.o parallelsubs.o
 m_hsx.o: precision.o siesta_geom.o sys.o
+m_initwf.o: alloc.o densematrix.o fermid.o m_diagon_opt.o m_eo.o m_memory.o
+m_initwf.o: parallel.o parallelsubs.o precision.o sparse_matrices.o sys.o
+m_initwf.o: wavefunctions.o
 m_integrate.o: precision.o
 m_io.o: sys.o
 m_io_s.o: class_OrbitalDistribution.o class_SpData1D.o class_SpData2D.o
 m_io_s.o: class_Sparsity.o parallel.o precision.o
+m_iodipol.o: files.o m_dipol.o m_steps.o
 m_iodm.o: class_OrbitalDistribution.o class_SpData2D.o class_Sparsity.o
 m_iodm.o: m_io_s.o parallel.o
 m_iodm_old.o: alloc.o files.o m_matio.o m_mpi_utils.o parallel.o parallelsubs.o
@@ -418,9 +430,13 @@
 m_iorho.o: alloc.o parallel.o parallelsubs.o precision.o sys.o
 m_iostruct.o: alloc.o files.o m_mpi_utils.o parallel.o precision.o
 m_iostruct.o: siesta_geom.o sys.o units.o
+m_iotddft.o: files.o m_dipol.o m_steps.o parallel.o siesta_options.o
+m_iotddft.o: wavefunctions.o
 m_kinetic.o: precision.o
 m_mat_invert.o: intrinsic_missing.o m_pivot_array.o precision.o
+m_matdiag.o: precision.o
 m_matio.o: alloc.o parallel.o
+m_matswinvers.o: precision.o
 m_memory.o: memoryinfo.o parallel.o sys.o
 m_mesh_node.o: intrinsic_missing.o parallel.o precision.o
 m_mixing.o: class_Data1D.o class_Fstack_Data1D.o parallel.o precision.o
@@ -468,6 +484,15 @@
 m_planewavematrix.o: atomlist.o dhscf.o m_planewavematrixvar.o precision.o
 m_planewavematrix.o: siesta_geom.o sparse_matrices.o
 m_planewavematrixvar.o: precision.o
+m_plot_1d_arr3d.o: precision.o
+m_poisson_fdf.o: m_char.o precision.o
+m_poisson_obc_comm.o: atm_types.o atmfuncs.o m_mpi_stop.o parallel.o
+m_poisson_obc_comm.o: precision.o siesta_geom.o sys.o
+m_poisson_obc_fastconv.o: cellsubs.o m_mpi_stop.o m_poisson_obc_fastconv_r2c.o
+m_poisson_obc_fastconv.o: parallel.o precision.o sys.o
+m_poisson_obc_fastconv_r2c.o: cellsubs.o m_fastconv_dims.o m_plot_1d_arr3d.o
+m_poisson_obc_fastconv_r2c.o: m_poisson_obc_comm.o precision.o sys.o
+m_poisson_switch.o: m_poisson_obc_fastconv.o precision.o sys.o
 m_pulay.o: alloc.o atomlist.o m_mpi_utils.o m_spin.o m_svd.o parallel.o
 m_pulay.o: precision.o siesta_options.o sparse_matrices.o sys.o
 m_redist_spmatrix.o: alloc.o class_Distribution.o
@@ -483,7 +508,8 @@
 m_sparsity_handling.o: class_Data1D.o class_Data2D.o
 m_sparsity_handling.o: class_OrbitalDistribution.o class_SpData1D.o
 m_sparsity_handling.o: class_SpData2D.o class_Sparsity.o geom_helper.o
-m_sparsity_handling.o: m_interpolate.o m_region.o precision.o
+m_sparsity_handling.o: intrinsic_missing.o m_interpolate.o m_region.o
+m_sparsity_handling.o: precision.o
 m_spin.o: alloc.o parallel.o precision.o sys.o units.o
 m_stress.o: precision.o
 m_supercell.o: atom_graph.o class_OrbitalDistribution.o class_SpData2D.o
@@ -595,9 +621,10 @@
 m_ts_options.o: m_ts_hartree.o m_ts_method.o m_ts_mumps_init.o m_ts_tdir.o
 m_ts_options.o: m_ts_weight.o parallel.o precision.o siesta_options.o units.o
 m_ts_pivot.o: alloc.o class_OrbitalDistribution.o class_Sparsity.o
-m_ts_pivot.o: create_Sparsity_SC.o create_Sparsity_Union.o intrinsic_missing.o
-m_ts_pivot.o: m_char.o m_pivot.o m_region.o m_sparsity_handling.o m_ts_debug.o
-m_ts_pivot.o: m_ts_electype.o parallel.o precision.o
+m_ts_pivot.o: create_Sparsity_SC.o create_Sparsity_Union.o geom_helper.o
+m_ts_pivot.o: intrinsic_missing.o m_char.o m_pivot.o m_region.o
+m_ts_pivot.o: m_sparsity_handling.o m_ts_debug.o m_ts_electype.o parallel.o
+m_ts_pivot.o: precision.o
 m_ts_rgn2trimat.o: alloc.o class_OrbitalDistribution.o class_Sparsity.o
 m_ts_rgn2trimat.o: create_Sparsity_Union.o geom_helper.o m_region.o
 m_ts_rgn2trimat.o: m_ts_electype.o m_ts_method.o m_ts_tri_common.o parallel.o
@@ -740,8 +767,8 @@
 rdiag.o: alloc.o m_diagon_opt.o memoryinfo.o parallel.o precision.o sys.o
 rdiag_elpa.o: alloc.o elpa1.o elpa2.o parallel.o precision.o sys.o
 rdiag_mrrr.o: alloc.o parallel.o precision.o sys.o
-read_options.o: files.o m_charge_add.o m_hartree_add.o m_mixing_scf.o m_spin.o
-read_options.o: m_target_stress.o memoryinfo.o parallel.o precision.o
+read_options.o: files.o m_charge_add.o m_cite.o m_hartree_add.o m_mixing_scf.o
+read_options.o: m_spin.o m_target_stress.o memoryinfo.o parallel.o precision.o
 read_options.o: siesta_cml.o siesta_options.o sys.o units.o
 read_xc_info.o: parallel.o precision.o sys.o
 readsp.o: parallel.o precision.o sys.o
@@ -808,30 +835,32 @@
 siesta_end.o: meshdscf.o meshphi.o moremeshsubs.o parallel.o siesta_cml.o
 siesta_end.o: siesta_cmlsubs.o siesta_master.o siesta_options.o
 siesta_end.o: sparse_matrices.o timestamp.o
-siesta_forces.o: atomlist.o class_Fstack_Data1D.o class_SpData2D.o compute_dm.o
-siesta_forces.o: compute_energies.o compute_max_diff.o files.o
-siesta_forces.o: final_H_f_stress.o flook_siesta.o kpoint_grid.o
-siesta_forces.o: m_convergence.o m_energies.o m_forces.o m_iodm_old.o
-siesta_forces.o: m_mixing.o m_mixing_scf.o m_mpi_utils.o m_ncdf_siesta.o
-siesta_forces.o: m_pexsi.o m_pexsi_driver.o m_rhog.o m_spin.o m_steps.o
-siesta_forces.o: m_stress.o m_transiesta.o m_ts_charge.o m_ts_electype.o
-siesta_forces.o: m_ts_global_vars.o m_ts_method.o m_ts_options.o mixer.o
-siesta_forces.o: parallel.o post_scf_work.o precision.o save_density_matrix.o
-siesta_forces.o: scfconvergence_test.o setup_H0.o setup_hamiltonian.o
-siesta_forces.o: siesta_cml.o siesta_dicts.o siesta_geom.o siesta_master.o
-siesta_forces.o: siesta_options.o sparse_matrices.o state_analysis.o
-siesta_forces.o: state_init.o sys.o units.o write_subs.o
+siesta_forces.o: alloc.o atomlist.o class_Fstack_Data1D.o class_SpData2D.o
+siesta_forces.o: compute_dm.o compute_energies.o compute_max_diff.o
+siesta_forces.o: densematrix.o files.o final_H_f_stress.o flook_siesta.o
+siesta_forces.o: kpoint_grid.o m_convergence.o m_energies.o m_eo.o m_evolve.o
+siesta_forces.o: m_forces.o m_gamma.o m_initwf.o m_iodm.o m_iodm_old.o
+siesta_forces.o: m_iotddft.o m_mixing.o m_mixing_scf.o m_mpi_utils.o
+siesta_forces.o: m_ncdf_siesta.o m_pexsi.o m_pexsi_driver.o m_rhog.o m_spin.o
+siesta_forces.o: m_steps.o m_stress.o m_transiesta.o m_ts_charge.o
+siesta_forces.o: m_ts_electype.o m_ts_global_vars.o m_ts_method.o
+siesta_forces.o: m_ts_options.o mixer.o overfsm.o parallel.o post_scf_work.o
+siesta_forces.o: precision.o save_density_matrix.o scfconvergence_test.o
+siesta_forces.o: setup_H0.o setup_hamiltonian.o siesta_cml.o siesta_dicts.o
+siesta_forces.o: siesta_geom.o siesta_master.o siesta_options.o
+siesta_forces.o: sparse_matrices.o state_analysis.o state_init.o sys.o units.o
+siesta_forces.o: wavefunctions.o write_subs.o
 siesta_geom.o: precision.o
 siesta_init.o: alloc.o atomlist.o bands.o bsc_xcmod.o
 siesta_init.o: class_Fstack_Pair_Geometry_SpData2D.o densematrix.o files.o
 siesta_init.o: flook_siesta.o ioxv.o kpoint_grid.o kpoint_pdos.o ksvinit.o
 siesta_init.o: m_cite.o m_energies.o m_eo.o m_fixed.o m_forces.o m_gamma.o
-siesta_init.o: m_iostruct.o m_mpi_utils.o m_new_dm.o m_rmaxh.o m_spin.o
-siesta_init.o: m_steps.o m_supercell.o m_timer.o m_wallclock.o metaforce.o
-siesta_init.o: molecularmechanics.o object_debug.o parallel.o parallelsubs.o
-siesta_init.o: projected_DOS.o siesta_cmlsubs.o siesta_dicts.o siesta_geom.o
-siesta_init.o: siesta_options.o sparse_matrices.o struct_init.o sys.o timer.o
-siesta_init.o: timestamp.o ts_init.o writewave.o zmatrix.o
+siesta_init.o: m_iostruct.o m_mpi_utils.o m_new_dm.o m_poisson_fdf.o m_rmaxh.o
+siesta_init.o: m_spin.o m_steps.o m_supercell.o m_timer.o m_wallclock.o
+siesta_init.o: metaforce.o molecularmechanics.o object_debug.o parallel.o
+siesta_init.o: parallelsubs.o projected_DOS.o siesta_cmlsubs.o siesta_dicts.o
+siesta_init.o: siesta_geom.o siesta_options.o sparse_matrices.o struct_init.o
+siesta_init.o: sys.o timer.o timestamp.o ts_init.o writewave.o zmatrix.o
 siesta_master.o: iopipes.o iosockets.o precision.o sys.o
 siesta_move.o: atomlist.o broyden_optim.o cell_broyden_optim.o
 siesta_move.o: cell_fire_optim.o dynamics.o fire_optim.o flook_siesta.o ioxv.o
@@ -863,18 +892,17 @@
 state_init.o: proximity_check.o siesta_cml.o siesta_geom.o siesta_options.o
 state_init.o: sparse_matrices.o sys.o units.o write_subs.o zmatrix.o
 struct_init.o: alloc.o atmfuncs.o atomlist.o files.o ioxv.o m_exp_coord.o
-struct_init.o: m_iostruct.o m_mpi_utils.o m_steps.o parallel.o periodic_table.o
-struct_init.o: siesta_cml.o siesta_geom.o siesta_master.o siesta_options.o
-struct_init.o: units.o zmatrix.o
+struct_init.o: m_iostruct.o m_mpi_utils.o m_poisson_fdf.o m_poisson_obc_comm.o
+struct_init.o: m_steps.o parallel.o periodic_table.o siesta_cml.o siesta_geom.o
+struct_init.o: siesta_master.o siesta_options.o units.o zmatrix.o
 sys.o: parallel.o
 timer.o: extrae_eventllist.o extrae_module.o m_timer.o parallel.o sys.o
 timer.o: timer_tree.o
 timer_tree.o: m_walltime.o
 transition_rate.o: alloc.o fermid.o parallel.o parallelsubs.o precision.o sys.o
-ts_init.o: m_cite.o m_fixed.o m_os.o m_ts_cctype.o m_ts_charge.o
-ts_init.o: m_ts_electrode.o m_ts_electype.o m_ts_gf.o m_ts_global_vars.o
-ts_init.o: m_ts_kpoints.o m_ts_method.o m_ts_options.o parallel.o
-ts_init.o: siesta_options.o
+ts_init.o: m_fixed.o m_os.o m_ts_cctype.o m_ts_charge.o m_ts_electrode.o
+ts_init.o: m_ts_electype.o m_ts_gf.o m_ts_global_vars.o m_ts_kpoints.o
+ts_init.o: m_ts_method.o m_ts_options.o parallel.o siesta_options.o
 ts_show_regions.o: m_region.o m_ts_electype.o m_ts_method.o parallel.o
 ts_show_regions.o: precision.o units.o
 typecell.o: precision.o
@@ -886,6 +914,7 @@
 vmatsp.o: alloc.o atm_types.o atmfuncs.o atomlist.o listsc.o mesh.o meshdscf.o
 vmatsp.o: meshphi.o precision.o
 vmb.o: m_fixed.o parallel.o precision.o sys.o
+wavefunctions.o: parallel.o precision.o
 write_inp_wannier.o: alloc.o atmfuncs.o atomlist.o m_ntm.o m_orderbands.o
 write_inp_wannier.o: mneighb.o parallel.o parallelsubs.o precision.o
 write_inp_wannier.o: siesta2wannier90.o siesta_geom.o
@@ -1007,6 +1036,7 @@
 m_ksv.o: ksv.o
 m_ksvinit.o: ksvinit.o
 m_ldau.o: ldau.o
+m_linpack.o: linpack.o
 m_local_dos.o: local_DOS.o
 m_matel_registry.o: matel_registry.o
 m_minvec.o: minvec.o
@@ -1051,6 +1081,8 @@
 m_vmat.o: vmat.o
 m_zm_broyden_optim.o: zm_broyden_optim.o
 m_zm_fire_optim.o: zm_fire_optim.o
+matdiagon.o: m_matdiag.o
+matswinversion.o: m_matswinvers.o
 mtprng.o: m_uuid.o
 neighbour.o: mneighb.o
 on_core.o: onmod.o

=== modified file 'Util/VCA/Makefile'
--- Util/VCA/Makefile	2016-08-30 14:02:33 +0000
+++ Util/VCA/Makefile	2016-11-22 08:16:45 +0000
@@ -142,7 +142,7 @@
 bsc_cellxc.o: alloc.o bsc_xcmod.o cellxc_mod.o mesh.o moremeshsubs.o parallel.o
 bsc_cellxc.o: parallelsubs.o precision.o sys.o
 bsc_xcmod.o: parallel.o precision.o sys.o
-cart2frac.o: sys.o
+cart2frac.o: precision.o sys.o
 cdiag.o: alloc.o m_diagon_opt.o memoryinfo.o parallel.o precision.o sys.o
 cell_broyden_optim.o: m_broyddj_nocomm.o parallel.o precision.o sys.o units.o
 cell_broyden_optim.o: zmatrix.o
@@ -154,6 +154,8 @@
 cgvc_zmatrix.o: alloc.o conjgr.o m_mpi_utils.o parallel.o precision.o sys.o
 cgvc_zmatrix.o: units.o zmatrix.o
 cgwf.o: egandd.o onmod.o parallel.o precision.o sys.o
+changebasis.o: alloc.o m_matdiag.o parallel.o parallelsubs.o precision.o
+changebasis.o: sparse_matrices.o sys.o wavefunctions.o
 chemical.o: parallel.o precision.o sys.o
 chempot.o: alloc.o mneighb.o parallel.o parallelsubs.o precision.o sys.o
 chkdim.o: sys.o
@@ -205,17 +207,17 @@
 denmat.o: alloc.o globalise.o m_mpi_utils.o onmod.o precision.o
 denmatlomem.o: alloc.o globalise.o onmod.o precision.o
 densematrix.o: alloc.o precision.o
-detover.o: alloc.o parallel.o parallelsubs.o precision.o
+detover.o: alloc.o linpack.o parallel.o parallelsubs.o precision.o
 dfscf.o: alloc.o atm_types.o atmfuncs.o atomlist.o listsc.o m_spin.o mesh.o
 dfscf.o: meshdscf.o meshphi.o parallel.o parallelsubs.o precision.o sys.o
 dhscf.o: alloc.o atmfuncs.o bsc_xcmod.o cellxc_mod.o delk.o dfscf.o
 dhscf.o: doping_uniform.o files.o forhar.o iogrid_netcdf.o m_charge_add.o
-dhscf.o: m_efield.o m_hartree_add.o m_iorho.o m_mesh_node.o m_ncdf_io.o
-dhscf.o: m_ncdf_siesta.o m_partial_charges.o m_rhog.o m_spin.o
-dhscf.o: m_ts_global_vars.o m_ts_hartree.o m_ts_options.o m_ts_tdir.o
-dhscf.o: m_ts_voltage.o mesh.o meshdscf.o meshsubs.o moremeshsubs.o parallel.o
-dhscf.o: parsing.o precision.o rhofft.o rhoofd.o siesta_options.o sys.o units.o
-dhscf.o: vmat.o
+dhscf.o: m_efield.o m_hartree_add.o m_iorho.o m_iotddft.o m_mesh_node.o
+dhscf.o: m_ncdf_io.o m_ncdf_siesta.o m_partial_charges.o m_poisson_fdf.o
+dhscf.o: m_poisson_switch.o m_rhog.o m_spin.o m_ts_global_vars.o m_ts_hartree.o
+dhscf.o: m_ts_options.o m_ts_tdir.o m_ts_voltage.o mesh.o meshdscf.o meshsubs.o
+dhscf.o: moremeshsubs.o parallel.o parsing.o precision.o rhofft.o rhoofd.o
+dhscf.o: siesta_options.o sys.o units.o vmat.o
 diag2g.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
 diag2k.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
 diag3g.o: fermid.o parallel.o parallelsubs.o precision.o sys.o
@@ -243,7 +245,7 @@
 domain_decom.o: sparse_matrices.o sys.o
 doping_uniform.o: alloc.o m_ntm.o mesh.o parallel.o precision.o sys.o
 dynamics.o: alloc.o atomlist.o files.o ioxv.o m_mpi_utils.o parallel.o
-dynamics.o: precision.o sys.o units.o
+dynamics.o: precision.o siesta_options.o sys.o units.o
 egandd.o: alloc.o denmat.o ener3.o globalise.o gradient.o on_subs.o onmod.o
 egandd.o: onmod.o precision.o sys.o
 eggbox.o: parallel.o precision.o
@@ -252,6 +254,10 @@
 elpa2.o: elpa1.o
 ener3.o: alloc.o globalise.o m_mpi_utils.o onmod.o precision.o
 ener3lomem.o: alloc.o globalise.o m_mpi_utils.o onmod.o precision.o
+evolg.o: alloc.o m_eo.o m_matswinvers.o m_steps.o parallel.o parallelsubs.o
+evolg.o: precision.o siesta_options.o sparse_matrices.o sys.o wavefunctions.o
+evolk.o: m_diagon_opt.o m_matswinvers.o parallel.o precision.o siesta_options.o
+evolk.o: sys.o wavefunctions.o
 extrapolateSpData2D.o: class_Data2D.o class_OrbitalDistribution.o
 extrapolateSpData2D.o: class_SpData2D.o class_Sparsity.o
 extrapolon.o: parallel.o precision.o sys.o
@@ -292,9 +298,9 @@
 initatom.o: atmparams.o atom.o atom_options.o basis_io.o basis_specs.o
 initatom.o: basis_types.o chemical.o electrostatic.o ldau_specs.o m_spin.o
 initatom.o: old_atmfuncs.o precision.o pseudopotential.o
-initparallel.o: alloc.o class_OrbitalDistribution.o domain_decom.o parallel.o
-initparallel.o: parallelsubs.o precision.o siesta_options.o sparse_matrices.o
-initparallel.o: spatial.o sys.o
+initparallel.o: alloc.o class_OrbitalDistribution.o domain_decom.o
+initparallel.o: m_poisson_fdf.o parallel.o parallelsubs.o precision.o
+initparallel.o: siesta_options.o sparse_matrices.o spatial.o sys.o
 io.o: sys.o
 iocg.o: files.o precision.o sys.o
 iodm.o: alloc.o files.o parallel.o parallelsubs.o precision.o sys.o
@@ -353,7 +359,9 @@
 m_efield.o: siesta_geom.o sys.o units.o
 m_energies.o: precision.o
 m_eo.o: precision.o
+m_evolve.o: alloc.o m_memory.o parallel.o parallelsubs.o precision.o sys.o
 m_exp_coord.o: files.o m_os.o parallel.o precision.o units.o
+m_fastconv_dims.o: precision.o
 m_filter.o: bessph.o precision.o radfft.o sorting.o sys.o
 m_fire.o: parallel.o precision.o
 m_fire_mixing.o: alloc.o m_fire_para.o m_mpi_utils.o parallel.o precision.o
@@ -375,10 +383,14 @@
 m_hs_matrix.o: alloc.o cellsubs.o geom_helper.o precision.o sys.o
 m_hsx.o: atm_types.o atmfuncs.o atomlist.o files.o parallel.o parallelsubs.o
 m_hsx.o: precision.o siesta_geom.o sys.o
+m_initwf.o: alloc.o densematrix.o fermid.o m_diagon_opt.o m_eo.o m_memory.o
+m_initwf.o: parallel.o parallelsubs.o precision.o sparse_matrices.o sys.o
+m_initwf.o: wavefunctions.o
 m_integrate.o: precision.o
 m_io.o: sys.o
 m_io_s.o: class_OrbitalDistribution.o class_SpData1D.o class_SpData2D.o
 m_io_s.o: class_Sparsity.o parallel.o precision.o
+m_iodipol.o: files.o m_dipol.o m_steps.o
 m_iodm.o: class_OrbitalDistribution.o class_SpData2D.o class_Sparsity.o
 m_iodm.o: m_io_s.o parallel.o
 m_iodm_old.o: alloc.o files.o m_matio.o m_mpi_utils.o parallel.o parallelsubs.o
@@ -386,9 +398,13 @@
 m_iorho.o: alloc.o parallel.o parallelsubs.o precision.o sys.o
 m_iostruct.o: alloc.o files.o m_mpi_utils.o parallel.o precision.o
 m_iostruct.o: siesta_geom.o sys.o units.o
+m_iotddft.o: files.o m_dipol.o m_steps.o parallel.o siesta_options.o
+m_iotddft.o: wavefunctions.o
 m_kinetic.o: precision.o
 m_mat_invert.o: intrinsic_missing.o m_pivot_array.o precision.o
+m_matdiag.o: precision.o
 m_matio.o: alloc.o parallel.o
+m_matswinvers.o: precision.o
 m_memory.o: memoryinfo.o parallel.o sys.o
 m_mesh_node.o: intrinsic_missing.o parallel.o precision.o
 m_mixing.o: class_Data1D.o class_Fstack_Data1D.o parallel.o precision.o
@@ -436,6 +452,15 @@
 m_planewavematrix.o: atomlist.o dhscf.o m_planewavematrixvar.o precision.o
 m_planewavematrix.o: siesta_geom.o sparse_matrices.o
 m_planewavematrixvar.o: precision.o
+m_plot_1d_arr3d.o: precision.o
+m_poisson_fdf.o: m_char.o precision.o
+m_poisson_obc_comm.o: atm_types.o atmfuncs.o m_mpi_stop.o parallel.o
+m_poisson_obc_comm.o: precision.o siesta_geom.o sys.o
+m_poisson_obc_fastconv.o: cellsubs.o m_mpi_stop.o m_poisson_obc_fastconv_r2c.o
+m_poisson_obc_fastconv.o: parallel.o precision.o sys.o
+m_poisson_obc_fastconv_r2c.o: cellsubs.o m_fastconv_dims.o m_plot_1d_arr3d.o
+m_poisson_obc_fastconv_r2c.o: m_poisson_obc_comm.o precision.o sys.o
+m_poisson_switch.o: m_poisson_obc_fastconv.o precision.o sys.o
 m_pulay.o: alloc.o atomlist.o m_mpi_utils.o m_spin.o m_svd.o parallel.o
 m_pulay.o: precision.o siesta_options.o sparse_matrices.o sys.o
 m_redist_spmatrix.o: alloc.o class_Distribution.o
@@ -451,7 +476,8 @@
 m_sparsity_handling.o: class_Data1D.o class_Data2D.o
 m_sparsity_handling.o: class_OrbitalDistribution.o class_SpData1D.o
 m_sparsity_handling.o: class_SpData2D.o class_Sparsity.o geom_helper.o
-m_sparsity_handling.o: m_interpolate.o m_region.o precision.o
+m_sparsity_handling.o: intrinsic_missing.o m_interpolate.o m_region.o
+m_sparsity_handling.o: precision.o
 m_spin.o: alloc.o parallel.o precision.o sys.o units.o
 m_stress.o: precision.o
 m_supercell.o: atom_graph.o class_OrbitalDistribution.o class_SpData2D.o
@@ -563,9 +589,10 @@
 m_ts_options.o: m_ts_hartree.o m_ts_method.o m_ts_mumps_init.o m_ts_tdir.o
 m_ts_options.o: m_ts_weight.o parallel.o precision.o siesta_options.o units.o
 m_ts_pivot.o: alloc.o class_OrbitalDistribution.o class_Sparsity.o
-m_ts_pivot.o: create_Sparsity_SC.o create_Sparsity_Union.o intrinsic_missing.o
-m_ts_pivot.o: m_char.o m_pivot.o m_region.o m_sparsity_handling.o m_ts_debug.o
-m_ts_pivot.o: m_ts_electype.o parallel.o precision.o
+m_ts_pivot.o: create_Sparsity_SC.o create_Sparsity_Union.o geom_helper.o
+m_ts_pivot.o: intrinsic_missing.o m_char.o m_pivot.o m_region.o
+m_ts_pivot.o: m_sparsity_handling.o m_ts_debug.o m_ts_electype.o parallel.o
+m_ts_pivot.o: precision.o
 m_ts_rgn2trimat.o: alloc.o class_OrbitalDistribution.o class_Sparsity.o
 m_ts_rgn2trimat.o: create_Sparsity_Union.o geom_helper.o m_region.o
 m_ts_rgn2trimat.o: m_ts_electype.o m_ts_method.o m_ts_tri_common.o parallel.o
@@ -708,8 +735,8 @@
 rdiag.o: alloc.o m_diagon_opt.o memoryinfo.o parallel.o precision.o sys.o
 rdiag_elpa.o: alloc.o elpa1.o elpa2.o parallel.o precision.o sys.o
 rdiag_mrrr.o: alloc.o parallel.o precision.o sys.o
-read_options.o: files.o m_charge_add.o m_hartree_add.o m_mixing_scf.o m_spin.o
-read_options.o: m_target_stress.o memoryinfo.o parallel.o precision.o
+read_options.o: files.o m_charge_add.o m_cite.o m_hartree_add.o m_mixing_scf.o
+read_options.o: m_spin.o m_target_stress.o memoryinfo.o parallel.o precision.o
 read_options.o: siesta_cml.o siesta_options.o sys.o units.o
 read_xc_info.o: parallel.o precision.o sys.o
 readsp.o: parallel.o precision.o sys.o
@@ -776,30 +803,32 @@
 siesta_end.o: meshdscf.o meshphi.o moremeshsubs.o parallel.o siesta_cml.o
 siesta_end.o: siesta_cmlsubs.o siesta_master.o siesta_options.o
 siesta_end.o: sparse_matrices.o timestamp.o
-siesta_forces.o: atomlist.o class_Fstack_Data1D.o class_SpData2D.o compute_dm.o
-siesta_forces.o: compute_energies.o compute_max_diff.o files.o
-siesta_forces.o: final_H_f_stress.o flook_siesta.o kpoint_grid.o
-siesta_forces.o: m_convergence.o m_energies.o m_forces.o m_iodm_old.o
-siesta_forces.o: m_mixing.o m_mixing_scf.o m_mpi_utils.o m_ncdf_siesta.o
-siesta_forces.o: m_pexsi.o m_pexsi_driver.o m_rhog.o m_spin.o m_steps.o
-siesta_forces.o: m_stress.o m_transiesta.o m_ts_charge.o m_ts_electype.o
-siesta_forces.o: m_ts_global_vars.o m_ts_method.o m_ts_options.o mixer.o
-siesta_forces.o: parallel.o post_scf_work.o precision.o save_density_matrix.o
-siesta_forces.o: scfconvergence_test.o setup_H0.o setup_hamiltonian.o
-siesta_forces.o: siesta_cml.o siesta_dicts.o siesta_geom.o siesta_master.o
-siesta_forces.o: siesta_options.o sparse_matrices.o state_analysis.o
-siesta_forces.o: state_init.o sys.o units.o write_subs.o
+siesta_forces.o: alloc.o atomlist.o class_Fstack_Data1D.o class_SpData2D.o
+siesta_forces.o: compute_dm.o compute_energies.o compute_max_diff.o
+siesta_forces.o: densematrix.o files.o final_H_f_stress.o flook_siesta.o
+siesta_forces.o: kpoint_grid.o m_convergence.o m_energies.o m_eo.o m_evolve.o
+siesta_forces.o: m_forces.o m_gamma.o m_initwf.o m_iodm.o m_iodm_old.o
+siesta_forces.o: m_iotddft.o m_mixing.o m_mixing_scf.o m_mpi_utils.o
+siesta_forces.o: m_ncdf_siesta.o m_pexsi.o m_pexsi_driver.o m_rhog.o m_spin.o
+siesta_forces.o: m_steps.o m_stress.o m_transiesta.o m_ts_charge.o
+siesta_forces.o: m_ts_electype.o m_ts_global_vars.o m_ts_method.o
+siesta_forces.o: m_ts_options.o mixer.o overfsm.o parallel.o post_scf_work.o
+siesta_forces.o: precision.o save_density_matrix.o scfconvergence_test.o
+siesta_forces.o: setup_H0.o setup_hamiltonian.o siesta_cml.o siesta_dicts.o
+siesta_forces.o: siesta_geom.o siesta_master.o siesta_options.o
+siesta_forces.o: sparse_matrices.o state_analysis.o state_init.o sys.o units.o
+siesta_forces.o: wavefunctions.o write_subs.o
 siesta_geom.o: precision.o
 siesta_init.o: alloc.o atomlist.o bands.o bsc_xcmod.o
 siesta_init.o: class_Fstack_Pair_Geometry_SpData2D.o densematrix.o files.o
 siesta_init.o: flook_siesta.o ioxv.o kpoint_grid.o kpoint_pdos.o ksvinit.o
 siesta_init.o: m_cite.o m_energies.o m_eo.o m_fixed.o m_forces.o m_gamma.o
-siesta_init.o: m_iostruct.o m_mpi_utils.o m_new_dm.o m_rmaxh.o m_spin.o
-siesta_init.o: m_steps.o m_supercell.o m_timer.o m_wallclock.o metaforce.o
-siesta_init.o: molecularmechanics.o object_debug.o parallel.o parallelsubs.o
-siesta_init.o: projected_DOS.o siesta_cmlsubs.o siesta_dicts.o siesta_geom.o
-siesta_init.o: siesta_options.o sparse_matrices.o struct_init.o sys.o timer.o
-siesta_init.o: timestamp.o ts_init.o writewave.o zmatrix.o
+siesta_init.o: m_iostruct.o m_mpi_utils.o m_new_dm.o m_poisson_fdf.o m_rmaxh.o
+siesta_init.o: m_spin.o m_steps.o m_supercell.o m_timer.o m_wallclock.o
+siesta_init.o: metaforce.o molecularmechanics.o object_debug.o parallel.o
+siesta_init.o: parallelsubs.o projected_DOS.o siesta_cmlsubs.o siesta_dicts.o
+siesta_init.o: siesta_geom.o siesta_options.o sparse_matrices.o struct_init.o
+siesta_init.o: sys.o timer.o timestamp.o ts_init.o writewave.o zmatrix.o
 siesta_master.o: iopipes.o iosockets.o precision.o sys.o
 siesta_move.o: atomlist.o broyden_optim.o cell_broyden_optim.o
 siesta_move.o: cell_fire_optim.o dynamics.o fire_optim.o flook_siesta.o ioxv.o
@@ -831,18 +860,17 @@
 state_init.o: proximity_check.o siesta_cml.o siesta_geom.o siesta_options.o
 state_init.o: sparse_matrices.o sys.o units.o write_subs.o zmatrix.o
 struct_init.o: alloc.o atmfuncs.o atomlist.o files.o ioxv.o m_exp_coord.o
-struct_init.o: m_iostruct.o m_mpi_utils.o m_steps.o parallel.o periodic_table.o
-struct_init.o: siesta_cml.o siesta_geom.o siesta_master.o siesta_options.o
-struct_init.o: units.o zmatrix.o
+struct_init.o: m_iostruct.o m_mpi_utils.o m_poisson_fdf.o m_poisson_obc_comm.o
+struct_init.o: m_steps.o parallel.o periodic_table.o siesta_cml.o siesta_geom.o
+struct_init.o: siesta_master.o siesta_options.o units.o zmatrix.o
 sys.o: parallel.o
 timer.o: extrae_eventllist.o extrae_module.o m_timer.o parallel.o sys.o
 timer.o: timer_tree.o
 timer_tree.o: m_walltime.o
 transition_rate.o: alloc.o fermid.o parallel.o parallelsubs.o precision.o sys.o
-ts_init.o: m_cite.o m_fixed.o m_os.o m_ts_cctype.o m_ts_charge.o
-ts_init.o: m_ts_electrode.o m_ts_electype.o m_ts_gf.o m_ts_global_vars.o
-ts_init.o: m_ts_kpoints.o m_ts_method.o m_ts_options.o parallel.o
-ts_init.o: siesta_options.o
+ts_init.o: m_fixed.o m_os.o m_ts_cctype.o m_ts_charge.o m_ts_electrode.o
+ts_init.o: m_ts_electype.o m_ts_gf.o m_ts_global_vars.o m_ts_kpoints.o
+ts_init.o: m_ts_method.o m_ts_options.o parallel.o siesta_options.o
 ts_show_regions.o: m_region.o m_ts_electype.o m_ts_method.o parallel.o
 ts_show_regions.o: precision.o units.o
 typecell.o: precision.o
@@ -854,6 +882,7 @@
 vmatsp.o: alloc.o atm_types.o atmfuncs.o atomlist.o listsc.o mesh.o meshdscf.o
 vmatsp.o: meshphi.o precision.o
 vmb.o: m_fixed.o parallel.o precision.o sys.o
+wavefunctions.o: parallel.o precision.o
 write_inp_wannier.o: alloc.o atmfuncs.o atomlist.o m_ntm.o m_orderbands.o
 write_inp_wannier.o: mneighb.o parallel.o parallelsubs.o precision.o
 write_inp_wannier.o: siesta2wannier90.o siesta_geom.o
@@ -974,6 +1003,7 @@
 m_ksv.o: ksv.o
 m_ksvinit.o: ksvinit.o
 m_ldau.o: ldau.o
+m_linpack.o: linpack.o
 m_local_dos.o: local_DOS.o
 m_matel_registry.o: matel_registry.o
 m_minvec.o: minvec.o
@@ -1018,6 +1048,8 @@
 m_vmat.o: vmat.o
 m_zm_broyden_optim.o: zm_broyden_optim.o
 m_zm_fire_optim.o: zm_fire_optim.o
+matdiagon.o: m_matdiag.o
+matswinversion.o: m_matswinvers.o
 mtprng.o: m_uuid.o
 neighbour.o: mneighb.o
 on_core.o: onmod.o

=== modified file 'version.info'
--- version.info	2016-10-20 12:18:03 +0000
+++ version.info	2016-11-22 08:16:45 +0000
@@ -1,1 +1,1 @@
-trunk-572--kovalp-manual-1
+trunk-572--kovalp-manual-1--npapior-1