void PostScriptView::Comments (ostream& out) {
PSVersion(out);
Creator(out);
FontNames(out);
Pages(out);
BoundingBox(out);
out << "%%EndComments\n\n";
}
void create(Triangulation& Tp) {
int N=simu.no_of_particles();
std::vector<Point> points;
// points.reserve(N);
if(simu.create_points()) {
if(simu.at_random()) {
points.reserve(N);
CGAL::Random_points_in_square_2<Point,Creator> g(LL/2.0-0.0001);
CGAL::cpp11::copy_n( g, N, std::back_inserter(points));
cout << N << " particles placed at random" << endl;
} else {
// if((plotting)&&(Nin%2==0)&&(spike)) {
// cout << "Please enter an odd number of particles" << endl;
// std::abort();
// }
int Nb=sqrt(N + 1e-12);
N=Nb*Nb;
simu.set_no_of_particles(N);
points.reserve(N);
cout << N << " particles placed on square lattice" << endl;
FT spacing=LL/FT(Nb+0);
FT side=LL-1*spacing;
points_on_square_grid_2(side/2.0, N, std::back_inserter(points),Creator());;
// for(int i = 0 ; i < Nb ; ++i )
if(simu.perturb()) {
CGAL::perturb_points_2(
points.begin(), points.end(),
simu.pert_rel()* spacing );//,Creator());
cout << "each particle perturbed about " << simu.pert_rel()* spacing << endl;
}
}
}
cout << "Inserting" << endl;
Tp.insert(points.begin(), points.end());
return;
}
ComponentStmOut* TranOutMgr::CreatorComponent(int nKey)
{
ComponentStmOut* pTranOut = LookUp(nKey);
if(pTranOut==NULL)
{
pTranOut = Creator();
pTranOut->Key(nKey);
m_outMap.SetAt(nKey,pTranOut);
}
return pTranOut;
}
void FuildDataMgr::Init(CStdioFile &File)
{
CString strTemp;
FuildDatabase *pFuild = NULL;
while(File.ReadString(strTemp))
{
if(!strTemp.IsEmpty())
{
pFuild = Creator();
pFuild->Init(strTemp);
}
}
}
bool Renderer::addInputObject(int sender, const std::string &senderPort, const std::string & portName,
vistle::Object::const_ptr object) {
int creatorId = object->getCreator();
CreatorMap::iterator it = m_creatorMap.find(creatorId);
if (it != m_creatorMap.end()) {
if (it->second.age < object->getExecutionCounter()) {
//std::cerr << "removing all created by " << creatorId << ", age " << object->getExecutionCounter() << ", was " << it->second.age << std::endl;
removeAllCreatedBy(creatorId);
} else if (it->second.age > object->getExecutionCounter()) {
std::cerr << "received outdated object created by " << creatorId << ", age " << object->getExecutionCounter() << ", was " << it->second.age << std::endl;
return false;
}
} else {
std::string name = getModuleName(object->getCreator());
it = m_creatorMap.insert(std::make_pair(creatorId, Creator(object->getCreator(), name))).first;
}
Creator &creator = it->second;
creator.age = object->getExecutionCounter();
std::shared_ptr<RenderObject> ro;
#if 0
std::cout << "++++++Renderer addInputObject " << object->getType()
<< " creator " << object->getCreator()
<< " exec " << object->getExecutionCounter()
<< " block " << object->getBlock()
<< " timestep " << object->getTimestep() << std::endl;
#endif
if (auto tex = vistle::Texture1D::as(object)) {
if (auto grid = vistle::Coords::as(tex->grid())) {
ro = addObject(sender, senderPort, object, grid, grid->normals(), nullptr, tex);
}
} else if (auto data = vistle::DataBase::as(object)) {
if (auto grid = vistle::Coords::as(data->grid())) {
ro = addObject(sender, senderPort, object, grid, grid->normals(), nullptr, nullptr);
}
}
if (!ro) {
ro = addObject(sender, senderPort, object, object, vistle::Object::ptr(), vistle::Object::ptr(), vistle::Object::ptr());
}
if (ro) {
vassert(ro->timestep >= -1);
if (m_objectList.size() <= size_t(ro->timestep+1))
m_objectList.resize(ro->timestep+2);
m_objectList[ro->timestep+1].push_back(ro);
}
return true;
}
BOOL MOCStationMgr::Init(MOCPipe *pPipe,int nSectionNum)
{
ASSERT(m_pPipe != NULL);
Empty();
if(nSectionNum < 1)
{
ASSERT(FALSE);
return FALSE;
}
m_pPipe = pPipe;
m_dStamR = pPipe->CalcR()/nSectionNum;
MOCStation *pStation = NULL;
for(int i=1;i<nSectionNum;i++)
{
pStation = Creator(i);
Init(pStation,i,nSectionNum);
}
m_dHPress = pPipe->GetDeltaG()/nSectionNum;
return TRUE;
}
void create(void) {
int N=simu.no_of_particles();
std::vector<Point> points;
// points.reserve(N);
if(simu.create_points()) {
if(simu.at_random()) {
points.reserve(N);
CGAL::Random_points_in_square_2<Point,Creator> g(LL/2.0-0.0001);
CGAL::cpp11::copy_n( g, N, std::back_inserter(points));
cout << N << " particles placed at random" << endl;
} else {
// if((plotting)&&(Nin%2==0)&&(spike)) {
// cout << "Please enter an odd number of particles" << endl;
// std::abort();
// }
int Nb=sqrt(N + 1e-12);
N=Nb*Nb;
simu.set_no_of_particles(N);
points.reserve(N);
cout << N << " particles placed on square lattice" << endl;
FT spacing=LL/FT(Nb+0);
FT side=LL-1*spacing;
points_on_square_grid_2(side/2.0, N, std::back_inserter(points),Creator());;
// for(int i = 0 ; i < Nb ; ++i )
if(simu.perturb()) {
CGAL::perturb_points_2(
points.begin(), points.end(),
simu.pert_rel()* spacing );//,Creator());
cout << "each particle perturbed about " << simu.pert_rel()* spacing << endl;
}
}
cout << "Inserting" << endl;
Tp.insert(points.begin(), points.end());
points.clear();
// int Nb = sqrt(N + 1e-12);
// int nm = Nb* simu.mesh_factor() + 1 ;
// int Nm = nm * nm;
int Nm=simu.no_of_nodes();
int nm=sqrt(Nm + 1e-12);
Nm= nm * nm;
simu.set_no_of_nodes(Nm);
points.reserve(Nm);
cout << Nm << " mesh on square lattice" << endl;
FT spacing=LL/FT( nm +0);
FT side=LL-1*spacing;
points_on_square_grid_2(side/2.0, Nm , std::back_inserter(points),Creator());;
// // TODO: perfectly regular square grids are not too good, in fact
// CGAL::perturb_points_2(
// points.begin(), points.end(),
// 0.001* spacing );
Tm.insert(points.begin(), points.end());
} else {
int N=simu.no_of_particles();
char part_file[]="particles.dat";
cout << "reading from file : " << part_file << endl;
std::ifstream main_data;
main_data.open(part_file );
for(int i=0;i<N;i++) {
FT x,y;
main_data >> x;
main_data >> y;
// cout << x << " " << y << endl;
Vertex_handle vh=Tp.insert(Point(x,y));
#include"readin.h"
}
cout << "particles' data read" << endl;
main_data.close();
char mesh_file[]="mesh.dat";
cout << "reading from file : " << mesh_file << endl;
main_data.open(mesh_file );
int Nm=simu.no_of_nodes();
for(int i=0;i<Nm;i++) {
FT x,y;
main_data >> x;
main_data >> y;
// cout << x << " " << y << endl;
Vertex_handle vh=Tm.insert(Point(x,y));
#include"readin.h"
}
cout << "mesh data read" << endl;
main_data.close();
}
// straight from the manual.-
Triangulation::Covering_sheets cs = Tp.number_of_sheets();
cout << "Original covering (particles): " << cs[0] << ' ' << cs[1] << endl;
// return ;
Tp.convert_to_1_sheeted_covering();
cs = Tp.number_of_sheets();
cout << "Current covering (particles): " << cs[0] << ' ' << cs[1] << endl;
return ;
if ( Tp.is_triangulation_in_1_sheet() ) // = true
{
bool is_extensible = Tp.is_extensible_triangulation_in_1_sheet_h1()
|| Tp.is_extensible_triangulation_in_1_sheet_h2(); // = false
Tp.convert_to_1_sheeted_covering();
cs = Tp.number_of_sheets();
cout << "Current covering: " << cs[0] << ' ' << cs[1] << endl;
if ( is_extensible ) // = false
cout << "It is safe to change the triangulation here." << endl;
else {
cout << "It is NOT safe to change the triangulation here!" << endl;
abort();
}
// T.convert_to_9_sheeted_covering();
// cs = T.number_of_sheets();
// cout << "Current covering: " << cs[0] << ' ' << cs[1] << endl;
} else {
cout << "Triangulation not on one sheet!" << endl;
abort();
}
// cout << "It is (again) safe to modify the triangulation." << endl;
return ;
}
void EquivalenceScope::CreateEquivalenceSets(ASTContext &C) {
SetCreator Creator(Connections);
while(Creator.CreateSet(C)) ;
}