<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;">The method of finite diferences frequently give us systens of sparse matrix. Are you using this one? I had some recent experience with it. When the sparse matrix has a 3 diagonal behavior the most eficient thing I tested was thomas method. But We can't use it when the non 0 elements has no pattern :/<br><br> I've used the finite difference method to solve a problem of the propagation of a wave in a wave guide. In this the sparse system to a "L" shaped wave guide gave me a sparse matrix with no pattern to non 0 elements. In the case of waveguide it lead me to a eigenvalue problem. Unfortunately I didn't found a function like eig() from Matlab, is there something similat in Scilab?<br><br><br>--- Em <b>sáb, 2/10/10, Wojciech Artichowicz <i><w_artichowicz@o2.pl></i></b> escreveu:<br><blockquote style="border-left: 2px solid rgb(16, 16, 255);
 margin-left: 5px; padding-left: 5px;"><br>De: Wojciech Artichowicz <w_artichowicz@o2.pl><br>Assunto: [scilab-Users] Sparse preprocessing<br>Para: users@lists.scilab.org<br>Data: Sábado, 2 de Outubro de 2010, 7:30<br><br><div id="yiv940806513">


 
 
<style>
<!--
#yiv940806513  
#yiv940806513 p.yiv940806513MsoNormal, #yiv940806513 li.yiv940806513MsoNormal, #yiv940806513 div.yiv940806513MsoNormal
        {margin:0cm;margin-bottom:.0001pt;font-size:12.0pt;font-family:"Times New Roman";}
#yiv940806513 a:link, #yiv940806513 span.yiv940806513MsoHyperlink
        {color:blue;text-decoration:underline;}
#yiv940806513 a:visited, #yiv940806513 span.yiv940806513MsoHyperlinkFollowed
        {color:purple;text-decoration:underline;}
#yiv940806513 span.yiv940806513Stylwiadomocie-mail17
        {font-family:Arial;color:windowtext;}
 _filtered #yiv940806513 {margin:70.85pt 70.85pt 70.85pt 70.85pt;}
#yiv940806513 div.yiv940806513Section1
        {}
-->
</style>

<div class="yiv940806513Section1">

<p class="yiv940806513MsoNormal"><font face="Arial" size="2"><span style="font-size: 10pt; font-family: Arial;" lang="EN-US">Hello,</span></font></p> 

<p class="yiv940806513MsoNormal"><font face="Arial" size="2"><span style="font-size: 10pt; font-family: Arial;" lang="EN-US">Fluid mechanics computations often lead to sparse
matrix systems.  Scilab has marvellous sparse functions, however</span></font></p> 

<p class="yiv940806513MsoNormal"><font face="Arial" size="2"><span style="font-size: 10pt; font-family: Arial;" lang="EN-US">I cannot find any function for sparse preprocessing. For
example:</span></font></p> 

<p class="yiv940806513MsoNormal"><font face="Arial" size="2"><span style="font-size: 10pt; font-family: Arial;" lang="EN-US">  </span></font></p> 

<p class="yiv940806513MsoNormal"><font face="Arial" size="2"><span style="font-size: 10pt; font-family: Arial;" lang="EN-US">I have built quite big system with sparse matrix with
no explicit pattern of non-zero elements. The best way to solve this system is
to use</span></font></p> 

<p class="yiv940806513MsoNormal"><font face="Arial" size="2"><span style="font-size: 10pt; font-family: Arial;" lang="EN-US">modified SOR method (my own implementation). My
question is:</span></font></p> 

<p class="yiv940806513MsoNormal"><font face="Arial" size="2"><span style="font-size: 10pt; font-family: Arial;" lang="EN-US">Is there any function in Scilab that would remove
zero entries from diagonal (which isn’t in this case very easy), and that
would do some additional stuff like</span></font></p> 

<p class="yiv940806513MsoNormal"><font face="Arial" size="2"><span style="font-size: 10pt; font-family: Arial;" lang="EN-US">placing at diagonal the biggest possible element.</span></font></p> 

<p class="yiv940806513MsoNormal"><font face="Arial" size="2"><span style="font-size: 10pt; font-family: Arial;" lang="EN-US">  </span></font></p> 

<p class="yiv940806513MsoNormal"><font face="Arial" size="2"><span style="font-size: 10pt; font-family: Arial;" lang="EN-US">If no, is there any well known not too complicated
algorithm for such a purpose?</span></font></p> 

<p class="yiv940806513MsoNormal"><font face="Arial" size="2"><span style="font-size: 10pt; font-family: Arial;" lang="EN-US">  </span></font></p> 

<p class="yiv940806513MsoNormal"><font face="Arial" size="2"><span style="font-size: 10pt; font-family: Arial;" lang="EN-US">Thank You in advance,</span></font></p> 

<p class="yiv940806513MsoNormal"><font face="Arial" size="2"><span style="font-size: 10pt; font-family: Arial;" lang="EN-US">Regards,</span></font></p> 

<p class="yiv940806513MsoNormal"><font face="Arial" size="2"><span style="font-size: 10pt; font-family: Arial;" lang="EN-US">W.A.</span></font></p> 

</div>

 

</div></blockquote></td></tr></table><br>