I noticed that dgeev returns real & imaginary parts of the complex eigenvectors -- reasonable, but unexpected. J9.04k Windows
load '~addons/math/lapack2/example/dgeev.ijs' ev 0.799482 _0.0994125j0.400792 _0.0994125j_0.400792 _0.100657 vr _0.655089 _0.193302 0.254632 0.125333 _0.523629 0.251857 _0.522405 0.332022 0.536218 0.0971825 _0.308384 0.593838 _0.0956068 0.675954 0 0.722087 'e1 e2 e3 e4'=:ev 'v1 v2 v3 v4'=:|:vr a 0.35 0.45 _0.14 _0.17 0.09 0.07 _0.54 0.35 _0.44 _0.33 _0.03 0.17 0.25 _0.32 _0.13 0.11 a +/ . * v1 _0.523732 _0.418632 0.428697 _0.0764359 e1*v1 _0.523732 _0.418632 0.428697 _0.0764359 a +/ . * v2 j. v3 _0.0828378j_0.102787 0.184338j0.152876 0.113937j0.0696072 _0.0671983j0.270917 e2 * v2 j. v3 _0.0828378j_0.102787 0.184338j0.152876 0.113937j0.0696072 _0.0671983j0.270917 On Wed, Feb 15, 2023 at 6:46 AM Raul Miller <rauldmil...@gmail.com> wrote: > For documentation, it looks like > https://code.jsoftware.com/wiki/Vocabulary/LAPACK > and maybe > https://code.jsoftware.com/wiki/Addons/math/lapack > > should be decent. > > As for the definitions themselves, it looks like > > require'math/lapack2' > > brings in all of the lapack2 the definitions, and that they're all in > the jlapack2 locale. > > So, dgeev_jlapack2_ should work. Or, if you prefer, you could use > > coinsert'jlapack2' > > after which dgeev should work. > > I hope this helps, > > -- > Raul > > On Wed, Feb 15, 2023 at 4:21 AM Erik Papp <erikmail...@gmail.com> wrote: > > > > Thank you for your answer! > > > > The example is running perfectly, but the eigenvalues, the eigenvectors > and > > SVD are not. > > I do not understand why? > > > > load'math/lapack2' > > > > load'math/lapack2/example/dgeev' > > > > 0.799482 _0.0994125j0.400792 _0.0994125j_0.400792 _0.100657 > > > > > ┌────────────────────────────────────────┬───────────────────────────────────────┐ > > > > │ _0.624471 0.533023 0 0.664141 │ _0.655089 _0.193302 0.254632 > > 0.125333 │ > > > > │ _0.599489 _0.266616 0.404136 _0.106815 │ _0.523629 0.251857 _0.522405 > > 0.332022 │ > > > > │ 0.499916 0.345526 0.315285 0.729325 │ 0.536218 0.0971825 _0.308384 > > 0.593838 │ > > > > │_0.0270862 _0.254081 _0.445113 0.124866 │_0.0956068 0.675954 0 > > 0.722087 │ > > > > > └────────────────────────────────────────┴───────────────────────────────────────┘ > > > > load'math/lapack2/example/dgesvd' > > > > > ┌──────────┬─────────────────────────────┬───────────────────────────────────────┐ > > > > │ 1.45189 │_0.819947 0.556294 0.134997 │_0.801549 _0.446592 _0.314262 > > _0.243545│ > > > > │ 0.143339 │_0.466218 _0.512125 _0.721366│ 0.572863 _0.391724 _0.512788 > > _0.505399│ > > > > │0.00422124│_0.332157 _0.65442 0.679269 │ 0.16912 _0.738115 0.120025 > > 0.642012 │ > > > > │ │ │0.0272991 _0.319835 0.789862 > > _0.522569 │ > > > > > └──────────┴─────────────────────────────┴───────────────────────────────────────┘ > > > > 1.45189 0 0 0 > > > > 0 0.143339 0 0 > > > > 0 0 0.00422124 0 > > > > ┌───────────────────────────┬───────────────────────────┐ > > > > │ 1 0.5 0.3333 0.25 │ 1 0.5 0.3333 0.25 │ > > > > │ 0.5 0.3333 0.25 0.2 │ 0.5 0.3333 0.25 0.2 │ > > > > │0.3333 0.25 0.2 0.1667 │0.3333 0.25 0.2 0.1667 │ > > > > └───────────────────────────┴───────────────────────────┘ > > > > coclass 'jlapack' > > > > > > > > dgeev=: geev > > > > tdgeev=: tgeev > > > > testdgeev=: testgeev > > > > > > > > require '~addons/math/lapack2/geev.ijs' > > > > not found: c:/users/tulajdonos/j903/addons/math/lapack2/geev.ijs > > > > |file name error: script > > > > | 0!:0 y[4!:55<'y' > > > > > > > > dgesvd=: gesvd > > > > tdgsvd=: tgsvd > > > > restdgsvd=: testgsvd > > > > > > > > require '~addons/math/lapack2/gesvd.ijs' > > > > not found: c:/users/tulajdonos/j903/addons/math/lapack2/gesvd.ijs > > > > |file name error: script > > > > | 0!:0 y[4!:55<'y' > > > > load'addons/math/lapack/gesvd.ijs' > > > > not found: C:/Users/Tulajdonos/addons/math/lapack/gesvd.ijs > > > > |file name error: script > > > > | 0!:0 y[4!:55<'y' > > > > > > > > bill lam <bbill....@gmail.com> ezt írta (időpont: 2023. febr. 14., K, > > 15:43): > > > > > The math/lapack addon was no longer supported on current J. > > > Instead, lapack is provided by math/lapack2. The windows binaries from > it > > > has optimized blas, > > > should be order of magnitude faster than that in the older lapack > addon. > > > However, usage of these 2 addons are different. > > > You need to study various scripts under the test folder. > > > The dgeev.ijs should be the script you need. > > > > > > FYI the left argument of dyad verb tdgeev there is the parameter for > data > > > type/precision. > > > 0 1 2 3 for double complex single-precision float and > single-precision > > > complex. > > > > > > > > > On Tue, Feb 14, 2023 at 10:06 PM Erik Papp <erikmail...@gmail.com> > wrote: > > > > > > > Dear J Jusers, > > > > > > > > Earlier I used *LAPACK* in the *J602a* and *J804* versions without > any > > > > problems. > > > > > > > > I installed the *J903* version along with *LAPACK* based on > > > > > > > > *System/Installation/J903/AIO*. I repeated the installation several > > > > times, but it does not run, > > > > > > > > it *always gives an error*. > > > > > > > > I tried running the *Addons/math/lapack* example, but I kept getting > > > > an error message. > > > > > > > > After > > > > > > > > load 'math/lapack math/lapack/dgeev' > > > > > > > > not found: c:/users/tulajdonos/j903/addons/math/lapack/lapack.ijs > > > > > > > > |file name error: script > > > > > > > > | 0!:0 y[4!:55<'y' > > > > > > > > I getting error message. > > > > > > > > It was the same problem after > > > > > > > > load 'math/lapack2 math/lapack2/dgeev' > > > > > > > > not found: c:/users/tulajdonos/j903/addons/math/lapack2/dgeev.ijs > > > > > > > > |file name error: script > > > > > > > > | 0!:0 y[4!:55<'y' > > > > > > > > > > > > Has anyone had a similar experience? > > > > > > > > How to install LAPACK on J903? > > > > > > > > Thank you in advance for your kind help. > > > > > > > > Sincerely, Erik > > > > > ---------------------------------------------------------------------- > > > > For information about J forums see > http://www.jsoftware.com/forums.htm > > > > > > > ---------------------------------------------------------------------- > > > For information about J forums see http://www.jsoftware.com/forums.htm > > > > > ---------------------------------------------------------------------- > > For information about J forums see http://www.jsoftware.com/forums.htm > ---------------------------------------------------------------------- > For information about J forums see http://www.jsoftware.com/forums.htm > ---------------------------------------------------------------------- For information about J forums see http://www.jsoftware.com/forums.htm
