List Archives >
Maple User Group List Archive >
Archive by date >
This Month By Date >
This Month By Topic
[MUG] Re: Pb with GramShmidt
| [MUG] Re: Pb with GramShmidt |
|
Author: Jmcornil
Posted: Fri, 13 Sep 2002 08:40:16 -0400
|
>> From: jmcornil "jmcornil"
Hello,
When working on an exercice for my students in order to reduce a symmetric matrix
with an orthogonal transition matrix, I fell on the following error :
> with(LinearAlgebra):
> A:=Matrix(3,3,[[2,-2,-2],[2,-2],[2]],scan=triangular,shape=symmetric);
> T:=Eigenvectors(A);
> T2:=T[2];
> B:=[Column(T2,1..3)];
> U:=GramSchmidt(B,'normalized');
Error, (in LinearAlgebra:-LA_Main:-Normalize) datatype of in-place Vector
(algebraic) and datatype of result (anything) do not agree
It is the option 'normalized' which causes the trouble.
I think that the problem comes from the "datatype" of the elements of the
resulting matrix T2 as we can see below but it is impossible to modify this
option.
> VectorOptions(B[1]);
shape = [], datatype = algebraic, orientation = column,
storage = rectangular, order = Fortran_order
> MatrixOptions(T2);
shape = [], datatype = algebraic, storage = rectangular,
order = Fortran_order
Is it a bug ?
Thanks in advance
Jack-Michel CORNIL
|
[View Complete Thread]
Previous by date: [MUG] Comparison of mathematical programs for data analysis, Thomas Schramm
Next by date: [MUG] dchange, Robert J Lopez
Previous thread: [MUG] nasty bug in odeplot with macos X, Roberto Sussman
Next thread: [MUG] dchange, Robert J Lopez
|
(NAG)
