read_input_parameters.F90 3.35 KB
Newer Older
Andreas Marek's avatar
Andreas Marek committed
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
!    This file is part of ELPA.
!
!    The ELPA library was originally created by the ELPA consortium,
!    consisting of the following organizations:
!
!    - Rechenzentrum Garching der Max-Planck-Gesellschaft (RZG),
!    - Bergische Universität Wuppertal, Lehrstuhl für angewandte
!      Informatik,
!    - Technische Universität München, Lehrstuhl für Informatik mit
!      Schwerpunkt Wissenschaftliches Rechnen ,
!    - Fritz-Haber-Institut, Berlin, Abt. Theorie,
!    - Max-Plack-Institut für Mathematik in den Naturwissenschaften,
!      Leipzig, Abt. Komplexe Strukutren in Biologie und Kognition,
!      and
!    - IBM Deutschland GmbH
!
!
!    More information can be found here:
!    http://elpa.rzg.mpg.de/
!
!    ELPA is free software: you can redistribute it and/or modify
!    it under the terms of the version 3 of the license of the
!    GNU Lesser General Public License as published by the Free
!    Software Foundation.
!
!    ELPA is distributed in the hope that it will be useful,
!    but WITHOUT ANY WARRANTY; without even the implied warranty of
!    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
!    GNU Lesser General Public License for more details.
!
!    You should have received a copy of the GNU Lesser General Public License
!    along with ELPA.  If not, see <http://www.gnu.org/licenses/>
!
!    ELPA reflects a substantial effort on the part of the original
!    ELPA consortium, and we ask you to respect the spirit of the
!    license that we chose: i.e., please contribute any changes you
!    may have back to the original ELPA library distribution, and keep
!    any derivatives of ELPA under the same license that we chose for
!    the original distribution, the GNU Lesser General Public License.
!
!
#include "config-f90.h"
module mod_read_input_parameters

  contains

    subroutine read_input_parameters(na, nev, nblk, write_to_file)
48
      use ELPA_utilities, only : error_unit
Andreas Marek's avatar
Andreas Marek committed
49
      implicit none
50
51
52
53
      include 'mpif.h'

      integer, intent(out) :: na, nev, nblk
      logical, intent(out) :: write_to_file
Andreas Marek's avatar
Andreas Marek committed
54

55
56
57
      ! Command line arguments
      character(len=128)   :: arg1, arg2, arg3, arg4
      integer :: mpierr
Andreas Marek's avatar
Andreas Marek committed
58

59
60
61
62
63
      ! default parameters
      na = 4000
      nev = 1500
      nblk = 16
      write_to_file = .false.
Andreas Marek's avatar
Andreas Marek committed
64

65
      if (.not. any(COMMAND_ARGUMENT_COUNT() == [0, 3, 4])) then
66
67
        write(error_unit, '(a,i0,a)') "Invalid number (", COMMAND_ARGUMENT_COUNT(), ") of command line arguments!"
        write(error_unit, *) "Expected: program [ [matrix_size num_eigenvalues block_size] ""output""]"
68
69
        stop 1
      endif
Andreas Marek's avatar
Andreas Marek committed
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89

      if (COMMAND_ARGUMENT_COUNT() == 3) then
        call GET_COMMAND_ARGUMENT(1, arg1)
        call GET_COMMAND_ARGUMENT(2, arg2)
        call GET_COMMAND_ARGUMENT(3, arg3)

        read(arg1, *) na
        read(arg2, *) nev
        read(arg3, *) nblk
      endif

      if (COMMAND_ARGUMENT_COUNT() == 4) then
        call GET_COMMAND_ARGUMENT(1, arg1)
        call GET_COMMAND_ARGUMENT(2, arg2)
        call GET_COMMAND_ARGUMENT(3, arg3)
        call GET_COMMAND_ARGUMENT(4, arg4)
        read(arg1, *) na
        read(arg2, *) nev
        read(arg3, *) nblk

90
91
92
        if (arg4 .eq. "output") then
          write_to_file = .true.
        else
93
          write(error_unit, *) "Invalid value for output flag! Must be ""output"" or omitted"
94
95
96
          stop 1
        endif

Andreas Marek's avatar
Andreas Marek committed
97
98
99
100
101
      endif

    end subroutine

end module