Problem with pdsvdtst.f

Post here if you want to report a bug to the LAPACK team

Problem with pdsvdtst.f

Postby gutheil » Thu Feb 23, 2012 8:44 am

The testing of SVD gives always
Test I = 2 has failed

Test I = 3 has failed

Test I = 6 has failed
and so on
The error message is
{ 0, 0}: On entry to PDGESVD parameter number 19 had an illegal value
{ 1, 1}: On entry to PDGESVD parameter number 19 had an illegal value
{ 1, 0}: On entry to PDGESVD parameter number 19 had an illegal value
{ 0, 1}: On entry to PDGESVD parameter number 19 had an illegal value
It can be seen that WPDGESVD changes its value from
{ 1 , 1 }: After dummy Call with -1, WPDGESVD= 1896
{ 1 , 0 }: After dummy Call with -1, WPDGESVD= 1952
{ 0 , 1 }: After dummy Call with -1, WPDGESVD= 1928
{ 0 , 0 }: After dummy Call with -1, WPDGESVD= 2009
to
{ 1 , 1 }: After PDSVDCHK WPDGESVD= 92
{ 1 , 0 }: After PDSVDCHK WPDGESVD= 162
{ 0 , 1 }: After PDSVDCHK WPDGESVD= 155
{ 0 , 0 }: After PDSVDCHK WPDGESVD= 274
and
{ 1 , 1 }: Before call, WPDGESVD= 13
{ 1 , 0 }: Before call, WPDGESVD= 20
{ 0 , 1 }: Before call, WPDGESVD= 13
{ 0 , 0 }: Before call, WPDGESVD= 20
compiled with
mpixlf95_r -qfixed -O3 -g -qmaxmem=-1 -qarch=450 -qtune=450
on a BlueGene/P system.

Inge Gutheil
gutheil
 
Posts: 3
Joined: Thu Feb 23, 2012 5:51 am

Re: Problem with pdsvdtst.f

Postby rodney » Thu Feb 23, 2012 8:35 pm

Hello,

Could you post the SLmake.inc file you are using? What LAPACK and BLAS libraries are you linking against?

Rodney
rodney
 
Posts: 49
Joined: Thu Feb 10, 2011 8:20 pm
Location: Colorado College

Re: Problem with pdsvdtst.f

Postby gutheil » Fri Mar 02, 2012 11:27 am

The SLmake.inc is the following:############################################################################
#
# Program: ScaLAPACK
#
# Module: SLmake.inc
#
# Purpose: Top-level Definitions
#
# Creation date: February 15, 2000
#
# Modified: October 13, 2011
#
# Send bug reports, comments or suggestions to scalapack@cs.utk.edu
#
############################################################################

JUGENEHOME = $(HOME)/Jugene
home = $(JUGENEHOME)/Scalapack/scalapack-2.0.1
SCATEST = $(home)/TESTING
BGP_SYS =/bgsys/drivers/ppcfloor/

#
# C preprocessor definitions: set CDEFS to one of the following:
#
# -DNoChange (fortran subprogram names are lower case without any suffix)
# -DUpCase (fortran subprogram names are upper case without any suffix)
# -DAdd_ (fortran subprogram names are lower case with "_" appended)

CDEFS = -DNoChange
#
# The fortran and C compilers, loaders, and their flags
#

FC = mpixlf95_r -qfixed
CC = mpixlc_r
NOOPT = -O0 -g -I$(BGP_SYS)/include -L$(BGP_SYS)/lib -qarch=450 -qtune=450
FCFLAGS = -O3 -g -qmaxmem=-1 -I$(BGP_SYS)/include -L$(BGP_SYS)/lib -qarch=450 -qtune=450
CCFLAGS = -O3 -g -I$(BGP_SYS)/include -L$(BGP_SYS)/lib -qarch=450 -qtune=450
SRCFLAG =
FCLOADER = $(FC)
CCLOADER = $(CC)
#FCLOADFLAGS =
FCLOADFLAGS = -L$(BGP_SYS)/lib -g -Wl,-allow-multiple-definition
#CCLOADFLAGS =
CCLOADFLAGS = -L$(BGP_SYS)/lib -g

#
# The archiver and the flag(s) to use when building archive (library)
# Also the ranlib routine. If your system has no ranlib, set RANLIB = echo
#

ARCH = ar
ARCHFLAGS = cr
RANLIB = ranlib

#
# The name of the ScaLAPACK library to be created
#

SCALAPACKLIB = libscalapack.a

#
# BLAS, LAPACK (and possibly other) libraries needed for linking test programs
#

BLASLIB = -L/opt/ibmmath/essl/4.4/lib -lesslbg
LAPACKLIB = -L$(LAPACK_LIB) -llapack
LIBS = $(LAPACKLIB) $(BLASLIB)



BLAS is from esslbg as can be seen in the -L/opt/ibmmath/essl/4.4/lib -lesslbg,
LAPACK_LIB is a directory, where a lapack-3.3.0 is installed
gutheil
 
Posts: 3
Joined: Thu Feb 23, 2012 5:51 am

Re: Problem with pdsvdtst.f

Postby rodney » Fri Mar 02, 2012 12:29 pm

I would try two things:

1. reduce the optimization levels to -O1 or -O2; it is possible -O3 is breaking something, it has happened before

2. try using LAPACK 3.4.0 instead of 3.3.0

Rodney
rodney
 
Posts: 49
Joined: Thu Feb 10, 2011 8:20 pm
Location: Colorado College


Return to Bug report

Who is online

Users browsing this forum: No registered users and 2 guests