Hi, thank you very much . there is no windows version for gsl>1.8 . can u guide me how can i compile it from source on windows.
thanks & regards, Siji > Send Bug-gsl mailing list submissions to > [email protected] > > To subscribe or unsubscribe via the World Wide Web, visit > http://lists.gnu.org/mailman/listinfo/bug-gsl > or, via email, send a message with subject or body 'help' to > [email protected] > > You can reach the person managing the list at > [email protected] > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of Bug-gsl digest..." > > > Today's Topics: > > 1. gsl_eigen_nonsymm_workspace (Siji P George) > 2. Re: gsl_eigen_nonsymm_workspace (Patrick Alken) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Wed, 22 Jul 2009 17:10:32 +0530 (IST) > From: "Siji P George" <[email protected]> > Subject: [Bug-gsl] gsl_eigen_nonsymm_workspace > To: [email protected] > Message-ID: <[email protected]> > Content-Type: text/plain;charset=iso-8859-1 > > Hi, > > We are using gsl in our program and the version installed is gsl-1.8.exe. > > it is giving compilation error ,which is given below: > > : error C2065: 'gsl_eigen_nonsymm_workspace' : undeclared identifier > : error C2065: 'workspace' : undeclared identifier > : error C3861: 'gsl_eigen_nonsymm_alloc': identifier not found > : error C3861: 'gsl_eigen_nonsymm': identifier not found > : error C3861: 'gsl_eigen_nonsymm_free': identifier not found > > i have included the following headers: > > #include <gsl/gsl_errno.h> > #include <gsl/gsl_vector.h> > #include <gsl/gsl_matrix.h> > #include <gsl/gsl_blas.h> > #include <gsl/gsl_linalg.h> > #include <gsl/gsl_multiroots.h> > #include <gsl/gsl_sort.h> > #include <gsl/gsl_sort_vector.h> > #include <gsl/gsl_eigen.h> > > and the piece of code is also given below: > ------------------------------------------------------------------------------- > void SteadyState::classifyState( const double* T ) > { > unsigned int nConsv = nVarMols_ - rank_; > gsl_matrix* J = gsl_matrix_calloc ( nVarMols_, nVarMols_ ); > // double* yprime = new double[ nVarMols_ ]; > // vector< double > yprime( nVarMols_, 0.0 ); > // Generate an approximation to the Jacobean by generating small > // increments to each of the molecules in the steady state, one > // at a time, and putting the resultant rate vector into a column > // of the J matrix. > // This needs a bit of heuristic to decide what is a 'small' increment. > // Use the totals for this. > double tot = 0.0; > for ( unsigned int i = 0; i < nConsv; ++i ) { > tot += T[i]; > } > tot *= DELTA; > > // Fill up Jacobian > for ( unsigned int i = 0; i < nVarMols_; ++i ) { > double orig = s_->S()[i]; > s_->S()[i] = orig + tot; > s_->updateV(); > s_->S()[i] = orig; > // yprime.assign( nVarMols_, 0.0 ) > // vector< double >::iterator y = yprime.begin(); > > // Compute the rates for each mol. > for ( unsigned int j = 0; j < nVarMols_; ++j ) { > // *y++ = N_.computeRowRate( j, s_->velocity() ); > double temp = s_->N().computeRowRate( j, s_->velocity() > ); > gsl_matrix_set( J, i, j, temp ); > } > } > > // Jacobian is now ready. Find eigenvalues. > gsl_vector_complex* vec = gsl_vector_complex_alloc( nVarMols_ ); > gsl_eigen_nonsymm_workspace* workspace = > gsl_eigen_nonsymm_alloc( nVarMols_ ); > int status = gsl_eigen_nonsymm( J, vec, workspace ); > if ( status != GSL_SUCCESS ) { > cout << "Warning: failed to find eigenvalues. Status = " << > status << endl; > solutionStatus_ = 2; // Steady state OK, eig classification > failed > } else { // Eigenvalues are ready. Classify state. > nNegEigenvalues_ = 0; > nPosEigenvalues_ = 0; > for ( unsigned int i = 0; i < nVarMols_; ++i ) { > gsl_complex z = gsl_vector_complex_get( vec, i ); > double r = GSL_REAL( z ); > nNegEigenvalues_ += ( r < -EPSILON ); > nPosEigenvalues_ += ( r > EPSILON ); > } > if ( nNegEigenvalues_ == rank_ ) > stateType_ = 0; // Stable > else if ( nPosEigenvalues_ == rank_ ) > stateType_ = 1; // Unstable > else if (nNegEigenvalues_ == rank_ - 1) > stateType_ = 2; // Saddle > else > stateType_ = 3; // Other > } > > gsl_vector_complex_free( vec ); > gsl_matrix_free ( J ); > gsl_eigen_nonsymm_free( workspace ); > } > ---------------------------------------------------------------------------------- > > the same program works fine on linux. no idea what is wrong on windows? > plz help. > > > thanks & regards, > Siji > > > > > > ------------------------------ > > Message: 2 > Date: Wed, 22 Jul 2009 07:44:06 -0600 > From: Patrick Alken <[email protected]> > Subject: Re: [Bug-gsl] gsl_eigen_nonsymm_workspace > To: [email protected] > Message-ID: <[email protected]> > Content-Type: text/plain; charset=us-ascii > >> >> > We are using gsl in our program and the version installed is gsl-1.8.exe. > > Nonsymmetric eigensystems were added in gsl 1.9 > > > > > ------------------------------ > > _______________________________________________ > Bug-gsl mailing list > [email protected] > http://lists.gnu.org/mailman/listinfo/bug-gsl > > > End of Bug-gsl Digest, Vol 78, Issue 12 > *************************************** > _______________________________________________ Bug-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/bug-gsl
