This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Re: gsl-discuss Digest 22 Nov 2002 09:55:44 -0000 Issue 685
- From: Trond Mathias Døskeland <trond dot doskeland at nhh dot no>
- To: gsl-discuss-digest-help at sources dot redhat dot com
- Cc: gsl-discuss at sources dot redhat dot com
- Date: Sun, 24 Nov 2002 14:54:14 +0100
- Subject: Re: gsl-discuss Digest 22 Nov 2002 09:55:44 -0000 Issue 685
- Organization: NHH
- References: <1037958944.10682.ezmlm@sources.redhat.com>
I have the same problem but can't fix it. When I change the settings in
Project..Settings..C/C++...Project options from /MLd to /MD I get many new errors...
gsl-discuss-digest-help@sources.redhat.com wrote:
> gsl-discuss Digest 22 Nov 2002 09:55:44 -0000 Issue 685
>
> Topics (messages 2869 through 2876):
>
> GSL_Wigner_6J error
> 2869 by: Gerard Jungman <jungman@lanl.gov>
> 2870 by: Peter Roche <P.J.P.Roche@damtp.cam.ac.uk>
> 2871 by: "iwamae atsushi" <iwamaeatsushi@hotmail.com>
> 2872 by: "iwamae atsushi" <iwamaeatsushi@hotmail.com>
> 2874 by: Brian Gough <bjg@network-theory.co.uk>
>
> GSL Function Availablity
> 2873 by: "K.LakshmiNarayanan" <narayan@cat.ernet.in>
> 2875 by: Manoj Warrier <mow@ipp.mpg.de>
>
> HELP : Can't link with MSVC.
> 2876 by: Brian Gough <bjg@network-theory.co.uk>
>
> Administrivia:
>
> To subscribe to the digest, e-mail:
> gsl-discuss-digest-subscribe@sources.redhat.com
>
> To unsubscribe from the digest, e-mail:
> gsl-discuss-digest-unsubscribe@sources.redhat.com
>
> To post to the list, e-mail:
> gsl-discuss@sources.redhat.com
>
> ----------------------------------------------------------------------
>
> ------------------------------------------------------------------------
>
> Subject: Re: GSL_Wigner_6J error
> Date: 20 Nov 2002 12:43:14 -0700
> From: Gerard Jungman <jungman@lanl.gov>
> To: Peter Roche <P.J.P.Roche@damtp.cam.ac.uk>
> CC: iwamae@kues.kyoto-u.ac.jp, gsl-discuss@sources.redhat.com
>
> Hi all. Sorry for the confusion here. I implemented this code,
> but I never looked closely at it. Is there a consensus about
> how it should work? It is definitely true that, as implemented,
> the function returns
>
> / Ja Jb Je \
> \ Jd Jc Jf /
>
> I'm not sure where this convention comes from, if anywhere.
>
> Should I change the interface or just change the documentation?
>
> Thanks.
>
> --
> G. Jungman
>
> On Mon, 2002-11-18 at 11:07, Peter Roche wrote:
> > Cheers, that solved the problem.
> >
> > Peter
> >
> > *********************************************************
> >
> > Peter Roche
> > Department of Applied Mathematics and Theoretical Physics,
> > University of Cambridge, Silver Street, Cambridge, CB3 9EW
> >
> > email: p.j.p.roche@damtp.cam.ac.uk
> >
> > *********************************************************
> >
> > On Thu, 14 Nov 2002, iwamae atsushi wrote:
> >
> > > Dear Peter
> > >
> > > One of my friend have found the same problem on the 6j symbol in GSL.
> > > In order to get the correct value, the order of J value should be written
> > > as
> > > follows,
> > >
> > > gsl_sf_coupling_6j(Ja, Jb, Jc, Jd, Je, Jf)
> > >
> > > { Ja, Jb, Je }
> > > Jd, Jc, Jf
> > >
> > > The manual discription is not consistent.
> > > It says
> > > {Ja, Jb, Jc}
> > > Jd, Je, Jf
> > > , but this order gives wrong value as you noticed
> > >
> > >
> > >
> > >
> > > Atsushi Iwamae
> > >
> > >
> > > >From: Peter Roche <P.J.P.Roche@damtp.cam.ac.uk>
> > > >To: gsl-discuss@sources.redhat.com
> > > >Subject: GSL_Wigner_6J error
> > > >Date: Wed, 13 Nov 2002 19:36:40 +0000 (GMT)
> > > >
> > > >
> > > >Dear All,
> > > >
> > > >I hope this isn't a repeat of an already dealt with problem. I have been
> > > >trying to use the GSL Wigner 6j coefficients routines
> > > >(gsl_sf_coupling_6j(2*ja, 2*jb, 2*jc, 2*jd,2*je,2*jf)) and have found,
> > > >I think, that some of the returned values from the GSL
> > > >routine are not correct, even for small input values.
> > > >
> > > >I have included a table showing, in the first column, the input
> > > >parameter values, the angular mometa (which are multiplied by 2 in my
> > > >actual program), the second column shows the 'correct' values, obtained
> > > >from another program/published tables, and the third column shows the
> > > >return values from the GSL routine.
> > > >
> > > > Input Correct values GSL return
> > > > (011100) 0.5773502691896257 0.0000000000000000
> > > > (101101) 0.3333333333333335 0.0000000000000000
> > > > (121101) 0.3333333333333334 0.0000000000000000
> > > > (101101) 0.3333333333333335 0.0000000000000000
> > > > (022200) 0.4472135954999579 0.0000000000000000
> > > > (132201) 0.2581988897471612 0.0000000000000000
> > > > (242202) 0.2000000000000000 0.0000000000000000
> > > > (213302) 0.1690308509457034 0.0000000000000000
> > > > (233302) 0.1690308509457033 0.0000000000000000
> > > > (253302) 0.1690308509457033 0.0000000000000000
> > > >
> > > >
> > > >Has anyone encountered the same problems, or can explain the differences.
> > > >
> > > >Cheers,
> > > >Peter Roche
> > > >
> > > >*********************************************************
> > > >
> > > >Peter Roche
> > > >Department of Applied Mathematics and Theoretical Physics,
> > > >University of Cambridge, Silver Street, Cambridge, CB3 9EW
> > > >
> > > >email: p.j.p.roche@damtp.cam.ac.uk
> > > >
> > > >*********************************************************
> > >
> > >
> > > _________________________________________________________________
> > > $B%&%$%k%9%a!<%k!"LBOG%a!<%kBP:v$J$i(B MSN Hotmail http://www.hotmail.com/
> > >
> > >
> --
> Gerard Jungman <jungman@lanl.gov>
> Los Alamos National Laboratory
>
> ------------------------------------------------------------------------
>
> Subject: Re: GSL_Wigner_6J error
> Date: Wed, 20 Nov 2002 20:30:18 +0000 (GMT)
> From: Peter Roche <P.J.P.Roche@damtp.cam.ac.uk>
> To: Gerard Jungman <jungman@lanl.gov>
> CC: iwamae@kues.kyoto-u.ac.jp, <gsl-discuss@sources.redhat.com>
>
> I think the convention is like that which relates the Racah W
> coefficients and Wigner 6J symbols:
>
> W(Ja Jb Jc Jd; Je Jf) = pow(-1,(Ja+Jb+Jc+Jd)) / Ja Jb Je \
> \ Jd Jc Jf /
>
> I think the simplest way to solve the problem is change the manual entry,
> rather than the source code?
>
> Regards,
> Peter
>
> *********************************************************
>
> Peter Roche
> Department of Applied Mathematics and Theoretical Physics,
> University of Cambridge, Silver Street, Cambridge, CB3 9EW
>
> email: p.j.p.roche@damtp.cam.ac.uk
>
> *********************************************************
>
> On 20 Nov 2002, Gerard Jungman wrote:
>
> >
> > Hi all. Sorry for the confusion here. I implemented this code,
> > but I never looked closely at it. Is there a consensus about
> > how it should work? It is definitely true that, as implemented,
> > the function returns
> >
> > / Ja Jb Je \
> > \ Jd Jc Jf /
> >
> > I'm not sure where this convention comes from, if anywhere.
> >
> > Should I change the interface or just change the documentation?
> >
> > Thanks.
> >
> > --
> > G. Jungman
> >
> >
> > On Mon, 2002-11-18 at 11:07, Peter Roche wrote:
> > > Cheers, that solved the problem.
> > >
> > > Peter
> > >
> > > *********************************************************
> > >
> > > Peter Roche
> > > Department of Applied Mathematics and Theoretical Physics,
> > > University of Cambridge, Silver Street, Cambridge, CB3 9EW
> > >
> > > email: p.j.p.roche@damtp.cam.ac.uk
> > >
> > > *********************************************************
> > >
> > > On Thu, 14 Nov 2002, iwamae atsushi wrote:
> > >
> > > > Dear Peter
> > > >
> > > > One of my friend have found the same problem on the 6j symbol in GSL.
> > > > In order to get the correct value, the order of J value should be written
> > > > as
> > > > follows,
> > > >
> > > > gsl_sf_coupling_6j(Ja, Jb, Jc, Jd, Je, Jf)
> > > >
> > > > { Ja, Jb, Je }
> > > > Jd, Jc, Jf
> > > >
> > > > The manual discription is not consistent.
> > > > It says
> > > > {Ja, Jb, Jc}
> > > > Jd, Je, Jf
> > > > , but this order gives wrong value as you noticed
> > > >
> > > >
> > > >
> > > >
> > > > Atsushi Iwamae
> > > >
> > > >
> > > > >From: Peter Roche <P.J.P.Roche@damtp.cam.ac.uk>
> > > > >To: gsl-discuss@sources.redhat.com
> > > > >Subject: GSL_Wigner_6J error
> > > > >Date: Wed, 13 Nov 2002 19:36:40 +0000 (GMT)
> > > > >
> > > > >
> > > > >Dear All,
> > > > >
> > > > >I hope this isn't a repeat of an already dealt with problem. I have been
> > > > >trying to use the GSL Wigner 6j coefficients routines
> > > > >(gsl_sf_coupling_6j(2*ja, 2*jb, 2*jc, 2*jd,2*je,2*jf)) and have found,
> > > > >I think, that some of the returned values from the GSL
> > > > >routine are not correct, even for small input values.
> > > > >
> > > > >I have included a table showing, in the first column, the input
> > > > >parameter values, the angular mometa (which are multiplied by 2 in my
> > > > >actual program), the second column shows the 'correct' values, obtained
> > > > >from another program/published tables, and the third column shows the
> > > > >return values from the GSL routine.
> > > > >
> > > > > Input Correct values GSL return
> > > > > (011100) 0.5773502691896257 0.0000000000000000
> > > > > (101101) 0.3333333333333335 0.0000000000000000
> > > > > (121101) 0.3333333333333334 0.0000000000000000
> > > > > (101101) 0.3333333333333335 0.0000000000000000
> > > > > (022200) 0.4472135954999579 0.0000000000000000
> > > > > (132201) 0.2581988897471612 0.0000000000000000
> > > > > (242202) 0.2000000000000000 0.0000000000000000
> > > > > (213302) 0.1690308509457034 0.0000000000000000
> > > > > (233302) 0.1690308509457033 0.0000000000000000
> > > > > (253302) 0.1690308509457033 0.0000000000000000
> > > > >
> > > > >
> > > > >Has anyone encountered the same problems, or can explain the differences.
> > > > >
> > > > >Cheers,
> > > > >Peter Roche
> > > > >
> > > > >*********************************************************
> > > > >
> > > > >Peter Roche
> > > > >Department of Applied Mathematics and Theoretical Physics,
> > > > >University of Cambridge, Silver Street, Cambridge, CB3 9EW
> > > > >
> > > > >email: p.j.p.roche@damtp.cam.ac.uk
> > > > >
> > > > >*********************************************************
> > > >
> > > >
> > > > _________________________________________________________________
> > > > $B%&%$%k%9%a!<%k!"LBOG%a!<%kBP:v$J$i(B MSN Hotmail http://www.hotmail.com/
> > > >
> > > >
> > --
> > Gerard Jungman <jungman@lanl.gov>
> > Los Alamos National Laboratory
> >
> >
>
> ------------------------------------------------------------------------
>
> Subject: Re: GSL_Wigner_6J error
> Date: Thu, 21 Nov 2002 11:50:36 +0900
> From: "iwamae atsushi" <iwamaeatsushi@hotmail.com>
> Reply-To: iwamae@kues.kyoto-u.ac.jp
> To: gsl-discuss@sources.redhat.com
> CC: P.J.P.Roche@damtp.cam.ac.uk, jungman@lanl.gov
>
> Hi all
> In order to avoid the confusion of the order of 6js, the description in
> appendix C
> of "Me'canique Quantique" by Albert Messiah is helpful, I think.
>
> Each side of a quadrangle is correspond to the j values and the face two
> sides in a column of 6j-symbol.
>
> see Fig C.2 in Messiah
>
> { ja, jb, je }
> jd, jc, jf
>
> __jb__
> / \
> ja/ \jd
> /___________\
> jc
>
> diagonal lines
> ja + jb = je
> jb + jd = jf
>
> Atsushi Iwamae
> Department of Engineering Physics and Mechanics
> Kyoto University, 606-8501, Kyoto, Japan
>
> _________________________________________________________________
> $B:G?7$N%U%!%$%J%s%9>pJs$H%i%$%U%W%i%s$N%"%I%P%$%9(B MSN $B%^%M!<(B
> http://money.msn.co.jp/
>
> ------------------------------------------------------------------------
>
> Subject: Re: GSL_Wigner_6J error
> Date: Thu, 21 Nov 2002 12:42:01 +0900
> From: "iwamae atsushi" <iwamaeatsushi@hotmail.com>
> Reply-To: iwamae@kues.kyoto-u.ac.jp
> To: gsl-discuss@sources.redhat.com
>
> Hi all
> In order to avoid the confusion of the order of 6js, the description
> in appendix C of "Me'canique Quantique" by Albert Messiah is helpful.
>
> Each side of a tetragon in 3D, with four surfaces, is correspond to the j
> values and the face two sides in a column of 6j-symbol. (not reactangle in
> 2D, sorry for miss using the last mail)
>
> see Fig C.2 in Messiah
> (It is not easy to draw the tetrapack using asci code)
>
> { ja, jb, je }
> jd, jc, jf
>
> __jb__
> / \
> ja/ \jd
> /___________\
> jc
>
> ja + jb = je
> jb + jd = jf
>
> Atsushi Iwamae
> Department of Engineering Physics and Mechanics
> Kyoto University, 606-8501, Kyoto, Japan
>
> _________________________________________________________________
> $B%M%C%H$r;H$&$R$H$KM-Mx$JFCE5$$$C$Q$$(B MSN $B%+!<%I(B http://card.msn.co.jp/
>
> ------------------------------------------------------------------------
>
> Subject: Re: GSL_Wigner_6J error
> Date: Thu, 21 Nov 2002 11:57:06 +0000 (GMT)
> From: Brian Gough <bjg@network-theory.co.uk>
> To: Gerard Jungman <jungman@lanl.gov>
> CC: gsl-discuss@sources.redhat.com
>
> Gerard Jungman writes:
> >
> > Hi all. Sorry for the confusion here. I implemented this code,
> > but I never looked closely at it. Is there a consensus about
> > how it should work? It is definitely true that, as implemented,
> > the function returns
> >
> > / Ja Jb Je \
> > \ Jd Jc Jf /
> >
> > I'm not sure where this convention comes from, if anywhere.
> >
> > Should I change the interface or just change the documentation?
>
> We should use the most definitive version of 6j and 9j. Here
> is my suggestion:
>
> http://wwwinfo.cern.ch/asdoc/psdir/shortwrups.dir/u111.ps.gz
>
> ------------------------------------------------------------------------
>
> Subject: GSL Function Availablity
> Date: Thu, 21 Nov 2002 03:51:41 GMT
> From: "K.LakshmiNarayanan" <narayan@cat.ernet.in>
> To: gsl-discuss@sources.redhat.com
>
> Hello People,
>
> Does GSL has routines to ..
>
> (1)Calculate Eigen values of Ax=?Bx, Where A is a real symmetric matrix and
> B is a real symmetric positive-definite matrix.(i.e equal to NAG c function
> f02adc).
>
> 2)Calculate Eigen Vectors of Ax=?Bx, Where A is a real symmetric matrix and
> B is a real symmetric positive-definite matrix.(i.e equal to NAG c function
> f02aec).
>
> (3)Calculate all eigen values for real unsymmetric matrix.(i.e equal to NAG
> c Function f02afc).
>
> (4)Calculate all eigen vectors for real unsymmetric matrix.(i.e equal to NAG
> c Function f02agc).
>
> (5)Calculate all the eigenvalues and, if required, all the eigenvectors of
> the generalized eigen problem Ax = ?Bx where A and B are real, square
> matrices(i.e equal to NAG c Function f02bjc).
>
> (6)Computes selected eigenvalues and eigenvectors of a real general matrix.
> (i.e equal to NAG c Function f02ecc).
>
> (7)Computes selected eigenvalues and eigenvectors of a complex general
> matrix.(i.e equal to NAG c Function f02gcc).
>
> (8)Returns all, or part, of the singular value decomposition of a general
> real matrix.(i.e equal to NAG c Function f02wec).
>
> (9)Returns all, or part, of the singular value decomposition of a general
> complex matrix.(i.e equal to NAG c Function f02xec).
>
> Thanks in Advance.
>
> Narayan.
>
> ------------------------------------------------------------------------
>
> Subject: Re: GSL Function Availablity
> Date: Thu, 21 Nov 2002 13:35:46 +0100 (CET)
> From: Manoj Warrier <mow@ipp.mpg.de>
> To: "K.LakshmiNarayanan" <narayan@cat.ernet.in>
> CC: <gsl-discuss@sources.redhat.com>
>
> Lakshminarayanan,
>
> Check out http://sources.redhat.com/gsl/ref/gsl-ref_toc.html
> which has "table of contents" of the GSL manual. OR do U prefer
> someone to check it out for you if you please and let you know
> which routines you should use and also how one should call the
> routines?
>
> Sorry ppl. Im in a particularly nasty mood at the presumption some
> people make about the sociology of freely distributed software
> (more so if it is a fellow countryman). Hope I dont get kicked out
> of the mailing list for this :-|
>
> Manoj
>
> On Thu, 21 Nov 2002, K.LakshmiNarayanan wrote:
>
> > Hello People,
> >
> > Does GSL has routines to ..
> >
> > (1)Calculate Eigen values of Ax=?Bx, Where A is a real symmetric matrix and
> > B is a real symmetric positive-definite matrix.(i.e equal to NAG c function
> > f02adc).
> >
> > 2)Calculate Eigen Vectors of Ax=?Bx, Where A is a real symmetric matrix and
> > B is a real symmetric positive-definite matrix.(i.e equal to NAG c function
> > f02aec).
> >
> > (3)Calculate all eigen values for real unsymmetric matrix.(i.e equal to NAG
> > c Function f02afc).
> >
> > (4)Calculate all eigen vectors for real unsymmetric matrix.(i.e equal to NAG
> > c Function f02agc).
> >
> > (5)Calculate all the eigenvalues and, if required, all the eigenvectors of
> > the generalized eigen problem Ax = ?Bx where A and B are real, square
> > matrices(i.e equal to NAG c Function f02bjc).
> >
> > (6)Computes selected eigenvalues and eigenvectors of a real general matrix.
> > (i.e equal to NAG c Function f02ecc).
> >
> > (7)Computes selected eigenvalues and eigenvectors of a complex general
> > matrix.(i.e equal to NAG c Function f02gcc).
> >
> > (8)Returns all, or part, of the singular value decomposition of a general
> > real matrix.(i.e equal to NAG c Function f02wec).
> >
> > (9)Returns all, or part, of the singular value decomposition of a general
> > complex matrix.(i.e equal to NAG c Function f02xec).
> >
> >
> > Thanks in Advance.
> >
> > Narayan.
>
> --
> Manoj Warrier (manoj.warrier@ipp.mpg.de)
>
> Stellaratortheorie, Max-Planck Institut Fur Plasmaphysik
> TeilInstitut Greifswald Wendelsteinstrasse 1
> D-17491 Greifswald Germany Tel: +49-3834-882434
>
> --------- History of Computing 10-11-3003 ---------------
> Then there used to be this great user friendly OS which
> overwrote your MBR whenever you installed it.
> ---------------------------------------------------------
>
> ------------------------------------------------------------------------
>
> Subject: Re: HELP : Can't link with MSVC.
> Date: Fri, 22 Nov 2002 09:52:07 +0000 (GMT)
> From: Brian Gough <bjg@network-theory.co.uk>
> To: vincent daanen <vincent.daanen@imag.fr>
> CC: gsl-discuss@sources.redhat.com
>
> vincent daanen writes:
> > Hi,
> > I'm trying to use the simulated annealing optimization function.
> > When linking, the linker hangs with these errors :
> > Linking... MyRigidRegistrationDMap.obj : error LNK2001:
> > unresolved external
> > symbol _gsl_rng_default
> > Debug/SpinePer.exe : fatal error LNK1120: 1 unresolved
> > externals Error executing link.exe.
> >
> > SpinePer.exe - 2 error(s), 0 warning(s)
> > I add the gsl.lib and gslcblas.lib to the libraries to link with
> > but nothing change !
> >
> > does anybody know what should I do ??
> > PS : i'm using the binary build find at : www.network-theory.co.uk
> > <http://www.network-theory.co.uk/gsl/freedownloads.html>.
> >
>
> The library is a DLL so you need to compile your code with /MD.
>
> See the online help for the compiler about runtime link settings, or
> use the project file for the demo program as a starting point.