On Mon, Jul 23, 2012 at 9:18 PM, Nicholas Pezolano <span dir="ltr"><<a href="mailto:npezolano@gmail.com" target="_blank">npezolano@gmail.com</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
@Matt, <br><br> thanks but this is just general matrix multiplication , not the dot product multiplication's for a matrix<br></blockquote><div><br></div><div>Please explain to me the difference.</div><div><br></div><div>
   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>@Jed Brown<br><br>I am trying to draw random numbers from a complex multivariate mixture distribution needed for a monte carlo simulation.  X:= mu + Lambda.* W + gamma.*sqrt(W) , where they are all NxN dense matrices of equal size  where N ranges from 100,000 to 2,000,000. I need to do the dot product multiplication and addition for every part of this distribution. Would petsc be the right tool for this or should i just stick to C and mpi. I've already used many petsc functions that have made life much eaiser<br>


</blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener<br>