<html>
  <head>
    <meta content="text/html; charset=ISO-8859-1"
      http-equiv="Content-Type">
  </head>
  <body bgcolor="#FFFFFF" text="#000000">
    <div class="moz-cite-prefix">Le 24/09/2012 22:46, Paul Carrico a
      écrit :<br>
    </div>
    <blockquote
      cite="mid:004001cd9a95$bbbb95d0$3332c170$@carrico@free.fr"
      type="cite">
      <meta http-equiv="Content-Type" content="text/html;
        charset=ISO-8859-1">
      <meta name="Generator" content="Microsoft Word 12 (filtered
        medium)">
      <style><!--
/* Font Definitions */
@font-face
        {font-family:Wingdings;
        panose-1:5 0 0 0 0 0 0 0 0 0;}
@font-face
        {font-family:"Cambria Math";
        panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
        {font-family:Calibri;
        panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0cm;
        margin-bottom:.0001pt;
        font-size:12.0pt;
        font-family:"Times New Roman","serif";}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:blue;
        text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
        {mso-style-priority:99;
        color:purple;
        text-decoration:underline;}
p.MsoListParagraph, li.MsoListParagraph, div.MsoListParagraph
        {mso-style-priority:34;
        margin-top:0cm;
        margin-right:0cm;
        margin-bottom:0cm;
        margin-left:36.0pt;
        margin-bottom:.0001pt;
        font-size:12.0pt;
        font-family:"Times New Roman","serif";}
span.EmailStyle17
        {mso-style-type:personal-compose;
        font-family:"Calibri","sans-serif";
        color:windowtext;}
.MsoChpDefault
        {mso-style-type:export-only;}
@page WordSection1
        {size:612.0pt 792.0pt;
        margin:70.85pt 70.85pt 70.85pt 70.85pt;}
div.WordSection1
        {page:WordSection1;}
/* List Definitions */
@list l0
        {mso-list-id:1890871484;
        mso-list-type:hybrid;
        mso-list-template-ids:-719579868 -1328888404 67895299 67895301 67895297 67895299 67895301 67895297 67895299 67895301;}
@list l0:level1
        {mso-level-start-at:0;
        mso-level-number-format:bullet;
        mso-level-text:-;
        mso-level-tab-stop:none;
        mso-level-number-position:left;
        text-indent:-18.0pt;
        font-family:"Calibri","sans-serif";
        mso-fareast-font-family:"Times New Roman";}
ol
        {margin-bottom:0cm;}
ul
        {margin-bottom:0cm;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
      <div class="WordSection1">
        <p class="MsoNormal"><span
            style="font-family:"Calibri","sans-serif"">Dear
            all,<o:p></o:p></span></p>
        <p class="MsoNormal"><span
            style="font-family:"Calibri","sans-serif""><o:p> </o:p></span></p>
        <p class="MsoNormal"><span
            style="font-family:"Calibri","sans-serif""><o:p> </o:p></span></p>
        <p class="MsoNormal"><span
            style="font-family:"Calibri","sans-serif""
            lang="EN-US">May I ask some advice in matrix calculation
            issue, in order to optimize the code ?<o:p></o:p></span></p>
        <p class="MsoNormal"><span
            style="font-family:"Calibri","sans-serif""
            lang="EN-US"><o:p> </o:p></span></p>
        <p class="MsoNormal"><span
            style="font-family:"Calibri","sans-serif""
            lang="EN-US">The matrix hereafter is the first trial in
            treating a (m x 6) matrix where :<o:p></o:p></span></p>
        <p class="MsoListParagraph"
          style="text-indent:-18.0pt;mso-list:l0 level1 lfo1"><!--[if !supportLists]--><span
style="font-family:"Calibri","sans-serif""
            lang="EN-US"><span style="mso-list:Ignore">-<span
                style="font:7.0pt "Times New Roman"">         
              </span></span></span><!--[endif]--><span
            style="font-family:"Calibri","sans-serif""
            lang="EN-US">The first line is the upper part of a
            symmetrical 6x6 tensor (the 3 first components are the
            diagonal terms)<o:p></o:p></span></p>
        <p class="MsoListParagraph"
          style="text-indent:-18.0pt;mso-list:l0 level1 lfo1"><!--[if !supportLists]--><span
style="font-family:"Calibri","sans-serif""
            lang="EN-US"><span style="mso-list:Ignore">-<span
                style="font:7.0pt "Times New Roman"">         
              </span></span></span><!--[endif]--><span
            style="font-family:"Calibri","sans-serif""
            lang="EN-US">m may correspond to millions of lines<o:p></o:p></span></p>
        <p class="MsoListParagraph"
          style="text-indent:-18.0pt;mso-list:l0 level1 lfo1"><!--[if !supportLists]--><span
style="font-family:"Calibri","sans-serif""
            lang="EN-US"><span style="mso-list:Ignore">-<span
                style="font:7.0pt "Times New Roman"">         
              </span></span></span><!--[endif]--><span
            style="font-family:"Calibri","sans-serif""
            lang="EN-US">(each line is independent from the other)<o:p></o:p></span></p>
        <p class="MsoNormal"><span
            style="font-family:"Calibri","sans-serif""
            lang="EN-US"><o:p> </o:p></span></p>
        <p class="MsoNormal"><span
            style="font-family:"Calibri","sans-serif""
            lang="EN-US">I had a look in some function (such as schur,
            spec, bdiag and so on) in order to see if I can use this
            symmetrical specificity … but I don’t see anything </span><span
            style="font-family:Wingdings" lang="EN-US">è</span><span
            style="font-family:"Calibri","sans-serif""
            lang="EN-US"> Am I wrong ?<o:p></o:p></span></p>
        <p class="MsoNormal"><span
            style="font-family:"Calibri","sans-serif""
            lang="EN-US"><o:p> </o:p></span></p>
      </div>
    </blockquote>
    The <span
      style="font-family:"Calibri","sans-serif""
      lang="EN-US">schur, spec, bdiag function internally check if the
      matrix is symmetric and if yes apply specific algorithm. But the
      matrix must be passed in its full form.<br>
      <br>
    </span>
    <blockquote
      cite="mid:004001cd9a95$bbbb95d0$3332c170$@carrico@free.fr"
      type="cite">
      <div class="WordSection1">
        <p class="MsoNormal"><span
            style="font-family:"Calibri","sans-serif""
            lang="EN-US">Several calculations will be done: Eigen values
            calculation for each submatrix, summation, multiplication
            and so on …<o:p></o:p></span></p>
        <p class="MsoNormal"><span
            style="font-family:"Calibri","sans-serif""
            lang="EN-US"><o:p> </o:p></span></p>
        <p class="MsoNormal"><span
            style="font-family:"Calibri","sans-serif""
            lang="EN-US">Vectorization can be useful to improve speedup
            ?</span></p>
      </div>
    </blockquote>
    I did not see any possible vectorization  (there is no way to
    compute several schur decompositions with only one call to a
    primitive function)<br>
    If you have a multicore processor you can try the parallel_run
    function.<br>
    <br>
    Serge Steer<br>
    INRIA<br>
    <blockquote
      cite="mid:004001cd9a95$bbbb95d0$3332c170$@carrico@free.fr"
      type="cite">
      <div class="WordSection1">
        <p class="MsoNormal"><span
            style="font-family:"Calibri","sans-serif""
            lang="EN-US"><o:p></o:p></span></p>
        <p class="MsoNormal"><span
            style="font-family:"Calibri","sans-serif""
            lang="EN-US"><o:p> </o:p></span></p>
        <p class="MsoNormal"><span
            style="font-family:"Calibri","sans-serif""
            lang="EN-US">Any advices ?<o:p></o:p></span></p>
        <p class="MsoNormal"><span
            style="font-family:"Calibri","sans-serif""
            lang="EN-US"><o:p> </o:p></span></p>
        <p class="MsoNormal"><span
            style="font-family:"Calibri","sans-serif""
            lang="EN-US">Thanks in advance<o:p></o:p></span></p>
        <p class="MsoNormal"><span
            style="font-family:"Calibri","sans-serif""
            lang="EN-US"><o:p> </o:p></span></p>
        <p class="MsoNormal"><span
            style="font-family:"Calibri","sans-serif""
            lang="EN-US">Paul<o:p></o:p></span></p>
        <p class="MsoNormal"><span
            style="font-family:"Calibri","sans-serif""
            lang="EN-US"><o:p> </o:p></span></p>
        <p class="MsoNormal"><span
            style="font-family:"Calibri","sans-serif""
            lang="EN-US">PS : M is made from a text file … its structure
            can be change in order to be as efficient as possible …<o:p></o:p></span></p>
        <p class="MsoListParagraph"><span
            style="font-family:"Calibri","sans-serif""
            lang="EN-US"><o:p> </o:p></span></p>
        <p class="MsoNormal"><span
            style="font-family:"Calibri","sans-serif""
            lang="EN-US">Mat = [<o:p></o:p></span></p>
        <p class="MsoNormal"><span
            style="font-family:"Calibri","sans-serif""
            lang="EN-US">    9.315        9.315        28.41  -
            0.0000999  - 3.169      - 3.169     ; <br>
                2.945        2.945        28.45    0.04826    -
            3.175      - 3.175      ;<br>
                8.486        9.341        25.73    0.0101     -
            3.179      - 3.189      ;<br>
                2.043        2.939        25.73    0.038      -
            3.185      - 3.182      ;<br>
                9.341        8.486        25.73    0.0101     -
            3.189      - 3.179      ;<br>
                2.939        2.043        25.73    0.038      -
            3.182      - 3.185      ;<br>
                8.492        8.492        23.      0.02036    -
            3.179      - 3.179      ;<br>
                2.018        2.018        22.97    0.02779    -
            3.172      - 3.172     ;<br>
                9.108        8.179        25.01    0.01326    -
            3.197      - 1.608     ;<br>
                4.695        2.302        25.04    0.0905     -
            3.199      - 1.61       ;<br>
                8.271        8.186        22.27    0.02824    -
            3.212      - 1.601     ; <br>
                3.842        2.302        22.29    0.07548    -
            3.214      - 1.599     ; <br>
                9.134        8.128        24.91    0.01554    -
            3.278      - 1.61      ;       <br>
                4.675        2.145        24.9     0.08819    -
            3.276      - 1.612     ;   <o:p></o:p></span></p>
        <p class="MsoNormal"><span
            style="font-family:"Calibri","sans-serif""
            lang="EN-US">    …<o:p></o:p></span></p>
        <p class="MsoNormal"><span
            style="font-family:"Calibri","sans-serif""
            lang="EN-US">    …   <br>
                8.291        8.129        22.16    0.03055    -
            3.263      - 1.598  ] ;<o:p></o:p></span></p>
        <p class="MsoNormal"><span
            style="font-family:"Calibri","sans-serif""
            lang="EN-US">    sub11        sub22       sub33   
            sub12        sub23       sub13 … where sub is a submatrix ==
            1 line<o:p></o:p></span></p>
        <p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri","sans-serif""
            lang="EN-US"><o:p> </o:p></span></p>
      </div>
      <br>
      <fieldset class="mimeAttachmentHeader"></fieldset>
      <br>
      <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>
    <br>
  </body>
</html>