std::shared_ptr<RenderObject> DisCOVERay::addObject(int sender, const std::string &senderPort,
vistle::Object::const_ptr container,
vistle::Object::const_ptr geometry,
vistle::Object::const_ptr normals,
vistle::Object::const_ptr texture) {
std::shared_ptr<RayRenderObject> ro(new RayRenderObject(m_device, sender, senderPort, container, geometry, normals, texture));
const int t = ro->timestep;
if (t == -1) {
static_geometry.push_back(ro);
} else {
if (anim_geometry.size() <= size_t(t))
anim_geometry.resize(t+1);
anim_geometry[t].push_back(ro);
}
auto rod = ro->data.get();
if (rod->scene) {
RTCGeometry geom_0 = rtcNewGeometry (m_device, RTC_GEOMETRY_TYPE_INSTANCE);
rtcSetGeometryInstancedScene(geom_0,rod->scene);
rtcSetGeometryTimeStepCount(geom_0,1);
rod->instID = rtcAttachGeometry(m_scene,geom_0);
rtcReleaseGeometry(geom_0);
if (instances.size() <= rod->instID)
instances.resize(rod->instID+1);
vassert(!instances[rod->instID]);
instances[rod->instID] = rod;
float transform[16];
auto geoTransform = geometry->getTransform();
for (int i=0; i<16; ++i) {
transform[i] = geoTransform(i%4, i/4);
}
auto inv = geoTransform.inverse().transpose();
for (int c=0; c<3; ++c) {
rod->normalTransform[c].x = inv(c,0);
rod->normalTransform[c].y = inv(c,1);
rod->normalTransform[c].z = inv(c,2);
}
rtcSetGeometryTransform(geom_0,0,RTC_FORMAT_FLOAT4X4_COLUMN_MAJOR,transform);
rtcCommitGeometry(geom_0);
if (t == -1 || t == m_timestep) {
rtcEnableGeometry(geom_0);
m_renderManager.setModified();
} else {
rtcDisableGeometry(geom_0);
}
rtcCommitScene(m_scene);
}
m_renderManager.addObject(ro);
return ro;
}
bool IsoSurface::work(vistle::Object::const_ptr grid,
vistle::Vec<vistle::Scalar>::const_ptr dataS,
vistle::DataBase::const_ptr mapdata) {
const int processorType = getIntParameter("processortype");
#ifdef CUTTINGSURFACE
const Scalar isoValue = 0.0;
#else
const Scalar isoValue = getFloatParameter("isovalue");
#endif
Leveller l(isocontrol, grid, isoValue, processorType);
l.setComputeNormals(m_computeNormals->getValue());
#ifndef CUTTINGSURFACE
l.setIsoData(dataS);
#endif
if(mapdata){
l.addMappedData(mapdata);
};
l.process();
#ifndef CUTTINGSURFACE
auto minmax = dataS->getMinMax();
if (minmax.first[0] < m_min)
m_min = minmax.first[0];
if (minmax.second[0] > m_max)
m_max = minmax.second[0];
#endif
Object::ptr result = l.result();
DataBase::ptr mapresult = l.mapresult();
if (result && !result->isEmpty()) {
#ifndef CUTTINGSURFACE
result->copyAttributes(dataS);
#endif
result->updateInternals();
result->copyAttributes(grid, false);
result->setTransform(grid->getTransform());
if (result->getTimestep() < 0) {
result->setTimestep(grid->getTimestep());
result->setNumTimesteps(grid->getNumTimesteps());
}
if (result->getBlock() < 0) {
result->setBlock(grid->getBlock());
result->setNumBlocks(grid->getNumBlocks());
}
if (mapdata && mapresult) {
mapresult->updateInternals();
mapresult->copyAttributes(mapdata);
mapresult->setGrid(result);
addObject(m_dataOut, mapresult);
}
#ifndef CUTTINGSURFACE
else {
addObject(m_dataOut, result);
}
#endif
}
return true;
}
