add linear constant coefficients equation test

Author: szaghi

src/lib/foodie.f90

@@ -219,18 +219,18 @@ subroutine foodie_integrator_factory(scheme, integrator, stages, tolerance, nu,
     type is(integrator_runge_kutta_emd)
       call integrator%initialize(scheme=scheme, tolerance=tolerance)
     endselect
-  elseif (index(trim(adjustl(scheme)), trim(int_runge_kutta_ls%class_name())) > 0) then
-    allocate(integrator_runge_kutta_ls :: integrator)
-    select type(integrator)
-    type is(integrator_runge_kutta_ls)
-      call integrator%initialize(scheme=scheme)
-    endselect
   elseif (index(trim(adjustl(scheme)), trim(int_runge_kutta_lssp%class_name())) > 0) then
     allocate(integrator_runge_kutta_lssp :: integrator)
     select type(integrator)
     type is(integrator_runge_kutta_lssp)
       call integrator%initialize(scheme=scheme, stages=stages)
     endselect
+  elseif (index(trim(adjustl(scheme)), trim(int_runge_kutta_ls%class_name())) > 0) then
+    allocate(integrator_runge_kutta_ls :: integrator)
+    select type(integrator)
+    type is(integrator_runge_kutta_ls)
+      call integrator%initialize(scheme=scheme)
+    endselect
   elseif (index(trim(adjustl(scheme)), trim(int_runge_kutta_ssp%class_name())) > 0) then
     allocate(integrator_runge_kutta_ssp :: integrator)
     select type(integrator)

src/lib/foodie_integrand_object.f90

@@ -16,7 +16,8 @@ module foodie_integrand_object
 #endif
   contains
     ! public deferred procedures that concrete integrand-field must implement
-    procedure(time_derivative), pass(self), deferred, public :: t !< Time derivative, residuals.
+    procedure(integrand_dimension_interface), pass(self), deferred, public :: integrand_dimension !< Return integrand dimension.
+    procedure(time_derivative),               pass(self), deferred, public :: t                   !< Time derivative, residuals.
     ! operators
     procedure(local_error_operator), pass(lhs), deferred, public :: local_error !< `||integrand - integrand||` operator.
     generic, public :: operator(.lterror.) => local_error !< Estimate local truncation error.
@@ -62,7 +63,14 @@ module foodie_integrand_object
 endtype integrand_object
 
 abstract interface
-  !< Abstract type bound procedures necessary for implementing a concrete extension of [[integrand_object]].
+   !< Abstract type bound procedures necessary for implementing a concrete extension of [[integrand_object]].
+
+   pure function integrand_dimension_interface(self) result(integrand_dimension)
+   !< Return integrand dimension.
+   import :: integrand_object, I_P
+   class(integrand_object), intent(in) :: self                !< Integrand.
+   integer(I_P)                        :: integrand_dimension !< Integrand dimension.
+   endfunction integrand_dimension_interface
 
   function time_derivative(self, t) result(dState_dt)
   !< Time derivative function of integrand class, i.e. the residuals function.

src/lib/foodie_integrator_backward_differentiation_formula.f90

@@ -279,7 +279,7 @@ subroutine integrate(self, U, previous, Dt, t, iterations, autoupdate)
   do s=1, self%steps - 1
     delta = delta + (previous(s) * (-self%a(s)))
   enddo
-  do s=1, iterations
+  do s=1, iterations_
     U = delta + (U%t(t=t(self%steps) + Dt) * (Dt * self%b))
   enddo
   if (autoupdate_) call self%update_previous(U=U, previous=previous)

src/tests/lcce/fobos

@@ -0,0 +1,127 @@
+[modes]
+modes = gnu gnu-debug gnu-coverage
+        intel intel-debug
+        errors-analysis
+        gnu-pure gnu-debug-pure gnu-coverage-pure
+        intel-pure intel-debug-pure
+
+[common-variables]
+$CSTATIC_GNU = -cpp -c -frealloc-lhs
+$CSTATIC_INT = -cpp -c -assume realloc_lhs
+$DEBUG_GNU   = -Og -g3 -Warray-bounds -Wcharacter-truncation -Wline-truncation -Wimplicit-interface -Wimplicit-procedure -Wunderflow -Wuninitialized -fcheck=all -fmodule-private -ffree-line-length-132 -fimplicit-none -fbacktrace -fdump-core -finit-real=nan -std=f2008 -fall-intrinsics
+$DEBUG_INT   = -O0 -debug all -check all -warn all -extend-source 132 -traceback -gen-interfaces -fp-stack-check -fstack-protector-all -ftrapuv -no-ftz -std08
+$OPTIMIZE    = -O2
+$EXDIRS      = src/tests/parallel/ src/tests/regression/
+               PENF/src/tests/
+               FACE/src/tests/ FACE/src/third_party/
+               FLAP/src/tests/ FLAP/src/third_party/
+               WenOOF/src/tests/ WenOOF/src/third_party/
+
+# main modes
+[gnu]
+help     = Build test in release mode with GNU gfortran
+compiler = gnu
+cflags   = $CSTATIC_GNU $OPTIMIZE
+lflags   = $OPTIMIZE
+template = template-common
+
+[gnu-debug]
+help     = Build test in debug mode with GNU gfortran
+compiler = gnu
+cflags   = $CSTATIC_GNU $DEBUG_GNU
+lflags   = $DEBUG_GNU
+template = template-common
+
+[gnu-coverage]
+help     = Build test in release coverage mode with GNU gfortran
+compiler = gnu
+cflags   = $CSTATIC_GNU $DEBUG_GNU
+lflags   = $DEBUG_GNU
+coverage = True
+template = template-common
+
+[intel]
+help     = Build test in release mode with Intel Fortran
+compiler = intel
+cflags   = $CSTATIC_INT $OPTIMIZE
+lflags   = $OPTIMIZE
+template = template-common
+
+[intel-debug]
+help     = Build test in debug mode with Intel Fortran
+compiler = intel
+cflags   = $CSTATIC_INT $DEBUG_INT
+lflags   = $DEBUG_INT
+template = template-common
+
+[errors-analysis]
+help            = Build test in release mode (for errors analysis) with GNU gfortran
+compiler        = gnu
+cflags          = $CSTATIC_GNU $OPTIMIZE
+lflags          = $OPTIMIZE
+cflags_heritage = True
+target          = src/tests/lcc/foodie_test_lcce.f90
+exclude_dirs    = $EXDIRS
+mod_dir         = mod
+obj_dir         = obj
+src             = src/
+colors          = True
+quiet           = False
+log             = True
+jobs            = 10
+
+# pure modes for testing FOODIE ADT with "pure" procedures
+[gnu-pure]
+help     = Build test in release mode with GNU gfortran
+compiler = gnu
+cflags   = $CSTATIC_GNU $OPTIMIZE
+lflags   = $OPTIMIZE
+preproc  = -DPURE
+template = template-common
+
+[gnu-debug-pure]
+help     = Build test in debug mode with GNU gfortran
+compiler = gnu
+cflags   = $CSTATIC_GNU $DEBUG_GNU
+lflags   = $DEBUG_GNU
+preproc  = -DPURE
+template = template-common
+
+[gnu-coverage-pure]
+help     = Build test in release coverage mode with GNU gfortran
+compiler = gnu
+cflags   = $CSTATIC_GNU $OPTIMIZE
+lflags   = $OPTIMIZE
+preproc  = -DPURE
+coverage = True
+template = template-common
+
+[intel-pure]
+help     = Build test in release mode with Intel Fortran
+compiler = intel
+cflags   = $CSTATIC_INT $OPTIMIZE
+lflags   = $OPTIMIZE
+preproc  = -DPURE
+template = template-common
+
+[intel-debug-pure]
+help     = Build test in debug mode with Intel Fortran
+compiler = intel
+cflags   = $CSTATIC_INT $DEBUG_INT
+lflags   = $DEBUG_INT
+preproc  = -DPURE
+template = template-common
+
+# templates
+[template-common]
+cflags_heritage = True
+build_dir       = build/tests/lcc/
+mod_dir         = mod
+obj_dir         = obj
+src             = src/
+target          = src/tests/lcc/foodie_test_lcce.f90
+exclude_dirs    = $EXDIRS
+colors          = True
+quiet           = False
+log             = True
+jobs            = 10