[prev in list] [next in list] [prev in thread] [next in thread]
List: scilab-users
Subject: Re: [Scilab-users] spec can crash scilab
From: Antoine Monmayrant <antoine.monmayrant () laas ! fr>
Date: 2021-05-05 8:13:39
Message-ID: f96b613e-f47c-e006-3eac-74aba8d662c8 () laas ! fr
[Download RAW message or body]
[Attachment #2 (multipart/alternative)]
Hello all,
Well, my workaround works up to a certain point.
When I increase the size of the matrix under study, I face another
blocking issue:
free(): corrupted unsorted chunks
Antoine
On 05/05/2021 09:23, Antoine Monmayrant wrote:
>
> Hello all,
>
> I might have found a workaround, but I still need to ensure that I get
> exactly the same end results in my code.
> It seems that I can replace:
>
> [W2, gammas] = spec(A); //crashes
>
> by
>
> [al,be,W2]=spec(A,eye(A)); //fine
> gammas=diag(al./be);
>
> If anyone has a solution with eigs, I'll take it.
> Could my workaround help investigate where the bug is in the gateway
> between scilab and lapack?
>
> Cheers,
>
> Antoine
>
> On 05/05/2021 08:22, Antoine Monmayrant wrote:
>>
>>
>> On 05/05/2021 00:09, Samuel Gougeon wrote:
>>> Le 04/05/2021 à 16:38, Antoine Monmayrant a écrit :
>>>>
>>>> Yes, it does the same
>>>>
>>>> So this bug is still a problem, at least on Ubuntu.
>>>> Can you guys try to run spec.tst on his machine and comment on
>>>> https://bugzilla.scilab.org/show_bug.cgi?id=15330 ?
>>>>
>>>
>> Hello Samuel,
>>>
>>> Besides, how confusing and messy the notations on the spec ()
>>> documentation page are! ...
>>>
>>> * R defined twice but differently among arguments
>>> * L defined among arguments but not used in syntaxes
>>> * Q and Z used in syntaxes but not described among arguments
>>> * al, be used in descriptions, instead of alpha and beta listed in
>>> arguments
>>> * E used in the description, but not in arguments nor in syntaxes...
>>>
>>> and so on..
>>> In addition, notations differ from a language version to another one...
>>> Some clarification is urgently required.
>> Indeed.
>> It is also a big problem for portability to get such a basic
>> functionality (solving eigenvalue problems in a linear algebra
>> software!) only working on some platforms!
>>>
>>> In next release, eigs() -- mainly based on the former Arnoldi
>>> module, that becomes internal -- is listed in the same help section.
>>> Have you tried it?
>> Well, I might be a bit far from my field of expertise, but it does
>> not seem that eigs can compute all the eigenvalues/eigenvectors, does it?
>> It seems that by default it computes the 6 biggest ones and can go up
>> to rank(A)-2.
>>
>> In my use case, I need all the eigenvalues...
>>
>> Antoine
>>
>>>
>>> Samuel
>>>
>>>
>>> _______________________________________________
>>> users mailing list
>>> users@lists.scilab.org
>>> http://lists.scilab.org/mailman/listinfo/users
>>
>> _______________________________________________
>> users mailing list
>> users@lists.scilab.org
>> http://lists.scilab.org/mailman/listinfo/users
>
> _______________________________________________
> users mailing list
> users@lists.scilab.org
> http://lists.scilab.org/mailman/listinfo/users
[Attachment #5 (text/html)]
<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
</head>
<body>
<p>Hello all,</p>
<p>Well, my workaround works up to a certain point.<br>
When I increase the size of the matrix under study, I face another
blocking issue:</p>
<p><tt>free(): corrupted unsorted chunks</tt><br>
</p>
<p>Antoine<br>
</p>
<div class="moz-cite-prefix">On 05/05/2021 09:23, Antoine Monmayrant
wrote:<br>
</div>
<blockquote type="cite"
cite="mid:40068fc0-3af3-e2ea-1d9b-795b46576d90@laas.fr">
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
<p>Hello all,</p>
<p>I might have found a workaround, but I still need to ensure
that I get exactly the same end results in my code.<br>
It seems that I can replace:<br>
</p>
<p><tt>[W2, gammas] = spec(A); //crashes</tt><br>
</p>
<p>by<br>
</p>
<p><tt>[al,be,W2]=spec(A,eye(A)); //fine</tt><tt><br>
</tt><tt>gammas=diag(al./be);</tt></p>
<p>If anyone has a solution with eigs, I'll take it.<br>
Could my workaround help investigate where the bug is in the
gateway between scilab and lapack?<br>
</p>
<p>Cheers,</p>
<p>Antoine<br>
</p>
<div class="moz-cite-prefix">On 05/05/2021 08:22, Antoine
Monmayrant wrote:<br>
</div>
<blockquote type="cite"
cite="mid:416a2084-0a27-8604-595d-9c06defcb7a3@laas.fr">
<meta http-equiv="Content-Type" content="text/html;
charset=UTF-8">
<p><br>
</p>
<div class="moz-cite-prefix">On 05/05/2021 00:09, Samuel Gougeon
wrote:<br>
</div>
<blockquote type="cite"
cite="mid:5400650e-5666-83a1-1cbc-1c33d6a447b6@free.fr">
<meta http-equiv="Content-Type" content="text/html;
charset=UTF-8">
<div class="moz-cite-prefix">Le 04/05/2021 Ã 16:38, Antoine
Monmayrant a écrit :<br>
</div>
<blockquote type="cite"
cite="mid:aafebdc4-984b-85a2-01d1-bb03da0df598@laas.fr">
<meta http-equiv="Content-Type" content="text/html;
charset=UTF-8">
<p>Yes, it does the same</p>
<p>So this bug is still a problem, at least on Ubuntu.<br>
Can you guys try to run spec.tst on his machine and
comment on <a class="moz-txt-link-freetext"
href="https://bugzilla.scilab.org/show_bug.cgi?id=15330"
moz-do-not-send="true">https://bugzilla.scilab.org/show_bug.cgi?id=15330</a>
?</p>
</blockquote>
<p><br>
</p>
</blockquote>
Hello Samuel,<br>
<blockquote type="cite"
cite="mid:5400650e-5666-83a1-1cbc-1c33d6a447b6@free.fr">
<p> </p>
<span class="VIiyi" style="display: inline; color: rgb(0, 0, 0); font-size: \
18px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; \
font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; \
text-transform: none; white-space: pre-wrap; word-spacing: 0px; \
-webkit-text-stroke-width: 0px; background-color: rgb(245, 245, 245); \
text-decoration-thickness: initial; text-decoration-style: initial; \
text-decoration-color: initial;" lang="en"><span class="JLqJ4b ChMk0b" \
data-language-for-alternatives="en" data-language-to-translate-into="fr" \
data-phrase-index="0" style="-webkit-tap-highlight-color: transparent; cursor: \
pointer;"><span style="-webkit-tap-highlight-color: transparent;">Besides, how \
confusing and messy the notations on the spec () documentation page are! \
...</span></span></span><span style="color: rgb(0, 0, 0); font-size: 18px; \
font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; \
font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; \
text-transform: none; white-space: pre-wrap; word-spacing: 0px; \
-webkit-text-stroke-width: 0px; background-color: rgb(245, 245, 245); \
text-decoration-thickness: initial; text-decoration-style: initial; \
text-decoration-color: initial; display: inline !important; float: none;"> \
</span><br> <span style="color: rgb(0, 0, 0); font-size: 18px; font-style: normal; \
font-variant-ligatures: normal; font-variant-caps: normal; font-weight: 400; \
letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; \
white-space: pre-wrap; word-spacing: 0px; -webkit-text-stroke-width: 0px; \
background-color: rgb(245, 245, 245); text-decoration-thickness: initial; \
text-decoration-style: initial; text-decoration-color: initial; display: inline \
!important; float: none;"></span> <ul>
<li>R defined twice but differently among arguments</li>
<li>L defined among arguments but not used in syntaxes</li>
<li>Q and Z used in syntaxes but not described among
arguments</li>
<li>al, be used in descriptions, instead of alpha and beta
listed in arguments</li>
<li>E used in the description, but not in arguments nor in
syntaxes...</li>
</ul>
and so on.. <br>
In addition, notations differ from a language version to
another one...<br>
Some clarification is urgently required.<br>
</blockquote>
Indeed.<br>
It is also a big problem for portability to get such a basic
functionality (solving eigenvalue problems in a linear algebra
software!) only working on some platforms!<br>
<blockquote type="cite"
cite="mid:5400650e-5666-83a1-1cbc-1c33d6a447b6@free.fr"> <br>
In next release, eigs() -- mainly based on the former Arnoldi
module, that becomes internal -- is listed in the same help
section. Have you tried it?<br>
</blockquote>
Well, I might be a bit far from my field of expertise, but it
does not seem that eigs can compute all the
eigenvalues/eigenvectors, does it?<br>
It seems that by default it computes the 6 biggest ones and can
go up to rank(A)-2.<br>
<p>In my use case, I need all the eigenvalues...</p>
<p>Antoine<br>
</p>
<blockquote type="cite"
cite="mid:5400650e-5666-83a1-1cbc-1c33d6a447b6@free.fr"> <br>
Samuel<br>
<br>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<pre class="moz-quote-pre" \
wrap="">_______________________________________________ users mailing list
<a class="moz-txt-link-abbreviated" href="mailto:users@lists.scilab.org" \
moz-do-not-send="true">users@lists.scilab.org</a> <a class="moz-txt-link-freetext" \
href="http://lists.scilab.org/mailman/listinfo/users" \
moz-do-not-send="true">http://lists.scilab.org/mailman/listinfo/users</a> </pre>
</blockquote>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<pre class="moz-quote-pre" \
wrap="">_______________________________________________ users mailing list
<a class="moz-txt-link-abbreviated" href="mailto:users@lists.scilab.org" \
moz-do-not-send="true">users@lists.scilab.org</a> <a class="moz-txt-link-freetext" \
href="http://lists.scilab.org/mailman/listinfo/users" \
moz-do-not-send="true">http://lists.scilab.org/mailman/listinfo/users</a> </pre>
</blockquote>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<pre class="moz-quote-pre" \
wrap="">_______________________________________________ users mailing list
<a class="moz-txt-link-abbreviated" \
href="mailto:users@lists.scilab.org">users@lists.scilab.org</a> <a \
class="moz-txt-link-freetext" \
href="http://lists.scilab.org/mailman/listinfo/users">http://lists.scilab.org/mailman/listinfo/users</a>
</pre>
</blockquote>
</body>
</html>
_______________________________________________
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users
[prev in list] [next in list] [prev in thread] [next in thread]
Configure |
About |
News |
Add a list |
Sponsored by KoreLogic