visit_handle
SimGetMesh(int domain, const char *name, void *cbdata)
{
visit_handle h = VISIT_INVALID_HANDLE;
if(strcmp(name, "mesh2d") == 0)
{
if(VisIt_CurvilinearMesh_alloc(&h) != VISIT_ERROR)
{
int nTuples;
visit_handle hxy;
nTuples = mesh_dims[domain][0] * mesh_dims[domain][1];
VisIt_VariableData_alloc(&hxy);
VisIt_VariableData_setDataF(hxy, VISIT_OWNER_SIM, 2, nTuples,
mesh_coords[domain]);
VisIt_CurvilinearMesh_setCoords2(h, mesh_dims[domain], hxy);
}
}
return h;
}
visit_handle
SimGetMesh(int domain, const char *name, void *cbdata)
{
visit_handle h = VISIT_INVALID_HANDLE;
if(strcmp(name, "mesh2d") == 0)
{
if(VisIt_RectilinearMesh_alloc(&h) != VISIT_ERROR)
{
visit_handle hxc, hyc;
VisIt_VariableData_alloc(&hxc);
VisIt_VariableData_alloc(&hyc);
VisIt_VariableData_setDataF(hxc, VISIT_OWNER_SIM, 1, rmesh_dims[0], rmesh_x);
VisIt_VariableData_setDataF(hyc, VISIT_OWNER_SIM, 1, rmesh_dims[1], rmesh_y);
VisIt_RectilinearMesh_setCoordsXY(h, hxc, hyc);
}
}
else if(strcmp(name, "mesh3d") == 0)
{
if(VisIt_CurvilinearMesh_alloc(&h) != VISIT_ERROR)
{
int nn;
visit_handle hxc, hyc, hzc;
nn = cmesh_dims[0] * cmesh_dims[1] * cmesh_dims[2];
VisIt_VariableData_alloc(&hxc);
VisIt_VariableData_alloc(&hyc);
VisIt_VariableData_alloc(&hzc);
VisIt_VariableData_setDataF(hxc, VISIT_OWNER_SIM, 1, nn, (float*)cmesh_x);
VisIt_VariableData_setDataF(hyc, VISIT_OWNER_SIM, 1, nn, (float*)cmesh_y);
VisIt_VariableData_setDataF(hzc, VISIT_OWNER_SIM, 1, nn, (float*)cmesh_z);
VisIt_CurvilinearMesh_setCoordsXYZ(h, cmesh_dims, hxc, hyc, hzc);
}
}
return h;
}
visit_handle
SimGetMesh(int domain, const char *name, void *cbdata)
{
visit_handle h = VISIT_INVALID_HANDLE;
simulation_data *sim = (simulation_data *)cbdata;
if(VisIt_CurvilinearMesh_alloc(&h) != VISIT_ERROR)
{
visit_handle hx, hy, hz, hxyz;
if(strcmp(name, "soa_s") == 0)
{
/* Set coordinates using a single variable data made from SOA data. */
serror(VisIt_VariableData_alloc(&hxyz));
serror(VisIt_VariableData_setArrayDataD(hxyz, 0, VISIT_OWNER_SIM, sim->soa.nnodes, 0, sizeof(double), sim->soa.x));
serror(VisIt_VariableData_setArrayDataD(hxyz, 1, VISIT_OWNER_SIM, sim->soa.nnodes, 0, sizeof(double), sim->soa.y));
serror(VisIt_VariableData_setArrayDataD(hxyz, 2, VISIT_OWNER_SIM, sim->soa.nnodes, 0, sizeof(double), sim->soa.z));
serror(VisIt_CurvilinearMesh_setCoords3(h, sim->dims, hxyz));
}
else if(strcmp(name, "soa_m") == 0)
{
/* Set coordinates using 3 variable datas made from SOA data. */
serror(VisIt_VariableData_alloc(&hx));
serror(VisIt_VariableData_alloc(&hy));
serror(VisIt_VariableData_alloc(&hz));
serror(VisIt_VariableData_setDataD(hx, VISIT_OWNER_SIM, 1, sim->soa.nnodes, sim->soa.x));
serror(VisIt_VariableData_setDataD(hy, VISIT_OWNER_SIM, 1, sim->soa.nnodes, sim->soa.y));
serror(VisIt_VariableData_setDataD(hz, VISIT_OWNER_SIM, 1, sim->soa.nnodes, sim->soa.z));
serror(VisIt_CurvilinearMesh_setCoordsXYZ(h, sim->dims, hx, hy, hz));
}
else if(strcmp(name, "soa_i") == 0)
{
/* Set coordinates using a single variable data made from SOA interleaved data. */
serror(VisIt_VariableData_alloc(&hxyz));
serror(VisIt_VariableData_setDataD(hxyz, VISIT_OWNER_SIM, 3, sim->soa.nnodes, sim->soa.xyz));
serror(VisIt_CurvilinearMesh_setCoords3(h, sim->dims, hxyz));
}
else if(strcmp(name, "aos_s") == 0)
{
/* Set coordinates using a single variable data made from AOS data. */
serror(VisIt_VariableData_alloc(&hxyz));
serror(VisIt_VariableData_setArrayDataD(hxyz, 0, VISIT_OWNER_SIM, sim->aos.nnodes, 0, sizeof(Node), &(sim->aos.nodes[0].x)));
serror(VisIt_VariableData_setArrayDataD(hxyz, 1, VISIT_OWNER_SIM, sim->aos.nnodes, 0, sizeof(Node), &(sim->aos.nodes[0].y)));
serror(VisIt_VariableData_setArrayDataD(hxyz, 2, VISIT_OWNER_SIM, sim->aos.nnodes, 0, sizeof(Node), &(sim->aos.nodes[0].z)));
serror(VisIt_CurvilinearMesh_setCoords3(h, sim->dims, hxyz));
}
else if(strcmp(name, "aos_m") == 0)
{
/* Set coordinates using 3 variable datas made from AOS data. */
serror(VisIt_VariableData_alloc(&hx));
serror(VisIt_VariableData_alloc(&hy));
serror(VisIt_VariableData_alloc(&hz));
serror(VisIt_VariableData_setArrayDataD(hx, 0, VISIT_OWNER_SIM, sim->aos.nnodes, 0, sizeof(Node), &(sim->aos.nodes[0].x)));
serror(VisIt_VariableData_setArrayDataD(hy, 0, VISIT_OWNER_SIM, sim->aos.nnodes, 0, sizeof(Node), &(sim->aos.nodes[0].y)));
serror(VisIt_VariableData_setArrayDataD(hz, 0, VISIT_OWNER_SIM, sim->aos.nnodes, 0, sizeof(Node), &(sim->aos.nodes[0].z)));
serror(VisIt_CurvilinearMesh_setCoordsXYZ(h, sim->dims, hx, hy, hz));
}
else if(strcmp(name, "aos_c") == 0)
{
/* Set coordinates using 3 variable datas made from AOS data. Let's COPY data */
serror(VisIt_VariableData_alloc(&hx));
serror(VisIt_VariableData_alloc(&hy));
serror(VisIt_VariableData_alloc(&hz));
serror(VisIt_VariableData_setArrayDataD(hx, 0, VISIT_OWNER_COPY, sim->aos.nnodes, 0, sizeof(Node), &(sim->aos.nodes[0].x)));
serror(VisIt_VariableData_setArrayDataD(hy, 0, VISIT_OWNER_COPY, sim->aos.nnodes, 0, sizeof(Node), &(sim->aos.nodes[0].y)));
serror(VisIt_VariableData_setArrayDataD(hz, 0, VISIT_OWNER_COPY, sim->aos.nnodes, 0, sizeof(Node), &(sim->aos.nodes[0].z)));
serror(VisIt_CurvilinearMesh_setCoordsXYZ(h, sim->dims, hx, hy, hz));
}
}
return h;
}
visit_handle
SimGetMesh(int domain, const char *name, void *cbdata)
{
visit_handle h = VISIT_INVALID_HANDLE;
simulation_data *sim = (simulation_data *)cbdata;
if(strcmp(name, "rect_blank") == 0)
{
if(VisIt_RectilinearMesh_alloc(&h) != VISIT_ERROR)
{
int ncells;
visit_handle hx, hy, gc;
VisIt_VariableData_alloc(&hx);
VisIt_VariableData_alloc(&hy);
VisIt_VariableData_setDataF(hx, VISIT_OWNER_SIM, 1, sim->blankRectMesh.dims[0], sim->blankRectMesh.x);
VisIt_VariableData_setDataF(hy, VISIT_OWNER_SIM, 1, sim->blankRectMesh.dims[1], sim->blankRectMesh.y);
VisIt_RectilinearMesh_setCoordsXY(h, hx, hy);
/* Attach ghost cells to blank out certain cells*/
VisIt_VariableData_alloc(&gc);
ncells = (sim->blankRectMesh.dims[0]-1) * (sim->blankRectMesh.dims[1]-1);
VisIt_VariableData_setDataI(gc, VISIT_OWNER_SIM, 1, ncells, sim->blankRectMesh.ghostCells);
VisIt_RectilinearMesh_setGhostCells(h, gc);
}
}
else if(strcmp(name, "curv_blank") == 0)
{
if(VisIt_CurvilinearMesh_alloc(&h) != VISIT_ERROR)
{
int ncells, nnodes;
visit_handle hx, hy, gc;
nnodes = sim->blankCurvMesh.dims[0] * sim->blankCurvMesh.dims[1];
ncells = (sim->blankCurvMesh.dims[0]-1) * (sim->blankCurvMesh.dims[1]-1);
VisIt_VariableData_alloc(&hx);
VisIt_VariableData_alloc(&hy);
VisIt_VariableData_setDataF(hx, VISIT_OWNER_SIM, 1, nnodes, sim->blankCurvMesh.x);
VisIt_VariableData_setDataF(hy, VISIT_OWNER_SIM, 1, nnodes, sim->blankCurvMesh.y);
VisIt_CurvilinearMesh_setCoordsXY(h, sim->blankCurvMesh.dims, hx, hy);
/* Attach ghost cells to blank out certain cells*/
VisIt_VariableData_alloc(&gc);
VisIt_VariableData_setDataI(gc, VISIT_OWNER_SIM, 1, ncells, sim->blankCurvMesh.ghostCells);
VisIt_CurvilinearMesh_setGhostCells(h, gc);
}
}
else if(strcmp(name, "ucd_blank") == 0)
{
if(VisIt_UnstructuredMesh_alloc(&h) != VISIT_ERROR)
{
visit_handle hxyz, hc, gc;
VisIt_VariableData_alloc(&hxyz);
VisIt_VariableData_setDataF(hxyz, VISIT_OWNER_SIM, 3, sim->blankUcdMesh.nnodes, sim->blankUcdMesh.xyz);
VisIt_UnstructuredMesh_setCoords(h, hxyz);
VisIt_VariableData_alloc(&hc);
VisIt_VariableData_setDataI(hc, VISIT_OWNER_SIM, 1, sim->blankUcdMesh.connectivityLen,
sim->blankUcdMesh.connectivity);
VisIt_UnstructuredMesh_setConnectivity(h, sim->blankUcdMesh.ncells, hc);
/* Attach ghost cells to blank out certain cells*/
VisIt_VariableData_alloc(&gc);
VisIt_VariableData_setDataI(gc, VISIT_OWNER_SIM, 1, sim->blankUcdMesh.ncells, sim->blankUcdMesh.ghostCells);
VisIt_UnstructuredMesh_setGhostCells(h, gc);
}
}
else if(strncmp(name, "multi_domain", 12) == 0)
{
if(VisIt_RectilinearMesh_alloc(&h) != VISIT_ERROR)
{
int ncells;
visit_handle hx, hy, gc;
VisIt_VariableData_alloc(&hx);
VisIt_VariableData_alloc(&hy);
VisIt_VariableData_setDataF(hx, VISIT_OWNER_SIM, 1, sim->multidomain[domain].m.dims[0], sim->multidomain[domain].m.x);
VisIt_VariableData_setDataF(hy, VISIT_OWNER_SIM, 1, sim->multidomain[domain].m.dims[1], sim->multidomain[domain].m.y);
VisIt_RectilinearMesh_setCoordsXY(h, hx, hy);
if(strcmp(name, "multi_domain") == 0)
{
/* Do ghost cells using a ghost cells array. */
VisIt_VariableData_alloc(&gc);
ncells = (sim->multidomain[domain].m.dims[0]-1) * (sim->multidomain[domain].m.dims[1]-1);
VisIt_VariableData_setDataI(gc, VISIT_OWNER_SIM, 1, ncells, sim->multidomain[domain].m.ghostCells);
VisIt_RectilinearMesh_setGhostCells(h, gc);
}
else
{
/* Do ghost cells using indices. Note: other mesh types can also
* designate ghost cells using indices.
*/
VisIt_RectilinearMesh_setRealIndices(h, sim->multidomain[domain].minReal,
sim->multidomain[domain].maxReal);
}
}
}
else if(strcmp(name, "multi_types") == 0)
{
if(domain < 4 && VisIt_RectilinearMesh_alloc(&h) != VISIT_ERROR)
{
int ncells;
visit_handle hx, hy, hz, gc;
VisIt_VariableData_alloc(&hx);
VisIt_VariableData_alloc(&hy);
VisIt_VariableData_alloc(&hz);
VisIt_VariableData_setDataF(hx, VISIT_OWNER_SIM, 1, sim->multipletypes[domain].dims[0], sim->multipletypes[domain].x);
VisIt_VariableData_setDataF(hy, VISIT_OWNER_SIM, 1, sim->multipletypes[domain].dims[1], sim->multipletypes[domain].y);
VisIt_VariableData_setDataF(hz, VISIT_OWNER_SIM, 1, sim->multipletypes[domain].dims[2], sim->multipletypes[domain].z);
VisIt_RectilinearMesh_setCoordsXYZ(h, hx, hy, hz);
/* Do ghost cells using a ghost cells array. */
VisIt_VariableData_alloc(&gc);
ncells = (sim->multipletypes[domain].dims[0]-1) *
(sim->multipletypes[domain].dims[1]-1) *
(sim->multipletypes[domain].dims[2]-1);
VisIt_VariableData_setDataI(gc, VISIT_OWNER_SIM, 1, ncells, sim->multipletypes[domain].ghostCells);
VisIt_RectilinearMesh_setGhostCells(h, gc);
}
}
return h;
}
visit_handle
SimGetMesh(int domain, const char *name, void *cbdata)
{
visit_handle h = VISIT_INVALID_HANDLE;
visit_handle c = VISIT_INVALID_HANDLE;
int owner = VISIT_OWNER_SIM;
int dt = VISIT_DATATYPE_FLOAT;
if(strstr(name, "double") != NULL)
dt = VISIT_DATATYPE_DOUBLE;
if(strstr(name, "dynamic") != NULL)
owner = VISIT_OWNER_VISIT;
if(strstr(name, "curv") != NULL)
{
if(VisIt_CurvilinearMesh_alloc(&h) == VISIT_OKAY &&
VisIt_VariableData_alloc(&c) == VISIT_OKAY)
{
if(strstr(name, "curv2d") != NULL)
{
int dims[2];
dims[0] = NX;
dims[1] = NY;
SetData(c, owner, dt, 2, NX*NY, (float*)coords2d);
VisIt_CurvilinearMesh_setCoords2(h, dims, c);
}
else
{
int dims[3];
dims[0] = NX;
dims[1] = NY;
dims[2] = NZ;
SetData(c, owner, dt, 3, NX*NY*NZ, (float*)coords3d);
VisIt_CurvilinearMesh_setCoords3(h, dims, c);
}
}
}
else if(strstr(name, "point") != NULL)
{
if(VisIt_PointMesh_alloc(&h) == VISIT_OKAY &&
VisIt_VariableData_alloc(&c) == VISIT_OKAY)
{
if(strstr(name, "point2d") != NULL)
SetData(c, owner, dt, 2, NX*NY, (float*)coords2d);
else
SetData(c, owner, dt, 3, NX*NY*NZ, (float*)coords3d);
VisIt_PointMesh_setCoords(h, c);
}
}
else if(strstr(name, "ucd") != NULL)
{
visit_handle hc;
if(VisIt_UnstructuredMesh_alloc(&h) == VISIT_OKAY &&
VisIt_VariableData_alloc(&c) == VISIT_OKAY &&
VisIt_VariableData_alloc(&hc) == VISIT_OKAY)
{
int nzones = 0;
if(strstr(name, "ucd2d") != NULL)
{
int i, j, lconnectivity;
int *connectivity = NULL, *conn = NULL;
nzones = (NX-1)*(NY-1);
lconnectivity = 5*nzones;
conn = connectivity = (int*)malloc(lconnectivity * sizeof(int));
for(j = 0; j < NY-1; ++j)
for(i = 0; i < NX-1; ++i)
{
*conn++ = VISIT_CELL_QUAD;
*conn++ = j*NX + i;
*conn++ = j*NX + (i+1);
*conn++ = (j+1)*NX + (i+1);
*conn++ = (j+1)*NX + i;
}
SetData(c, owner, dt, 2, NX*NY, (float*)coords2d);
VisIt_VariableData_setDataI(hc, VISIT_OWNER_VISIT, 1, lconnectivity,
connectivity);
}
else
{
int i, j, k, lconnectivity;
int *connectivity = NULL, *conn = NULL;
nzones = (NX-1)*(NY-1)*(NZ-1);
lconnectivity = 9*nzones;
conn = connectivity = (int*)malloc(lconnectivity * sizeof(int));
for(k = 0; k < NZ-1; ++k)
for(j = 0; j < NY-1; ++j)
for(i = 0; i < NX-1; ++i)
{
*conn++ = VISIT_CELL_HEX;
*conn++ = k*NX*NY + j*NX + i;
*conn++ = k*NX*NY + j*NX + (i+1);
*conn++ = k*NX*NY + (j+1)*NX + (i+1);
*conn++ = k*NX*NY + (j+1)*NX + i;
*conn++ = (k+1)*NX*NY + j*NX + i;
*conn++ = (k+1)*NX*NY + j*NX + (i+1);
*conn++ = (k+1)*NX*NY + (j+1)*NX + (i+1);
*conn++ = (k+1)*NX*NY + (j+1)*NX + i;
}
SetData(c, owner, dt, 3, NX*NY*NZ, (float*)coords3d);
VisIt_VariableData_setDataI(hc, VISIT_OWNER_VISIT, 1, lconnectivity,
connectivity);
}
VisIt_UnstructuredMesh_setCoords(h, c);
VisIt_UnstructuredMesh_setConnectivity(h, nzones, hc);
}
}
return h;
}
