benchmark: Add netlib-hpl.

* guix-science/packages/benchmark.scm (netlib-hpl): New variable.

2 files changed, 71 insertions, 0 deletions

diff --git a/guix-science/packages/benchmark.scm b/guix-science/packages/benchmark.scm
index f0d7f26..4242c8e 100644
--- a/guix-science/packages/benchmark.scm
+++ b/guix-science/packages/benchmark.scm
@@ -7,7 +7,9 @@
   #:use-module ((guix licenses)
                 #:prefix license:)
   #:use-module (guix build-system gnu)
+  #:use-module (gnu packages)
   #:use-module (gnu packages gcc)
+  #:use-module (gnu packages maths)
   #:use-module (gnu packages mpi)
   #:use-module (gnu packages perl)
   #:use-module (ice-9 match))
@@ -208,3 +210,52 @@ indicated as different classes.")
                   "NPB" #$(version-major+minor (package-version this-package))
                   "-MPI"))))))))
     (synopsis "NAS Parallel Benchmarks (NPB), MPI variant")))
+
+(define-public netlib-hpl
+  (package
+    (name "netlib-hpl")
+    (version "2.3")
+    (source
+     (origin
+       (method url-fetch)
+       (uri (string-append "https://netlib.org/benchmark/hpl/hpl-" version
+                           ".tar.gz"))
+       (sha256
+        (base32 "0c18c7fzlqxifz1bf3izil0bczv3a7nsv0dn6winy3ik49yw3i9j"))
+       (patches (search-patches
+                 "patches/netlib-hpl-use-default-hpl-dat.patch"))))
+    (build-system gnu-build-system)
+    (inputs (list openblas openmpi))
+    (arguments
+     (list
+      #:configure-flags
+      #~(list "CFLAGS=-DHPL_DETAILED_TIMING -DHPL_CALL_CBLAS")
+      #:phases
+      #~(modify-phases %standard-phases
+          (add-after 'unpack 'set-hpl-dat-path
+            (lambda _
+              (substitute* "testing/ptest/HPL_pdinfo.c"
+                (("@@GUIX_HPL_DAT@@")
+                 (string-append #$output
+                                "/share/"
+                                #$(package-name this-package)
+                                "-"
+                                #$(package-version this-package)
+                                "/examples")))))
+          (add-after 'install 'install-hpl-dat
+            (lambda _
+              (install-file "testing/ptest/HPL.dat"
+                            (string-append #$output
+                                           "/share/"
+                                           #$(package-name this-package)
+                                           "-"
+                                           #$(package-version this-package)
+                                           "/examples")))))))
+    (home-page "https://netlib.org/benchmark/hpl/")
+    (synopsis "High-Performance Linpack Benchmark for Distributed-Memory Computers")
+    (description "HPL is a software package that solves a (random) dense linear system
+in double precision (64 bits) arithmetic on distributed-memory
+computers.  It can thus be regarded as a portable as well as freely
+available implementation of the High Performance Computing Linpack
+Benchmark.")
+    (license (license:non-copyleft "https://netlib.org/benchmark/hpl/copyright.html"))))
diff --git a/patches/netlib-hpl-use-default-hpl-dat.patch b/patches/netlib-hpl-use-default-hpl-dat.patch
new file mode 100644
index 0000000..716eca3
--- /dev/null
+++ b/patches/netlib-hpl-use-default-hpl-dat.patch
@@ -0,0 +1,20 @@
+# Use example file for benchmark.
+--- hpl-2.3/testing/ptest/HPL_pdinfo.c  2024-12-20 14:26:45.296458712 +0100
++++ hpl-2.3/testing/ptest/HPL_pdinfo.c  2024-12-20 14:41:54.994292754 +0100
+@@ -297,8 +297,13 @@
+       if( ( infp = fopen( "HPL.dat", "r" ) ) == NULL )
+       { 
+-         HPL_pwarn( stderr, __LINE__, "HPL_pdinfo",
+-                    "cannot open file HPL.dat" );
+-         error = 1; goto label_error;
++         fprintf( stderr, "no file HPL.dat in current working directory\n"
++                          "using @@GUIX_HPL_DAT@@/HPL.dat\n\n");
++         if( ( infp = fopen( "@@GUIX_HPL_DAT@@/HPL.dat", "r" ) ) == NULL )
++         {
++             HPL_pwarn( stderr, __LINE__, "HPL_pdinfo",
++                        "cannot open file @@GUIX_HPL_DAT@@/HPL.dat");
++             error = 1; goto label_error;
++         }
+       }
+ 
+       (void) fgets( line, HPL_LINE_MAX - 2, infp );