[Jprogramming] Lapack on 64-bit Macs

Tue Aug 20 07:30:31 UTC 2013

Thanks. Updated again. 

Please run the test script test0.ijs and report if it run ok on a 
64-bit Mac.

Вт, 20 авг 2013, Marc Simpson писал(а):
> Thanks Bill.
> 
> Line 268 of lapack.ijs still needs updating though.
> 
> On Tue, Aug 20, 2013 at 2:42 AM, bill lam <bbill.lam at gmail.com> wrote:
> > Thanks for info. JAL updated.
> >
> > Пн, 19 авг 2013, Marc Simpson писал(а):
> >> Lapack is disabled on 64-bit Macs due to the following guard in
> >> 'math/lapack.ijs' (added in r4400):
> >>
> >>   if. IF64 > UNAME-: 'Linux' do.
> >>     '64-bit not supported' 13!:8[10
> >>   end.
> >>
> >> This signal is overly general as the OS X VecLib framework provides a
> >> 64-bit implementation. Since the 64-bit check is only relevant on
> >> Windows, here's a possible fix:
> >>
> >>   3 : 0''
> >>   if. UNAME-:'Linux' do.
> >>     if. IF64 do.
> >>       dll=: '"',path,'lapack64.so" '
> >>     else.
> >>       dll=: '"',path,'lapack.so" '
> >>     end.
> >>   elseif. UNAME-:'Darwin' do.
> >>     dll=: '/System/Library/Frameworks/vecLib.framework/vecLib '
> >>   elseif. UNAME-:'Win' do.
> >>     if. IF64 do.
> >>       '64-bit not supported' 13!:8[10
> >>     end.
> >>     dll=: '"',path,'jlapack.dll" '
> >>   elseif. do.
> >>     'platform not supported' 13!:8[10
> >>   end.
> >>   ''
> >>   )
> >>
> >> Best,
> >> Marc
> >> ----------------------------------------------------------------------
> >> For information about J forums see http://www.jsoftware.com/forums.htm
> >
> > --
> > regards,
> > ====================================================
> > GPG key 1024D/4434BAB3 2008-08-24
> > gpg --keyserver subkeys.pgp.net --recv-keys 4434BAB3
> > ----------------------------------------------------------------------
> > For information about J forums see http://www.jsoftware.com/forums.htm
> ----------------------------------------------------------------------
> For information about J forums see http://www.jsoftware.com/forums.htm

-- 
regards,
====================================================
GPG key 1024D/4434BAB3 2008-08-24
gpg --keyserver subkeys.pgp.net --recv-keys 4434BAB3
-------------- next part --------------
NB. examples

dbg 1
NB. dbstops'need'
dbstops''

mp=: +/ .*

need_jlapack_ 'gebal gees geev gehrd gels gelss geqrf gerqf gesv gesvd gesvx getrf heev lartg potrf trtrs'

gebal=: gebal_jlapack_
gees=: gees_jlapack_
geev=: geev_jlapack_
gehrd=: gehrd_jlapack_
gels=: gels_jlapack_
gelss=: gelss_jlapack_
geqrf=: geqrf_jlapack_
gerqf=: gerqf_jlapack_
gesv=: gesv_jlapack_
gesvd=: gesvd_jlapack_
gesvx=: gesvx_jlapack_
getrf=: getrf_jlapack_
heev=: heev_jlapack_
lartg=: lartg_jlapack_
potrf=: potrf_jlapack_
trtrs=: trtrs_jlapack_

m=: ?.4 4$100
m1=: ?.4 3$30
m2=: utri_jlapack_ m1
m3=: ?.4 5$30
mz=: 1j3+?.4 4$100
m1z=: 1j3+?.4 3$100
m2z=: utri_jlapack_ m1z
m3z=: 1j3+?.4 5$100
v=: ?.4 4$100
vz=: 1j3+?.4$100
m4=: 7 7 $ 6 0 0 0 0 1 0 0 4 0 0.00025 0.0125 0.02 0.125 1 128 64 0 0 _2 16 0 16384 0 1 _400 256 _4000 _2 _256 0 0.0125 2 2 32 0 0 0 0 0 0 0 0 8 0 0.004 0.125 _0.2 3
m4z=: 6 6 $ 1j1 1j1 0 1j1 1j1 1j1 1j1 1j1 0 1j1 1j1 1j1 1j1 1j1 1j1 1j1 1j1 1j1 0 0 0 1j1 0 0 1j1 1j1 0 1j1 1j1 1j1 1j1 1j1 0 1j1 1j1 1j1
v2=: _2 2
v2z=: _2j_2 2j_2

smoutput gebal m4             NB. balance a general square double matrix
smoutput gebal m4z            NB. balance a general square complex matrix
smoutput gees m               NB. real Schur factorization of a general double matrix
smoutput gees mz              NB. Schur factorization of a general complex matrix
smoutput geev m               NB. eigenvalues of a general double matrix
smoutput geev mz              NB. eigenvalues of a general complex matrix
smoutput gehrd m;1;(#m)       NB. reduce a general square double matrix to upper Hessenberg form
smoutput gehrd mz;1;(#mz)     NB. reduce a general square complex matrix to upper Hessenberg form
smoutput gels m3;v            NB. solves underdetermined double system m3 * x = v
smoutput gels m3z;vz          NB. solves underdetermined complex system m3z * xz = vz
smoutput gelss m3;v           NB. solves underdetermined double system m3 * x = v via SVD
smoutput gelss m3z;vz         NB. solves underdetermined complex system m3z * xz = vz via SVD
smoutput geqrf m1             NB. QR factorization of a general double matrix
smoutput geqrf m1z            NB. QR factorization of a general complex matrix
smoutput gerqf m1             NB. RQ factorization of a general double matrix
smoutput gerqf m1z            NB. RQ factorization of a general complex matrix
smoutput gesv m;m1            NB. solves double system m * x = m1
smoutput gesv mz;m1z          NB. solves complex system mz * xz = m1z
smoutput gesvd m              NB. SVD of a general double matrix
smoutput gesvd mz             NB. SVD of a general complex matrix
smoutput gesvx m;m1           NB. solves double system m * x = m1
smoutput gesvx mz;m1z         NB. solves complex system mz * xz = m1z
smoutput getrf m              NB. LU decomposition of a general double matrix
smoutput getrf mz             NB. LU decomposition of a general complex matrix
NB. the following test failed
NB. smoutput heev m               NB. eigenvalues and eigenvectors of a symmetric double matrix
NB. smoutput heev mz              NB. eigenvalues and eigenvectors of a Hermitian complex matrix
smoutput lartg v2             NB. generates a plane rotation of a double 2-vector
smoutput lartg v2z            NB. generates a plane rotation of a complex 2-vector
smoutput potrf m mp |: m      NB. Cholesky decomposition of symmetric positive definite double matrix
smoutput potrf mz mp +|: mz   NB. Cholesky decomposition of Hermitian positive definite complex matrix
smoutput trtrs m;v            NB. solves double triangular system m * x = v
smoutput trtrs mz;vz          NB. solves complex triangular system mz * xz = vz