<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Dear Dune developer team,<div><br></div><div>I just figured out that the reason can be the change in size of the matrix and the need for allocating memory. Let me ask my question in a different way. I have two rectangular dynamic matrices (Dune::DynamicMatrix) and want to multiply them with each other. How can I do that using the available functions (leftmultiply and leftmultiplyany) in dune?</div><div><br></div><div>Thanks,</div><div>Arya</div><div><br></div><div><br><div><div>On Jan 28, 2011, at 2:19 PM, arya fallahi wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Dear Dune developer team,<div><br></div><div>I am trying to multiply to rectangular matrices using the Dune::DenseMatrix class. However, I saw in the source code that one can only use leftmultiply and rightmultiply functions for square matrices. </div><div><span class="Apple-style-span" style="font-family: monospace, fixed; font-size: 15px; "><pre class="fragment" style="font-family: monospace, fixed; font-size: 14px; border-top-width: 1px; border-right-width: 1px; border-bottom-width: 1px; border-left-width: 1px; border-top-style: solid; border-right-style: solid; border-bottom-style: solid; border-left-style: solid; border-top-color: rgb(204, 204, 204); border-right-color: rgb(204, 204, 204); border-bottom-color: rgb(204, 204, 204); border-left-color: rgb(204, 204, 204); background-color: rgb(245, 245, 245); margin-top: 4px; margin-bottom: 4px; margin-left: 2px; margin-right: 8px; padding-left: 6px; padding-right: 6px; padding-top: 4px; padding-bottom: 4px; ">assert(M.<a class="code" href="http://www.dune-project.org/doc/doxygen/html/classDune_1_1DenseMatrix.html#2cfe6247ab4dc74124654344f7d5464e" title="number of rows" style="color: rgb(0, 0, 255); text-decoration: none; font-weight: normal; ">rows</a>() == M.<a class="code" href="http://www.dune-project.org/doc/doxygen/html/classDune_1_1DenseMatrix.html#fb1f8e85eb500e2e823d88fc6a5377a3" title="number of columns" style="color: rgb(0, 0, 255); text-decoration: none; font-weight: normal; ">cols</a>() && M.<a class="code" href="http://www.dune-project.org/doc/doxygen/html/classDune_1_1DenseMatrix.html#2cfe6247ab4dc74124654344f7d5464e" title="number of rows" style="color: rgb(0, 0, 255); text-decoration: none; font-weight: normal; ">rows</a>() == <a class="code" href="http://www.dune-project.org/doc/doxygen/html/classDune_1_1DenseMatrix.html#2cfe6247ab4dc74124654344f7d5464e" title="number of rows" style="color: rgb(0, 0, 255); text-decoration: none; font-weight: normal; ">rows</a>());</pre></span><div>Is there a particular reason for that? why don't you remove the first condition and let dune fulfill the multiplication also for rectangular matrices??</div></div><div><br></div><div>Thanks,</div><div><br><div>
<span class="Apple-style-span" style="border-collapse: separate; font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; font-size: medium; "><span class="Apple-style-span" style="border-collapse: separate; font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; font-size: medium; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Arya Fallahi, <div>Dr. sc. techn. ETH, </div><div>Computational Accelerator Scientist </div><div>Paul Scherrer Institute (PSI), </div><div>CH-5232 Villigen, Switzerland, </div><div><a href="mailto:arya.fallahi@psi.ch">arya.fallahi@psi.ch</a>, </div><div>+41(0)56 310 50 17<br><br></div></div></span></span>
</div>
<br></div></div>_______________________________________________<br>Dune mailing list<br><a href="mailto:Dune@dune-project.org">Dune@dune-project.org</a><br>http://lists.dune-project.org/mailman/listinfo/dune<br></blockquote></div><br>
<br></div></body></html>