MAINT: Add makefile

Makefile

@@ -0,0 +1,14 @@
+##
+# featom
+#
+# @file
+# @version 0.1
+
+app/.f90.o:
+	$(FC) $(FFLAGS) -c $<
+conv: conv.o
+	$(FC) $(FFLAGS) -o $@ $^
+
+include Makefile.dep
+
+# end

Makefile.dep

@@ -0,0 +1,203 @@
+# This file is generated automatically. DO NOT EDIT!
+
+app/conv :  \
+	app/conv.o \
+	src/constants.o \
+	src/dirac.o \
+	src/energies.o \
+	src/fe.o \
+	src/feutils.o \
+	src/gjp_gw.o \
+	src/graphs.o \
+	src/hartree_screening.o \
+	src/lapack.o \
+	src/linalg.o \
+	src/mesh.o \
+	src/mixings.o \
+	src/quadrature.o \
+	src/schroed_glob.o \
+	src/solvers.o \
+	src/states.o \
+	src/string_utils.o \
+	src/types.o \
+	src/xc.o
+
+app/conv_potential :  \
+	app/conv_potential.o \
+	src/constants.o \
+	src/dirac.o \
+	src/energies.o \
+	src/fe.o \
+	src/feutils.o \
+	src/gjp_gw.o \
+	src/graphs_potential.o \
+	src/hartree_screening.o \
+	src/lapack.o \
+	src/linalg.o \
+	src/mesh.o \
+	src/mixings.o \
+	src/quadrature.o \
+	src/schroed_dirac_solver.o \
+	src/schroed_glob.o \
+	src/solvers.o \
+	src/states.o \
+	src/string_utils.o \
+	src/types.o \
+	src/xc.o
+
+app/gpd_coulomb_schroed_nelements :  \
+	app/gpd_coulomb_schroed_nelements.o \
+	src/constants.o \
+	src/dirac.o \
+	src/energies.o \
+	src/fe.o \
+	src/feutils.o \
+	src/gjp_gw.o \
+	src/graphs_potential.o \
+	src/hartree_screening.o \
+	src/lapack.o \
+	src/linalg.o \
+	src/mesh.o \
+	src/mixings.o \
+	src/quadrature.o \
+	src/schroed_dirac_solver.o \
+	src/schroed_glob.o \
+	src/solvers.o \
+	src/states.o \
+	src/string_utils.o \
+	src/types.o \
+	src/xc.o
+
+src/main :  \
+	src/main.o
+
+src/conv.o :  \
+	src/graphs.o
+
+src/conv_potential.o :  \
+	src/graphs_potential.o
+
+src/gpd_coulomb_schroed_nelements.o :  \
+	src/graphs_potential.o
+
+src/main.o : 
+
+src/Makefile.o : 
+
+src/constants.o :  \
+	src/types.o
+
+src/dirac.o :  \
+	src/constants.o \
+	src/energies.o \
+	src/fe.o \
+	src/feutils.o \
+	src/gjp_gw.o \
+	src/hartree_screening.o \
+	src/linalg.o \
+	src/mesh.o \
+	src/mixings.o \
+	src/states.o \
+	src/types.o \
+	src/xc.o
+
+src/energies.o :  \
+	src/constants.o \
+	src/types.o
+
+src/fe.o :  \
+	src/constants.o \
+	src/feutils.o \
+	src/types.o
+
+src/feutils.o :  \
+	src/quadrature.o \
+	src/solvers.o \
+	src/types.o
+
+src/gjp_gw.o :  \
+	src/lapack.o \
+	src/types.o
+
+src/graphs.o :  \
+	src/dirac.o \
+	src/feutils.o \
+	src/mesh.o \
+	src/schroed_glob.o \
+	src/string_utils.o \
+	src/types.o
+
+src/graphs_potential.o :  \
+	src/dirac.o \
+	src/fe.o \
+	src/feutils.o \
+	src/linalg.o \
+	src/mesh.o \
+	src/schroed_dirac_solver.o \
+	src/schroed_glob.o \
+	src/string_utils.o \
+	src/types.o
+
+src/hartree_screening.o :  \
+	src/feutils.o \
+	src/solvers.o \
+	src/types.o
+
+src/lapack.o : 
+
+src/linalg.o :  \
+	src/lapack.o \
+	src/types.o
+
+src/mesh.o :  \
+	src/types.o
+
+src/mixings.o :  \
+	src/linalg.o \
+	src/types.o
+
+src/quadrature.o :  \
+	src/types.o
+
+src/schroed_dirac_solver.o :  \
+	src/constants.o \
+	src/dirac.o \
+	src/fe.o \
+	src/feutils.o \
+	src/gjp_gw.o \
+	src/linalg.o \
+	src/mesh.o \
+	src/schroed_glob.o \
+	src/string_utils.o \
+	src/types.o
+
+src/schroed_glob.o :  \
+	src/constants.o \
+	src/energies.o \
+	src/fe.o \
+	src/feutils.o \
+	src/hartree_screening.o \
+	src/lapack.o \
+	src/linalg.o \
+	src/mesh.o \
+	src/mixings.o \
+	src/solvers.o \
+	src/states.o \
+	src/types.o \
+	src/xc.o
+
+src/solvers.o :  \
+	src/lapack.o \
+	src/types.o
+
+src/states.o :  \
+	src/types.o
+
+src/string_utils.o :  \
+	src/types.o
+
+src/types.o : 
+
+src/xc.o :  \
+	src/constants.o \
+	src/types.o

app/conv.f90

@@ -37,4 +37,4 @@ program conv
 
 call run_convergence(study_type, dirac_int, p_or_Ne, ".")
 
-end program
+end program conv

app/conv_potential.f90

@@ -48,4 +48,4 @@ program conv_potential
 call run_convergence_potential(study_type, dirac_int, p_or_Ne, &
      potential_type, alpha_int, ".")
 
-end program
+end program conv_potential

src/Makefile.dep

@@ -0,0 +1,119 @@
+# This file is generated automatically. DO NOT EDIT!
+
+constants.o :  \
+	types.o
+
+dirac.o :  \
+	constants.o \
+	energies.o \
+	fe.o \
+	feutils.o \
+	gjp_gw.o \
+	hartree_screening.o \
+	linalg.o \
+	mesh.o \
+	mixings.o \
+	states.o \
+	types.o \
+	xc.o
+
+energies.o :  \
+	constants.o \
+	types.o
+
+fe.o :  \
+	constants.o \
+	feutils.o \
+	types.o
+
+feutils.o :  \
+	quadrature.o \
+	solvers.o \
+	types.o
+
+gjp_gw.o :  \
+	lapack.o \
+	types.o
+
+graphs.o :  \
+	dirac.o \
+	feutils.o \
+	mesh.o \
+	schroed_glob.o \
+	string_utils.o \
+	types.o
+
+graphs_potential.o :  \
+	dirac.o \
+	fe.o \
+	feutils.o \
+	linalg.o \
+	mesh.o \
+	schroed_dirac_solver.o \
+	schroed_glob.o \
+	string_utils.o \
+	types.o
+
+hartree_screening.o :  \
+	feutils.o \
+	solvers.o \
+	types.o
+
+lapack.o : 
+
+linalg.o :  \
+	lapack.o \
+	types.o
+
+mesh.o :  \
+	types.o
+
+mixings.o :  \
+	linalg.o \
+	types.o
+
+quadrature.o :  \
+	types.o
+
+schroed_dirac_solver.o :  \
+	constants.o \
+	dirac.o \
+	fe.o \
+	feutils.o \
+	gjp_gw.o \
+	linalg.o \
+	mesh.o \
+	schroed_glob.o \
+	string_utils.o \
+	types.o
+
+schroed_glob.o :  \
+	constants.o \
+	energies.o \
+	fe.o \
+	feutils.o \
+	hartree_screening.o \
+	lapack.o \
+	linalg.o \
+	mesh.o \
+	mixings.o \
+	solvers.o \
+	states.o \
+	types.o \
+	xc.o
+
+solvers.o :  \
+	lapack.o \
+	types.o
+
+states.o :  \
+	types.o
+
+string_utils.o :  \
+	types.o
+
+types.o : 
+
+xc.o :  \
+	constants.o \
+	types.o

src/constants.f90

@@ -27,4 +27,4 @@ module constants
 real(dp), parameter :: c_2002 = 137.03599911_dp
 real(dp), parameter :: c_1998 = 137.03599976_dp
 real(dp), parameter :: c_1986 = 137.0359895_dp ! To compare with dftatom
-end module
+end module constants

src/dirac.f90

@@ -315,4 +315,4 @@ subroutine csolve_dirac(Z, p, xiq, wtq, xe, eps, Nq, Ne, nenergies, energies, &
 energies = energies2
 end subroutine
 
-end module
+end module dirac

src/energies.f90

@@ -101,4 +101,4 @@ function get_hydrogen_energies(Z, no) result(E)
 end do
 end function
 
-end module
+end module energies

src/fe.f90

@@ -384,4 +384,4 @@ subroutine assemble_radial_dirac_SH(V, kappa, xin, xe, ib, xiq, wtq, xiq1, wtq1,
 end subroutine
 
 
-end module
+end module fe

src/feutils.f90

@@ -717,4 +717,4 @@ real(dp) function integrate2(xe, xiq, xiq_gj, wtq, wtq_gj, zeta, uq) result(r)
 end do
 end function
 
-end module
+end module feutils

src/graphs.f90

@@ -158,4 +158,4 @@ subroutine total_energy(Z, rmax, Ne, a, p, Nq, DOFs, Etot)
 
 end subroutine
 
-end module
+end module graphs

src/graphs_potential.f90

@@ -170,4 +170,4 @@ subroutine run_convergence_potential(study_type, dirac_int, p_or_Ne, &
 close(u)
 end subroutine
 
-end module
+end module graphs_potential

src/hartree_screening.f90

@@ -328,4 +328,4 @@ function hartree_potential_gj(p, xe, xiq, wtq, xiq1, wtq1, zeta, f0) result(uq)
 uq = Yq/xq
 end function
 
-end module
+end module hartree_screening