diff options
Diffstat (limited to 'gnu')
| -rw-r--r-- | gnu/packages/maths.scm | 169 | ||||
| -rw-r--r-- | gnu/packages/patches/mumps-build-parallelism.patch | 13 |
2 files changed, 182 insertions, 0 deletions
diff --git a/gnu/packages/maths.scm b/gnu/packages/maths.scm index e5da8c2ea2..a70e9dea03 100644 --- a/gnu/packages/maths.scm +++ b/gnu/packages/maths.scm @@ -32,6 +32,7 @@ #:use-module (guix download) #:use-module (guix svn-download) #:use-module (guix utils) + #:use-module (guix build utils) #:use-module (guix build-system cmake) #:use-module (guix build-system gnu) #:use-module (gnu packages algebra) @@ -628,6 +629,174 @@ scientific applications modeled by partial differential equations.") ,@(delete "--with-mpi=0" ,cf))))) (synopsis "Library to solve PDEs (with complex scalars and MPI support)"))) +(define-public mumps + (package + (name "mumps") + (version "5.0.0") + (source + (origin + (method url-fetch) + (uri (string-append "http://mumps.enseeiht.fr/MUMPS_" + version ".tar.gz")) + (sha256 + (base32 + "0690yp73sqk8zn2jnrzdr5swnjdyd7j0774s4xamjjwcxarw87hr")) + (patches (list (search-patch "mumps-build-parallelism.patch"))))) + (build-system gnu-build-system) + (inputs + `(("fortran" ,gfortran) + ;; These are required for linking against mumps, but we let the user + ;; declare the dependency. + ("blas" ,openblas) + ("metis" ,metis) + ("scotch" ,scotch))) + (arguments + `(#:modules ((ice-9 match) + (ice-9 popen) + (srfi srfi-1) + ,@%gnu-build-system-modules) + #:phases + (modify-phases %standard-phases + (replace + 'configure + (lambda* (#:key inputs #:allow-other-keys) + (call-with-output-file "Makefile.inc" + (lambda (port) + (format port " +PLAT = +LIBEXT = .a +OUTC = -o +OUTF = -o +RM = rm -f~:[ +CC = gcc +FC = gfortran +FL = gfortran +INCSEQ = -I$(topdir)/libseq +LIBSEQ = -L$(topdir)/libseq -lmpiseq +LIBSEQNEEDED = libseqneeded~; +CC = mpicc +FC = mpifort +FL = mpifort~] +AR = ar vr # rules require trailing space, ugh... +RANLIB = ranlib +LIBBLAS = -L~a -lopenblas~@[ +SCALAP = -L~a -lscalapack~] +LIBOTHERS = -pthread +CDEFS = -DAdd_ +PIC = -fPIC +OPTF = -O2 -DALLOW_NON_INIT $(PIC) +OPTL = -O2 $(PIC) +OPTC = -O2 $(PIC) +INCS = $(INCSEQ) +LIBS = $(SCALAP) $(LIBSEQ) +LPORDDIR = $(topdir)/PORD/lib +IPORD = -I$(topdir)/PORD/include +LPORD = -L$(LPORDDIR) -lpord +ORDERINGSF = -Dpord~@[ +METISDIR = ~a +IMETIS = -I$(METISDIR)/include +LMETIS = -L$(METISDIR)/lib -lmetis +ORDERINGSF += -Dmetis~]~@[~:{ +SCOTCHDIR = ~a +ISCOTCH = -I$(SCOTCHDIR)/include +LSCOTCH = -L$(SCOTCHDIR)/lib ~a-lesmumps -lscotch -lscotcherr +ORDERINGSF += ~a~}~] +ORDERINGSC = $(ORDERINGSF) +LORDERINGS = $(LPORD) $(LMETIS) $(LSCOTCH) +IORDERINGSF = $(ISCOTCH) +IORDERINGSC = $(IPORD) $(IMETIS) $(ISCOTCH)" + (assoc-ref inputs "mpi") + (assoc-ref inputs "blas") + (assoc-ref inputs "scalapack") + (assoc-ref inputs "metis") + (match (list (assoc-ref inputs "pt-scotch") + (assoc-ref inputs "scotch")) + ((#f #f) + #f) + ((#f scotch) + `((,scotch "" "-Dscotch"))) + ((ptscotch _) + `((,ptscotch + "-lptesmumps -lptscotch -lptscotcherr " + "-Dptscotch"))))))))) + (replace + 'build + ;; By default only the d-precision library is built. Make with "all" + ;; target so that all precision libraries and examples are built. + (lambda _ + (zero? (system* "make" "all" + (format #f "-j~a" (parallel-job-count)))))) + (replace + 'check + ;; Run the simple test drivers, which read test input from stdin: + ;; from the "real" input for the single- and double-precision + ;; testers, and from the "cmplx" input for complex-precision + ;; testers. The EXEC-PREFIX key is used by the mumps-openmpi + ;; package to prefix execution with "mpirun". + (lambda* (#:key (exec-prefix '()) #:allow-other-keys) + (with-directory-excursion "examples" + (every + (lambda (prec type) + (let ((tester (apply open-pipe* + `(,OPEN_WRITE + ,@exec-prefix + ,(string-append "./" prec + "simpletest")))) + (input (open-input-file + (string-append "input_simpletest_" type)))) + (begin + (dump-port input tester) + (close-port input) + (zero? (close-pipe tester))))) + '("s" "d" "c" "z") + '("real" "real" "cmplx" "cmplx"))))) + (replace + 'install + (lambda* (#:key outputs #:allow-other-keys) + (let ((out (assoc-ref outputs "out"))) + (copy-recursively "lib" (string-append out "/lib")) + (copy-recursively "include" (string-append out "/include")) + (when (file-exists? "libseq/libmpiseq.a") + (copy-file "libseq/libmpiseq.a" + (string-append out "/lib/libmpiseq.a"))))))))) + (home-page "http://mumps.enseeiht.fr") + (synopsis "Multifrontal sparse direct solver") + (description + "MUMPS (MUltifrontal Massively Parallel sparse direct Solver) solves a +sparse system of linear equations A x = b using Guassian elimination.") + (license license:cecill-c))) + +(define-public mumps-metis + (package (inherit mumps) + (name "mumps-metis") + (inputs + (alist-delete "scotch" (package-inputs mumps))))) + +(define-public mumps-openmpi + (package (inherit mumps) + (name "mumps-openmpi") + (inputs + `(("mpi" ,openmpi) + ("scalapack" ,scalapack) + ("pt-scotch" ,pt-scotch) + ,@(alist-delete "scotch" (package-inputs mumps)))) + (arguments + (substitute-keyword-arguments (package-arguments mumps) + ((#:phases phases) + `(modify-phases ,phases + (replace + 'check + (lambda _ + ((assoc-ref ,phases 'check) + #:exec-prefix '("mpirun" "-n" "2")))))))) + (synopsis "Multifrontal sparse direct solver (with MPI)"))) + +(define-public mumps-metis-openmpi + (package (inherit mumps-openmpi) + (name "mumps-metis-openmpi") + (inputs + (alist-delete "pt-scotch" (package-inputs mumps-openmpi))))) + (define-public superlu (package (name "superlu") diff --git a/gnu/packages/patches/mumps-build-parallelism.patch b/gnu/packages/patches/mumps-build-parallelism.patch new file mode 100644 index 0000000000..26dd330311 --- /dev/null +++ b/gnu/packages/patches/mumps-build-parallelism.patch @@ -0,0 +1,13 @@ +Building sublibraries in parallel can lead to race conditions on the libseq +and lipord targets. + +--- MUMPS_5.0.0/Makefile.orig 2015-06-15 10:08:54.523146562 -0500 ++++ MUMPS_5.0.0/Makefile 2015-06-15 10:19:44.074367512 -0500 +@@ -11,6 +11,7 @@ + sexamples dexamples cexamples zexamples \ + mumps_lib requiredobj libseqneeded clean + ++.NOTPARALLEL: + alllib: c z s d + all: cexamples zexamples sexamples dexamples + |