36 | 36 | It turns out that computing this simple global diagnostic is very tricky in parallel applications. It required some thinking to do this in shared-memory and serial with partitions. The issue here is that we are computing something equivalent to the intrinsic F90 function '''MAXLOC'''. As such, we need to get the '''(i,j,k)''' indices of the first occurrence of the maximum value of '''C''' and its qualifiers '''Cu''', '''Cv''', '''Cw''', and '''(i,j,k)'''. It took some serious debugging to get identical values for this new information in serial, serial with partitions, shared-memory, and distributed-memory. Many thanks to Sasha Shchepetkin for his suggestions to compute this. |