<div dir="ltr">On Mon, May 20, 2013 at 11:35 AM, Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov" target="_blank">jedbrown@mcs.anl.gov</a>></span> wrote:<br><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"><div class="im">Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> writes:<br>
<br>
> I think we have to allow an argument to the strided operations, but the arg<br>
> should be a factor of the bs of the Vec I think.<br>
<br>
</div>What if the Vec has bs=1?<br>
<br>
Maybe it should either be a factor or a multiple.  So you can't use<br>
stride=3 on a bs=5 Vec or vice versa, but you can use stride=2 on a<br>
stride=6 Vec and vice versa.<br>
</blockquote></div><br></div><div class="gmail_extra" style>Actually I now think I was wrong. Let me try to go through the use case, and you tell</div><div class="gmail_extra" style>me where I am wrong:</div><div class="gmail_extra" style>
<br></div><div class="gmail_extra" style>  1) We start with a system with block size B</div><div class="gmail_extra" style><br></div><div class="gmail_extra" style>  2) We split it into k parts, so that \sum_k b_k = B, and b_k is the blocksize of each</div>
<div class="gmail_extra" style>      component. Isn't this how FS for structured grids works?</div><div class="gmail_extra" style><br></div><div class="gmail_extra" style>So the condition should be any bs <= Vec bs, and we also check that it divides the size.</div>
<div class="gmail_extra" style><br></div><div class="gmail_extra" style>    Matt</div><div class="gmail_extra"><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
</div></div>