<div dir="ltr">Some points well made here. I have seen in the past job scripts passed on from graduate student to graduate student - the case I am thinking on was an Abaqus script for 8 core systems, being run on a new 32 core system. Why WOULD a graduate student question a script given to them - which works. They should be getting on with their science. I guess this is where Research Software Engineers come in.<div><br></div><div>Another point I would make is about modern processor architectures, for instance AMD Rome/Milan. You can have different Numa Per Socket options, which affect performance. We set the preferred IO path - which I have seen myself to have an effect on latency of MPI messages. IF you are not concerned about your hardware layout you would just go ahead and run, missing  a lot of performance.</div><div><br></div><div>I am now going to be controversial and common that over in Julia land the pattern seems to be these days people develop on their own laptops, or maybe local GPU systems. There is a lot of microbenchmarking going on. But there seems to be not a lot of thought given to CPU pinning or shat happens with hyperthreading. I guess topics like that are part of HPC 'Black Magic' - though I would imagine the low latency crowd are hot on them.</div><div><br></div><div>I often introduce people to the excellent lstopo/hwloc utilities which show the layout of a system. Most people are pleasantly surprised to find this.</div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br><div><br></div><div><br></div><div><br></div><div><br></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, 20 Sept 2021 at 19:28, Lux, Jim (US 7140) via Beowulf <<a href="mailto:beowulf@beowulf.org">beowulf@beowulf.org</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">





<div lang="EN-US" style="overflow-wrap: break-word;">
<div class="gmail-m_-3863071481103716459WordSection1">
<p class="MsoNormal">The recent comments on compilers, caches, etc., are why HPC isn’t a bigger deal.  The infrastructure today is reminiscent of what I used in the 1970s on a big CDC or Burroughs or IBM machine, perhaps with a FPS box attached.
<u></u><u></u></p>
<p class="MsoNormal">I prepare a job, with some sort of job control structure, submit it to a batch queue, and get my results some time later.  Sure, I’m not dropping off a deck or tapes, and I’m not getting green-bar paper or a tape back, but really, it’s
 not much different – I drop a file and get files back either way. <u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">And just like back then, it’s up to me to figure out how best to arrange my code to run fastest (or me, wall clock time, but others it might be CPU time or cost or something else)<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">It would be nice if the compiler (or run-time or infrastructure) figured out the whole “what’s the arrangement of cores/nodes/scratch storage for this application on this particular cluster”. 
<u></u><u></u></p>
<p class="MsoNormal">I also acknowledge that this is a “hard” problem and one that doesn’t have the commercial value of, say, serving the optimum ads to me when I read the newspaper on line.
<u></u><u></u></p>
<p class="MsoNormal"><br>
Yeah, it’s not that hard to call library routines for matrix operations, and to put my trust in the library writers – I trust them more than I trust me to find the fastest linear equation solver, fft, etc. – but so far, the next level of abstraction up – “how
 many cores/nodes” is still left to me, and that means doing instrumentation, figuring out what the results mean, etc.<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<div style="border-right:none;border-bottom:none;border-left:none;border-top:1pt solid rgb(181,196,223);padding:3pt 0in 0in">
<p class="MsoNormal"><b><span style="font-size:12pt;color:black">From: </span></b><span style="font-size:12pt;color:black">Beowulf <<a href="mailto:beowulf-bounces@beowulf.org" target="_blank">beowulf-bounces@beowulf.org</a>> on behalf of "<a href="mailto:beowulf@beowulf.org" target="_blank">beowulf@beowulf.org</a>" <<a href="mailto:beowulf@beowulf.org" target="_blank">beowulf@beowulf.org</a>><br>
<b>Reply-To: </b>Jim Lux <<a href="mailto:james.p.lux@jpl.nasa.gov" target="_blank">james.p.lux@jpl.nasa.gov</a>><br>
<b>Date: </b>Monday, September 20, 2021 at 10:42 AM<br>
<b>To: </b>Lawrence Stewart <<a href="mailto:stewart@serissa.com" target="_blank">stewart@serissa.com</a>>, Jim Cownie <<a href="mailto:jcownie@gmail.com" target="_blank">jcownie@gmail.com</a>><br>
<b>Cc: </b>Douglas Eadline <<a href="mailto:deadline@eadline.org" target="_blank">deadline@eadline.org</a>>, "<a href="mailto:beowulf@beowulf.org" target="_blank">beowulf@beowulf.org</a>" <<a href="mailto:beowulf@beowulf.org" target="_blank">beowulf@beowulf.org</a>><br>
<b>Subject: </b>Re: [Beowulf] [EXTERNAL] Re: Deskside clusters<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<div style="border-right:none;border-bottom:none;border-left:none;border-top:1pt solid rgb(181,196,223);padding:3pt 0in 0in">
<p class="MsoNormal"><b><span style="font-size:12pt;color:black">From: </span></b><span style="font-size:12pt;color:black">Beowulf <<a href="mailto:beowulf-bounces@beowulf.org" target="_blank">beowulf-bounces@beowulf.org</a>> on behalf of Lawrence Stewart <<a href="mailto:stewart@serissa.com" target="_blank">stewart@serissa.com</a>><br>
<b>Date: </b>Monday, September 20, 2021 at 9:17 AM<br>
<b>To: </b>Jim Cownie <<a href="mailto:jcownie@gmail.com" target="_blank">jcownie@gmail.com</a>><br>
<b>Cc: </b>Lawrence Stewart <<a href="mailto:stewart@serissa.com" target="_blank">stewart@serissa.com</a>>, Douglas Eadline <<a href="mailto:deadline@eadline.org" target="_blank">deadline@eadline.org</a>>, "<a href="mailto:beowulf@beowulf.org" target="_blank">beowulf@beowulf.org</a>" <<a href="mailto:beowulf@beowulf.org" target="_blank">beowulf@beowulf.org</a>><br>
<b>Subject: </b>Re: [Beowulf] [EXTERNAL] Re: Deskside clusters<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<p class="MsoNormal">Well said.  Expanding on this, caches work because of both temporal locality and
<u></u><u></u></p>
<div>
<p class="MsoNormal">spatial locality.  Spatial locality is addressed by having cache lines be substantially<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">larger than a byte or word.  These days, 64 bytes is pretty common.  Some prefetch schemes, <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">like the L1D version that fetches the VA ^ 64 clearly affect spatial locality.  Streaming <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">prefetch has an expanded notion of “spatial” I suppose!<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">What puzzles me is why compilers seem not to have evolved much notion of cache management. It <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">seems like something a smart compiler could do.  Instead, it is left to Prof. Goto and the folks<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">at ATLAS and BLIS to figure out how to rewrite algorithms for efficient cache behavior. To my<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">limited knowledge, compilers don’t make much use of PREFETCH or any non-temporal loads and stores<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">either. It seems to me that once the programmer helps with RESTRICT and so forth, then compilers could perfectly well dynamically move parts of arrays around to maximize cache use.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">-L<u></u><u></u></p>
<p class="gmail-m_-3863071481103716459MsoListParagraph" style="margin-left:0in"><u></u> <u></u></p>
<p class="gmail-m_-3863071481103716459MsoListParagraph" style="margin-left:0in">I suspect that there’s enough variability among cache implementation and the wide variety of algorithms that might use it that writing a smart-enough compiler is “hard” and “expensive”.<u></u><u></u></p>
<p class="gmail-m_-3863071481103716459MsoListParagraph" style="margin-left:0in"><u></u> <u></u></p>
<p class="gmail-m_-3863071481103716459MsoListParagraph" style="margin-left:0in">Leaving it to the library authors is probably the best “bang for the buck”. 
<u></u><u></u></p>
<p class="gmail-m_-3863071481103716459MsoListParagraph" style="margin-left:0in"><u></u> <u></u></p>
<p class="gmail-m_-3863071481103716459MsoListParagraph" style="margin-left:0in"><u></u> <u></u></p>
<p class="gmail-m_-3863071481103716459MsoListParagraph" style="margin-left:0in"><u></u> <u></u></p>
</div>
</div>
</div>

_______________________________________________<br>
Beowulf mailing list, <a href="mailto:Beowulf@beowulf.org" target="_blank">Beowulf@beowulf.org</a> sponsored by Penguin Computing<br>
To change your subscription (digest mode or unsubscribe) visit <a href="https://beowulf.org/cgi-bin/mailman/listinfo/beowulf" rel="noreferrer" target="_blank">https://beowulf.org/cgi-bin/mailman/listinfo/beowulf</a><br>
</blockquote></div>