Skip to content

Instantly share code, notes, and snippets.

Avatar

Pablo San-Jose pablosanjose

View GitHub Profile
@pablosanjose
pablosanjose / openblas-exshift.patch
Last active Jan 30, 2021
Patch for OpenBLAS@0.3.13
View openblas-exshift.patch
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 Jan 26, 2021
test Reference Lapack results, by steelt
View result.txt
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
@pablosanjose
pablosanjose / testlapack.jl
Last active Jan 22, 2021
Testing the zggev function of reference LAPACK from Julia
View testlapack.jl
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)