
<html><head></head><body><div style="color:#000; background-color:#fff; font-family:標楷體, dfkai-sb;font-size:16px"><div id="yui_3_16_0_ym19_1_1492910496303_2678" dir="ltr">Hello everyone,</div><div id="yui_3_16_0_ym19_1_1492910496303_2678" dir="ltr"><br></div><div id="yui_3_16_0_ym19_1_1492910496303_2678" dir="ltr">I want to manipulate a global vector that has values only on the subgrid of a larger grid obtained by DMDACreate3d(). For example, I'd like to extract a 6x6x6 chunk from a 7x7x7 grid. </div><div id="yui_3_16_0_ym19_1_1492910496303_2678" dir="ltr"><br></div><div id="yui_3_16_0_ym19_1_1492910496303_2678" dir="ltr">At some point, I need to put that 6x6x6 chunk from the logical xyz cartesian coordinates into a natural vector "b" to solve Ax=b at every iteration, then map it back to xyz coordinate system to do further calculations with ghost values.</div><div id="yui_3_16_0_ym19_1_1492910496303_2678" dir="ltr"><br></div><div id="yui_3_16_0_ym19_1_1492910496303_2678" dir="ltr">The closest thing I can compare it to is MPI_Type_create_subarray. Do I just use this or is there a better way?</div><div id="yui_3_16_0_ym19_1_1492910496303_2678" dir="ltr"><br></div><div id="yui_3_16_0_ym19_1_1492910496303_2678" dir="ltr"><br></div><div id="yui_3_16_0_ym19_1_1492910496303_2678" dir="ltr">Thanks for your time,</div><div id="yui_3_16_0_ym19_1_1492910496303_2678" dir="ltr">Tina</div></div></body></html>