close Warning: Can't synchronize with repository "(default)" (/var/svn/tolp does not appear to be a Subversion repository.). Look in the Trac log for more information.

Opened 21 years ago

Closed 21 years ago

Last modified 18 years ago

#152 closed defect (fixed)

TOL Estimation seg faults when there is less data than variables.

Reported by: manuelb Owned by: manuelb
Priority: highest Milestone:
Component: Math Version: head
Severity: major Keywords:
Cc:

Description

When I execute the given estimation, TOL gives this message:

gsl: svd.c:227: ERROR: svd of MxN matrix, M<N, is not implemented
Default GSL error handler invoked.

And then Seg faults.

Attachments (6)

estimation.tol (5.8 KB) - added by manuelb 21 years ago.
TOL source file to test this issue
input.bdt (19.0 KB) - added by manuelb 21 years ago.
Input file for estimation test
input.bst (809 bytes) - added by manuelb 21 years ago.
Input file for estimation test
output.bdt (355 bytes) - added by manuelb 21 years ago.
Output file for estimation test
parest.bst (204 bytes) - added by manuelb 21 years ago.
Input file for estimation test
estimation.tar.gz (2.8 KB) - added by manuelb 21 years ago.
gziped file including all sources and data files

Download all attachments as: .zip

Change History (11)

comment:1 Changed 21 years ago by manuelb

Owner: changed from Víctor de Buen Remiro to manuelb

comment:2 Changed 21 years ago by manuelb

Status: newassigned

Changed 21 years ago by manuelb

Attachment: estimation.tol added

TOL source file to test this issue

Changed 21 years ago by manuelb

Attachment: input.bdt added

Input file for estimation test

Changed 21 years ago by manuelb

Attachment: input.bst added

Input file for estimation test

Changed 21 years ago by manuelb

Attachment: output.bdt added

Output file for estimation test

Changed 21 years ago by manuelb

Attachment: parest.bst added

Input file for estimation test

comment:3 Changed 21 years ago by manuelb

Resolution: fixed
Status: assignedclosed

The problem is that there are more variables than data.

Eventhough it is statistically wrong, that is actually possible, but must be avoided
because the estimation computes the singular value decomposition using the GSL
modified Golub-Reinsch algorithm. This algorithm is very optimized for N>M, (that
is, many more data fields than variables), and has not even been implemented for
N<M.

A fix for this is on the CVS: Prints an error and aborts the estimation.

Changed 21 years ago by manuelb

Attachment: estimation.tar.gz added

gziped file including all sources and data files

comment:4 Changed 21 years ago by manuelb

attachments.isobsolete: 01

comment:5 Changed 18 years ago by Víctor de Buen Remiro

bug_file_loc: http://www.tol-project.org
Note: See TracTickets for help on using tickets.