<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
<br class="">
<div><br class="">
<blockquote type="cite" class="">
<div class="">On 25 Oct 2019, at 09:23, Erik Rodrígo Jiménez Vázquez <<a href="mailto:erjive@ciencias.unam.mx" class="">erjive@ciencias.unam.mx</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div dir="auto" class="">
<div style="font-family:sans-serif;font-size:12.8px" dir="auto" class="">Dear all</div>
<div style="font-family:sans-serif;font-size:12.8px" dir="auto" class=""><br class="">
</div>
<div style="font-family:sans-serif;font-size:12.8px" dir="auto" class="">I want to subtract an array at t-dt and t times only at the base grid (refinement level 0), following <a href="http://lists.einsteintoolkit.org/pipermail/users/2013-May/003004.html" style="text-decoration-line:none;color:rgb(66,133,244)" class="">http://lists.einsteintoolkit.org/pipermail/users/2013-May/003004.html</a>,
I implemented the suggestion in my code as<br class="">
</div>
<div style="font-family:sans-serif;font-size:12.8px" dir="auto" class="">
<pre style="white-space:pre-wrap" class="">BEGIN_GLOBAL_MODE(cctkGH) {
ENTER_LEVEL_MODE(cctkGH, 0) {
for(int k=0;k<cctk_lsh[2];k++) {
for(int j=0;j<cctk_lsh[1];j++) {
for(int i=0;i<cctk_lsh[0];i++) {
int index= CCTK_GFINDEX3D(cctkGH,i,j,k);
lres += abs(phi[index]-phi_p[index]);
}
}
}
} LEAVE_LEVEL_MODE;
} END_GLOBAL_MODE;</pre>
</div>
<div style="font-family:sans-serif;font-size:12.8px" dir="auto" class="">However I got a segmentation fault. I think is because of the array index, if I comment<br class="">
</div>
<div style="font-family:sans-serif;font-size:12.8px" dir="auto" class="">the line<br class="">
</div>
<div style="font-family:sans-serif;font-size:12.8px" dir="auto" class=""><br class="">
</div>
<div style="font-family:sans-serif;font-size:12.8px" dir="auto" class="">lres += abs(phi[index]-phi_p[index]);<br class="">
</div>
<div style="font-family:sans-serif;font-size:12.8px" dir="auto" class=""><br class="">
</div>
<div style="font-family:sans-serif;font-size:12.8px" dir="auto" class=""> the code works fine.</div>
<div style="font-family:sans-serif;font-size:12.8px" dir="auto" class="">My questions are: 1) How can I perform this operation only at the base grid?, 2) How can I get the correct index at the base level? .. any help is very welcome. Thanks.</div>
</div>
</div>
</blockquote>
<div><br class="">
</div>
<div>Hi Erik,</div>
<div><br class="">
</div>
<div>Can you check whether the phi and phi_p pointers are null?</div>
<div><br class="">
</div>
<div>I'm not sure whether you can access grid data in level mode. </div>
<div><br class="">
</div>
<div>I would have done this by scheduling the function in local mode in schedule.ccl and putting</div>
<div><br class="">
</div>
</div>
<blockquote style="margin: 0 0 0 40px; border: none; padding: 0px;" class="">
<div>
<div>if (cctk_levfac[0] == 1) {</div>
</div>
<div>
<div><br class="">
</div>
</div>
<div>
<div><br class="">
</div>
</div>
<div>
<div>}</div>
</div>
</blockquote>
<br class="">
<div>
<div>around the loop over grid points. I wouldn't have used the global and level mode macros.</div>
<div><br class="">
</div>
<div>cctk_levfac is described in the documentation (<a href="https://einsteintoolkit.org/usersguide/UsersGuidech9.html" class="">https://einsteintoolkit.org/usersguide/UsersGuidech9.html</a>, search for cctk_levfac):</div>
<div><br class="">
</div>
</div>
<blockquote style="margin: 0 0 0 40px; border: none; padding: 0px;" class="">
<div>
<div>cctk_levfac</div>
</div>
<div>
<div>An array of cctk_dim integer factors by which the local grid is refined in the corresponding direction with respect to the base grid.</div>
</div>
</blockquote>
<div>
<div><br class="">
</div>
</div>
<div class="">
<div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
<div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
-- <br class="">
Ian<b class=""><span class="Apple-converted-space"> </span></b>Hinder<br class="">
Research Software Engineer<br class="">
University of Manchester, UK</div>
</div>
</div>
<br class="">
</body>
</html>