[Csdp] Using cblas and clapack with CSDP
Henri Kjellberg
henrikjellberg at gmail.com
Sun Nov 7 21:17:57 EST 2010
Hi,
I have been able to get the example code to execute however it seems to not
be perfectly happy here is the output when I tried to execute ./example
root at ubuntu:/opt/Csdp-6.1.1/example# ./example
Iter: 0 Ap: 0.00e+00 Pobj: 3.0434379e+02 Ad: 0.00e+00 Dobj:
0.0000000e+00
Iter: 1 Ap: 5.57e-01 Pobj: 1.2392573e+02 Ad: 8.57e-01 Dobj:
3.3873125e+01
Iter: 2 Ap: 5.91e-01 Pobj: 4.4688559e+01 Ad: 1.00e+00 Dobj:
7.2128795e+01
Iter: 3 Ap: 1.65e-02 Pobj: 4.3979696e+01 Ad: 1.00e+00 Dobj:
9.0347493e+01
Iter: 4 Ap: 1.57e-02 Pobj: 4.3315315e+01 Ad: 1.00e+00 Dobj:
1.1978791e+02
Iter: 5 Ap: 4.98e-01 Pobj: 2.2831612e+01 Ad: 1.00e+00 Dobj:
2.1875630e+02
Iter: 6 Ap: 6.91e-01 Pobj: 8.0668171e+00 Ad: 1.00e+00 Dobj:
2.8249267e+02
Iter: 7 Ap: 2.00e-02 Pobj: 7.9342969e+00 Ad: 1.00e+00 Dobj:
2.6557930e+02
Iter: 8 Ap: 1.92e-02 Pobj: 7.8099083e+00 Ad: 1.00e+00 Dobj:
2.5983512e+02
Iter: 9 Ap: 7.21e-01 Pobj: 3.2291368e+00 Ad: 1.00e+00 Dobj:
2.5117984e+02
Iter: 10 Ap: 9.50e-01 Pobj: 1.5471455e+00 Ad: 1.00e+00 Dobj:
1.3933178e+02
Iter: 11 Ap: 1.00e+00 Pobj: 1.4724451e+00 Ad: 9.68e-01 Dobj:
4.6961726e+01
Iter: 12 Ap: 1.00e+00 Pobj: 1.5142991e+00 Ad: 1.00e+00 Dobj:
1.6000624e+01
Iter: 13 Ap: 1.00e+00 Pobj: 1.6405417e+00 Ad: 1.00e+00 Dobj:
6.8112636e+00
Iter: 14 Ap: 9.00e-01 Pobj: 1.9200607e+00 Ad: 1.00e+00 Dobj:
4.0502644e+00
Iter: 15 Ap: 8.44e-01 Pobj: 2.3090128e+00 Ad: 1.00e+00 Dobj:
3.2977094e+00
Iter: 16 Ap: 7.12e-01 Pobj: 2.5196314e+00 Ad: 8.44e-01 Dobj:
2.9513785e+00
Iter: 17 Ap: 1.00e+00 Pobj: 2.6624776e+00 Ad: 1.00e+00 Dobj:
2.8157451e+00
Iter: 18 Ap: 9.72e-01 Pobj: 2.7216848e+00 Ad: 8.32e-01 Dobj:
2.7723998e+00
Iter: 19 Ap: 1.00e+00 Pobj: 2.7394969e+00 Ad: 9.59e-01 Dobj:
2.7553363e+00
Iter: 20 Ap: 9.78e-01 Pobj: 2.7467939e+00 Ad: 1.00e+00 Dobj:
2.7512488e+00
Iter: 21 Ap: 9.06e-01 Pobj: 2.7490423e+00 Ad: 1.00e+00 Dobj:
2.7502901e+00
Iter: 22 Ap: 1.00e+00 Pobj: 2.7493223e+00 Ad: 1.00e+00 Dobj:
2.7502696e+00
Iter: 23 Ap: 1.00e+00 Pobj: 2.7494842e+00 Ad: 1.00e+00 Dobj:
2.7502046e+00
Iter: 24 Ap: 1.00e+00 Pobj: 2.7495981e+00 Ad: 1.00e+00 Dobj:
2.7501589e+00
Iter: 25 Ap: 1.00e+00 Pobj: 2.7498852e+00 Ad: 8.39e-01 Dobj:
2.7500624e+00
Iter: 26 Ap: 1.00e+00 Pobj: 2.7499676e+00 Ad: 9.34e-01 Dobj:
2.7500143e+00
Iter: 27 Ap: 9.67e-01 Pobj: 2.7499915e+00 Ad: 1.00e+00 Dobj:
2.7500019e+00
Iter: 28 Ap: 9.26e-01 Pobj: 2.7499975e+00 Ad: 1.00e+00 Dobj:
2.7500007e+00
Iter: 29 Ap: 9.49e-01 Pobj: 2.7499992e+00 Ad: 1.00e+00 Dobj:
2.7500002e+00
Iter: 30 Ap: 9.34e-01 Pobj: 2.7499998e+00 Ad: 1.00e+00 Dobj:
2.7500001e+00
Iter: 31 Ap: 9.74e-01 Pobj: 2.7499999e+00 Ad: 1.00e+00 Dobj:
2.7500000e+00
Iter: 32 Ap: 6.99e-01 Pobj: 2.7500000e+00 Ad: 8.53e-01 Dobj:
2.7500000e+00
Success: SDP solved
Primal objective value: 2.7500000e+00
Dual objective value: 2.7500000e+00
Relative primal infeasibility: 9.70e-17
Relative dual infeasibility: 7.20e-11
Real Relative Gap: 3.79e-09
XZ Relative Gap: 3.82e-09
DIMACS error measures: 1.05e-16 0.00e+00 2.12e-10 0.00e+00 3.79e-09 3.82e-09
*** glibc detected *** ./example: free(): invalid next size (fast):
0x000000000143dd90 ***
======= Backtrace: =========
/lib/libc.so.6(+0x775b6)[0x7f3c165b45b6]
/lib/libc.so.6(cfree+0x73)[0x7f3c165bae53]
./example[0x402c75]
./example[0x401436]
/lib/libc.so.6(__libc_start_main+0xfd)[0x7f3c1655bc4d]
./example[0x400d89]
======= Memory map: ========
00400000-00421000 r-xp 00000000 07:00 828533
/opt/Csdp-6.1.1/example/example
00620000-00621000 r--p 00020000 07:00 828533
/opt/Csdp-6.1.1/example/example
00621000-00622000 rw-p 00021000 07:00 828533
/opt/Csdp-6.1.1/example/example
0143d000-0145e000 rw-p 00000000 00:00 0
[heap]
7f3c10000000-7f3c10021000 rw-p 00000000 00:00 0
7f3c10021000-7f3c14000000 ---p 00000000 00:00 0
7f3c16326000-7f3c1633c000 r-xp 00000000 07:00 261712
/lib/libgcc_s.so.1
7f3c1633c000-7f3c1653b000 ---p 00016000 07:00 261712
/lib/libgcc_s.so.1
7f3c1653b000-7f3c1653c000 r--p 00015000 07:00 261712
/lib/libgcc_s.so.1
7f3c1653c000-7f3c1653d000 rw-p 00016000 07:00 261712
/lib/libgcc_s.so.1
7f3c1653d000-7f3c166b7000 r-xp 00000000 07:00 261678
/lib/libc-2.11.1.so
7f3c166b7000-7f3c168b6000 ---p 0017a000 07:00 261678
/lib/libc-2.11.1.so
7f3c168b6000-7f3c168ba000 r--p 00179000 07:00 261678
/lib/libc-2.11.1.so
7f3c168ba000-7f3c168bb000 rw-p 0017d000 07:00 261678
/lib/libc-2.11.1.so
7f3c168bb000-7f3c168c0000 rw-p 00000000 00:00 0
7f3c168c0000-7f3c16942000 r-xp 00000000 07:00 261727
/lib/libm-2.11.1.so
7f3c16942000-7f3c16b41000 ---p 00082000 07:00 261727
/lib/libm-2.11.1.so
7f3c16b41000-7f3c16b42000 r--p 00081000 07:00 261727
/lib/libm-2.11.1.so
7f3c16b42000-7f3c16b43000 rw-p 00082000 07:00 261727
/lib/libm-2.11.1.so
7f3c16b43000-7f3c16b79000 r-xp 00000000 07:00 825452
/usr/local/lib/libgslcblas.so.0.0.0
7f3c16b79000-7f3c16d78000 ---p 00036000 07:00 825452
/usr/local/lib/libgslcblas.so.0.0.0
7f3c16d78000-7f3c16d79000 r--p 00035000 07:00 825452
/usr/local/lib/libgslcblas.so.0.0.0
7f3c16d79000-7f3c16d7a000 rw-p 00036000 07:00 825452
/usr/local/lib/libgslcblas.so.0.0.0
7f3c16d7a000-7f3c16d9a000 r-xp 00000000 07:00 261654
/lib/ld-2.11.1.so
7f3c16f79000-7f3c16f7c000 rw-p 00000000 00:00 0
7f3c16f97000-7f3c16f9a000 rw-p 00000000 00:00 0
7f3c16f9a000-7f3c16f9b000 r--p 00020000 07:00 261654
/lib/ld-2.11.1.so
7f3c16f9b000-7f3c16f9c000 rw-p 00021000 07:00 261654
/lib/ld-2.11.1.so
7f3c16f9c000-7f3c16f9d000 rw-p 00000000 00:00 0
7fff2569d000-7fff256b2000 rw-p 00000000 00:00 0
[stack]
7fff2573c000-7fff2573d000 r-xp 00000000 00:00 0
[vdso]
ffffffffff600000-ffffffffff601000 r-xp 00000000 00:00 0
[vsyscall]
Aborted
Any clues as to what this all means?
Best,
Henri Kjellberg
On Wed, Nov 3, 2010 at 7:07 AM, Brian Borchers <borchers at nmt.edu> wrote:
> On Tue, Nov 2, 2010 at 9:08 PM, Henri Kjellberg
> <henrikjellberg at gmail.com> wrote:
>
> > I am trying to figure out if there is some significant change that I
> > am not seeing or if its purely chaning the names of things. For
> > example:
> >
> > In calc_dobj.c we have
> >
> > int k;
> > double *a;
> > double *y;
> > int incx;
> >
> > ddot(&k,a+1,&incx,y+1,&incx);
> >
> > Is it as straight forward as to change it to
> > cblas_ddot(&k,a+1,&incx,y+1,&incx);
> >
>
> I believe that it's just that simple. The CBLAS was defined by a
> group called the BLAS technical forum. There is online documentation
> of the CBLAS interface that would explain any differences.
>
>
>
> --
> Brian Borchers borchers at nmt.edu
> Department of Mathematics http://www.nmt.edu/~borchers/
> New Mexico Tech Phone: (575) 322-2592
> Socorro, NM 87801 FAX: (575) 835-5366
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://list.coin-or.org/pipermail/csdp/attachments/20101107/39b42dd2/attachment.html
More information about the Csdp
mailing list