Le forum de XCAS

Xcas: un logiciel libre de calcul formel
Nous sommes actuellement le Lun Mai 21, 2018 10:52 am

Heures au format UTC




Publier un nouveau sujet Répondre au sujet  [ 4 messages ] 
Auteur Message
 Sujet du message: row reduction without division
MessagePublié: Mar Fév 13, 2018 5:12 pm 
Hors-ligne

Inscrit le: Lun Avr 17, 2017 4:32 pm
Messages: 80
Hi,

Is there an option (e.g. in rref or ref) where I can row reduce a matrix without diivsion? I do sometimes have matrices with symbolic entries and division are not allowed in this case. This is especially important for me when computing the kernel of a matrix with some symbolic entries.

Jose


Haut
 Profil  
 
MessagePublié: Mer Fév 14, 2018 6:33 am 
Hors-ligne

Inscrit le: Mar Déc 20, 2005 4:02 pm
Messages: 4319
Yes, for example rref([[a,b,f],[c,d,g]],keep_pivot)


Haut
 Profil  
 
MessagePublié: Mer Fév 14, 2018 1:35 pm 
Hors-ligne

Inscrit le: Lun Avr 17, 2017 4:32 pm
Messages: 80
Thanks Bernard. Yes you are right. I guess though, kernel computation should still be done manually? The reason is because even after gauss jordan without division there are multiplications involved, and there is a problem if the multiplication is done symbolically because there will be specific cases where the multiplication by zero takes place (increasing the kernel dimension).

My idea is just to get the an nxm matrix (n>=m) M, transpose it, concatenate the transpose with id(n), apply rref with keep_pivot to produce a row-echelon form (not reduced) and then check at each first non-zero symbolic entry of the row and see when the values could possibly be zero. Use the zero conditions to substitute for the variable and apply ker on M for each of these conditions. I will soon provide a minimal example to illustrate.


Haut
 Profil  
 
MessagePublié: Mer Fév 14, 2018 3:11 pm 
Hors-ligne

Inscrit le: Mar Déc 20, 2005 4:02 pm
Messages: 4319
You could perhaps run rref(M,keep_pivot) and look at the diagonal coefficient, if none of them vanish then ker(M) is correct, otherwise you can perhaps express one symbol in terms of other and re-run the same method.


Haut
 Profil  
 
Afficher les messages publiés depuis:  Trier par  
Publier un nouveau sujet Répondre au sujet  [ 4 messages ] 

Heures au format UTC


Qui est en ligne ?

Utilisateurs parcourant actuellement ce forum : Google [Bot] et 2 invités


Vous ne pouvez pas publier de nouveaux sujets dans ce forum
Vous ne pouvez pas répondre aux sujets dans ce forum
Vous ne pouvez pas éditer vos messages dans ce forum
Vous ne pouvez pas supprimer vos messages dans ce forum
Vous ne pouvez pas insérer de pièces jointes dans ce forum

Rechercher pour:
Sauter vers:  
cron
Powered by phpBB® Forum Software © phpBB Group
Traduction réalisée par Maël Soucaze © 2009 phpBB.fr