<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
<br>
My real concern with randomness is that results are not reproducible.<br></blockquote><div><br></div><div>My feeling is that this will be small but it will finite in exact arithmetic, which is not ideal.</div><div><br></div><div>If the solve is sensitive to this then I think the you are in trouble anyway (if a plane craters you don't really care how deep the hole is).</div><div><br></div><div>If it did indeed lead to an irreplaceable error then you can look at the verbose output and see that the failures had lower estimates.  At least it is debuggable and since it has never actually happened, has it, then I think we could live with the possibility.</div><div><br></div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
If convergence is sensitive to the RHS used for the estimate, the last<br>
thing I want is to be unable to reproduce that failure.  One could<br>
simply hash the index, but then answers would be different if the<br>
variables are permuted.<br>
<br>
Using the RHS also needs one less vector (not a memory concern for GAMG,<br>
but the cost is not trivial for HPGMG, for example).<br></blockquote><div><br></div><div>You can compute it before you create the RHS, no?  Perhaps we add a SetVectos(Vec[],numvec) that you call before setup() to if you really care about memory?</div><div><br></div></div></div></div>