Complete first skeleton

src/bsemol/GNUmakefile

@@ -1,4 +1,5 @@
-OBJ = bsemol.o bse_input.o
+OBJ = bsemol.o bse_input.o bse_init.o bse_analytic.o \
+      bse_davidson.o bse_lanczos.o bse_defaults.o
 
 USES_BLAS =
 OBJ_OPTIMIZE =

src/bsemol/bse_analytic.F

@@ -0,0 +1,9 @@
+      subroutine bse_analytic(pars)
+
+      implicit none
+
+#include "bse.fh"
+
+      type(bse_params_t) :: pars
+
+      end subroutine bse_analytic

src/bsemol/bse_davidson.F

@@ -0,0 +1,9 @@
+      subroutine bse_davidson(pars)
+
+      implicit none
+
+#include "bse.fh"
+
+      type(bse_params_t) :: pars
+
+      end subroutine bse_davidson

src/bsemol/bse_defaults.F

@@ -0,0 +1,6 @@
+      subroutine bse_defaults(rtdb)
+      implicit none
+
+      integer rtdb
+
+      end subroutine bse_defaults

src/bsemol/bse_init.F

@@ -0,0 +1,9 @@
+      subroutine bse_init(pars)
+
+      implicit none
+
+#include "bse.fh"
+
+      type(bse_params_t) :: pars
+
+      end subroutine bse_init

src/bsemol/bse_lanczos.F

@@ -0,0 +1,9 @@
+      subroutine bse_lanczos(pars)
+
+      implicit none
+
+#include "bse.fh"
+
+      type(bse_params_t) :: pars
+
+      end subroutine bse_lanczos

src/bsemol/bsemol.F

@@ -55,6 +55,8 @@ c
         
         if (params%davidson) then
           call bse_davidson(params)
+        elseif (params%lanczos) then
+          call bse_lanczos(params)
         else
           call bse_analytic(params)
         endif