473,320 Members | 1,535 Online
Bytes | Software Development & Data Engineering Community
Post Job

Home Posts Topics Members FAQ

Join Bytes to post your question to a community of 473,320 software developers and data experts.

solve linear systems in boost::numerics::ublas

Hi.

I'm using boost::numerics::ublas to solve linear symmetric and
hermitian systems (cholesky, qr, etc.). I'm aware that ublas only provide
basic algebra operations (blas1, blas2 and blas3) so i've tried some
other libraries to work with ublas, with no success so far. In order to
do so i downloaded boost::numerics::bindings (which aren't available in
www.boost.org right now but in the list which boost developers have in
yahoo), those bindings enable you to use libraries developed in c or
fortran to do more advanced algebra operations, but i'm stuck:

+ cLAPACK, when i tried to compile cLAPACK (a c library for which
there are bindings in boost), the build procedure in windows is quite
tricky and seems to be out of date. I was unable to build the library and
the prebuild libraries are outdated (1999).

+ Atlas, any of the prebuild libraries in sourceforge suits my
microprocessor and i can't build them since i don't have cygwin (i have a
56k. modem and i can't download such a big package), evenmore, i intend
to use the libraries in msvc and the object format differs. There's an
alternative procedure to build in mscv, but i think the 'nmake' tool in
msvc doesn't behave as the 'make' in the -ix world. ¿am i wrong?.

+ Bayes++, this library is intended to be used to calculate karman
filters, but they provide uLAPACK compliance, when i compile it's ok, but
when i try to link them, the functions which provide LAPACK functionality
are missing (in namespace bayesian_filters::LAPACK::rawLAPACK, dgetrs_,
dgetrf_, etc.)

I've also tried some other libraries like blitz++ and tnt (formerly
lapack++), but i didn't suceed in building the libraries.

So, i'm really stuck, i think i should look for prebuild libraries
at this point, but my micro is relatively odd (amd duron). Please, ¿what
step would you take next to perform cholesky factorization?. I don't have
the expertise nor the time to write the algorithms hand-craft

Any help would be appreciated, thank you very much in advance and
excuse my bad english.
Jul 22 '05 #1
3 4049
Manzanita wrote:
Hi.

I'm using boost::numerics::ublas to solve linear symmetric and
hermitian systems (cholesky, qr, etc.). I'm aware that ublas only provide
basic algebra operations (blas1, blas2 and blas3) so i've tried some
other libraries to work with ublas, with no success so far. In order to
do so i downloaded boost::numerics::bindings (which aren't available in
www.boost.org right now but in the list which boost developers have in
yahoo), those bindings enable you to use libraries developed in c or
fortran to do more advanced algebra operations, but i'm stuck:

+ cLAPACK, when i tried to compile cLAPACK (a c library for which
there are bindings in boost), the build procedure in windows is quite
tricky and seems to be out of date. I was unable to build the library and
the prebuild libraries are outdated (1999).

+ Atlas, any of the prebuild libraries in sourceforge suits my
microprocessor and i can't build them since i don't have cygwin (i have a
56k. modem and i can't download such a big package), evenmore, i intend
to use the libraries in msvc and the object format differs. There's an
alternative procedure to build in mscv, but i think the 'nmake' tool in
msvc doesn't behave as the 'make' in the -ix world. ¿am i wrong?.

+ Bayes++, this library is intended to be used to calculate karman
filters, but they provide uLAPACK compliance, when i compile it's ok, but
when i try to link them, the functions which provide LAPACK functionality
are missing (in namespace bayesian_filters::LAPACK::rawLAPACK, dgetrs_,
dgetrf_, etc.)

I've also tried some other libraries like blitz++ and tnt (formerly
lapack++), but i didn't suceed in building the libraries.

So, i'm really stuck, i think i should look for prebuild libraries
at this point, but my micro is relatively odd (amd duron). Please, ¿what
step would you take next to perform cholesky factorization?. I don't have
the expertise nor the time to write the algorithms hand-craft


You should persist in trying to get [C]LAPACK to build, independently of
other libraries (ublas, etc.). [C]LAPACK does depend on [C]BLAS, but a
version of the latter is usually included in the former (see SRC/BLAS
directory). ATLAS is a replacement for [C]BLAS - you may improve
performance by linking with it, but it's not necessary.

If you use the Fortran versions, you will also need G2C.

See README and INSTALL files in LAPACK for info about exactly how to
install.

Cheers,
Oleg

--
Oleg Trott <ol********@columbia.edu>
Jul 22 '05 #2
Oleg Trott <ol********@columbia.edu> wrote in message news:<c1*********@newsmaster.cc.columbia.edu>...
Manzanita wrote:
Hi.

I'm using boost::numerics::ublas to solve linear symmetric and
hermitian systems (cholesky, qr, etc.). I'm aware that ublas only provide
basic algebra operations (blas1, blas2 and blas3) so i've tried some
other libraries to work with ublas, with no success so far. In order to
do so i downloaded boost::numerics::bindings (which aren't available in
www.boost.org right now but in the list which boost developers have in
yahoo), those bindings enable you to use libraries developed in c or
fortran to do more advanced algebra operations, but i'm stuck:

+ cLAPACK, when i tried to compile cLAPACK (a c library for which
there are bindings in boost), the build procedure in windows is quite
tricky and seems to be out of date. I was unable to build the library and
the prebuild libraries are outdated (1999).

+ Atlas, any of the prebuild libraries in sourceforge suits my
microprocessor and i can't build them since i don't have cygwin (i have a
56k. modem and i can't download such a big package), evenmore, i intend
to use the libraries in msvc and the object format differs. There's an
alternative procedure to build in mscv, but i think the 'nmake' tool in
msvc doesn't behave as the 'make' in the -ix world. ¿am i wrong?.

+ Bayes++, this library is intended to be used to calculate karman
filters, but they provide uLAPACK compliance, when i compile it's ok, but
when i try to link them, the functions which provide LAPACK functionality
are missing (in namespace bayesian_filters::LAPACK::rawLAPACK, dgetrs_,
dgetrf_, etc.)

I've also tried some other libraries like blitz++ and tnt (formerly
lapack++), but i didn't suceed in building the libraries.

So, i'm really stuck, i think i should look for prebuild libraries
at this point, but my micro is relatively odd (amd duron). Please, ¿what
step would you take next to perform cholesky factorization?. I don't have
the expertise nor the time to write the algorithms hand-craft
You should persist in trying to get [C]LAPACK to build, independently of
other libraries (ublas, etc.). [C]LAPACK does depend on [C]BLAS, but a
version of the latter is usually included in the former (see SRC/BLAS
directory). ATLAS is a replacement for [C]BLAS - you may improve
performance by linking with it, but it's not necessary.

If you use the Fortran versions, you will also need G2C.


Or you could just call the Fortran libraries directly from
C/C++...just remember in Fortran, everything is passed by pointer. I
have a class library that interfaces to LAPACK that I can show you as
an example.

Tino

See README and INSTALL files in LAPACK for info about exactly how to
install.

Cheers,
Oleg

Jul 22 '05 #3
Tino wrote:
Oleg Trott wrote:

If you use the Fortran versions, you will also need G2C.


Or you could just call the Fortran libraries directly from
C/C++...just remember in Fortran, everything is passed by pointer. I
have a class library that interfaces to LAPACK that I can show you as
an example.


It probably depends on the compiler actually, but none of the compilers I
use (GCC and Intel) will link to LAPACK without G2C. I think the reason is
that some of the intrinsic Fortran functions (like lg10) do not exist in C,
and so they have to be provided in libraries.

Cheers,
Oleg

--
Oleg Trott <ol********@columbia.edu>
Jul 22 '05 #4

This thread has been closed and replies have been disabled. Please start a new discussion.

Similar topics

15
by: christopher diggins | last post by:
Here is some code I wrote for Matrix multiplication for arbitrary dimensionality known at compile-time. I am curious how practical it is. For instance, is it common to know the dimensionality of...
0
by: ryjfgjl | last post by:
ExcelToDatabase: batch import excel into database automatically...
0
isladogs
by: isladogs | last post by:
The next Access Europe meeting will be on Wednesday 6 Mar 2024 starting at 18:00 UK time (6PM UTC) and finishing at about 19:15 (7.15PM). In this month's session, we are pleased to welcome back...
1
isladogs
by: isladogs | last post by:
The next Access Europe meeting will be on Wednesday 6 Mar 2024 starting at 18:00 UK time (6PM UTC) and finishing at about 19:15 (7.15PM). In this month's session, we are pleased to welcome back...
0
by: Vimpel783 | last post by:
Hello! Guys, I found this code on the Internet, but I need to modify it a little. It works well, the problem is this: Data is sent from only one cell, in this case B5, but it is necessary that data...
0
by: jfyes | last post by:
As a hardware engineer, after seeing that CEIWEI recently released a new tool for Modbus RTU Over TCP/UDP filtering and monitoring, I actively went to its official website to take a look. It turned...
1
by: PapaRatzi | last post by:
Hello, I am teaching myself MS Access forms design and Visual Basic. I've created a table to capture a list of Top 30 singles and forms to capture new entries. The final step is a form (unbound)...
1
by: CloudSolutions | last post by:
Introduction: For many beginners and individual users, requiring a credit card and email registration may pose a barrier when starting to use cloud servers. However, some cloud server providers now...
1
by: Defcon1945 | last post by:
I'm trying to learn Python using Pycharm but import shutil doesn't work
1
by: Shællîpôpï 09 | last post by:
If u are using a keypad phone, how do u turn on JavaScript, to access features like WhatsApp, Facebook, Instagram....

By using Bytes.com and it's services, you agree to our Privacy Policy and Terms of Use.

To disable or enable advertisements and analytics tracking please visit the manage ads & tracking page.