Skip to content

Instantly share code, notes, and snippets.

View pablosanjose's full-sized avatar

Pablo San-Jose pablosanjose

16 followers · 3 following
View GitHub Profile
@pablosanjose
pablosanjose / # mpich-scalapack - 2022-02-22_08-15-46.txt
Created February 22, 2022 15:56
mpich-scalapack (dpo/mumps-jl/mpich-scalapack) on macOS 12.2 - Homebrew build logs
Homebrew build logs for dpo/mumps-jl/mpich-scalapack on macOS 12.2
Build date: 2022-02-22 08:15:46
@pablosanjose
pablosanjose / openblas-exshift.patch
Last active January 30, 2021 16:12
Patch for OpenBLAS@0.3.13
diff --git a/lapack-netlib/SRC/chgeqz.f b/lapack-netlib/SRC/chgeqz.f
index 73d35621c..0351c79e1 100644
--- a/lapack-netlib/SRC/chgeqz.f
+++ b/lapack-netlib/SRC/chgeqz.f
@@ -320,12 +320,13 @@ SUBROUTINE CHGEQZ( JOB, COMPQ, COMPZ, N, ILO, IHI, H, LDH, T, LDT,
$ C, SAFMIN, TEMP, TEMP2, TEMPR, ULP
COMPLEX ABI22, AD11, AD12, AD21, AD22, CTEMP, CTEMP2,
$ CTEMP3, ESHIFT, RTDISC, S, SHIFT, SIGNBC, T1,
- $ U12, X
+ $ U12, X, ABI12, Y
@pablosanjose
pablosanjose / result.txt
Created January 26, 2021 12:13
test Reference Lapack results, by steelt
This file has been truncated, but you can view the full file.
STARTING ITERATION 1
NO CONVERGENCE IN 0 ITERATIONS
A
0.00E+00 0.00E+00 -6.59E-01 0.00E+00 4.89E-01 0.00E+00 -3.72E-01 0.00E+00 -1.23E-01 0.00E+00 -6.95E-01 0.00E+00 -7.04E-01 0.00E+00 2.62E+00 0.00E+00 -8.26E-01 0.00E+00 1.33E+00 0.00E+00 2.16E-01 0.00E+00 -2.48E+00 0.00E+00 1.52E-01 0.00E+00 1.55E+00 0.00E+00 4.76E-01 0.00E+00 -2.29E+00 0.00E+00
1.00E+00 0.00E+00 -4.60E-01 0.00E+00 1.47E+00 0.00E+00 -7.74E-01 0.00E+00 -2.48E-01 0.00E+00 -1.42E+00 0.00E+00 -1.45E+00 0.00E+00 5.39E+00 0.00E+00 -1.66E+00 0.00E+00 2.54E+00 0.00E+00 3.67E-01 0.00E+00 -4.21E+00 0.00E+00 2.58E-01 0.00E+00 2.64E+00 0.00E+00 8.09E-01 0.00E+00 -3.90E+00 0.00E+00
0.00E+00 0.00E+00 1.25E+00 0.00E+00 -6.12E-01 0.00E+00 5.87E-01 0.00E+00 -4.51E-01 0.00E+00 1.29E+00 0.00E+00 2.14E+00 0.00E+00 -6.70E+00 0.00E+00 1.48E+00 0.00E+00 -2.39E-01 0.00E+00 -1.95E-01 0.00E+00 2.22E+00 0.00E+00 -3.92E-01 0.00E+00 -3.96E+00 0.00E+00 5.28E-01 0.00E+00 -2
@pablosanjose
pablosanjose / testlapack.jl
Last active January 22, 2021 10:56
Testing the zggev function of reference LAPACK from Julia
const liblapack = "/usr/local/Cellar/lapack/3.9.0_1/lib/liblapacke.3.9.0.dylib"
import LinearAlgebra: LAPACK, BlasFloat, BlasInt, LAPACKException,
DimensionMismatch, SingularException, PosDefException, chkstride1, checksquare
using Base: require_one_based_indexing
function myggev!(jobvl::AbstractChar, jobvr::AbstractChar, A::AbstractMatrix{ComplexF64}, B::AbstractMatrix{ComplexF64})
require_one_based_indexing(A, B)
chkstride1(A, B)