예제 #1
0
파일: Rake.cpp 프로젝트: yyr/vapor
//////////////////////////////////////////////////////////////////////////
// generate seeds
//////////////////////////////////////////////////////////////////////////
bool SeedGenerator::GetSeeds(int timestep, VaporFlow* vFlow,
							 float *pSeeds, 
							 const bool bRandom, 
							 const unsigned int randomSeed,
							 int stride)
{
	Rake* pRake = 0;

	if (bRandom) pRake = new SolidRake();
	else {
		if(rakeDimension == POINT)
			pRake = new PointRake();
		else if(rakeDimension == LINE)
			pRake = new LineRake();
		else if(rakeDimension == PLANE)
			pRake = new PlaneRake();
		else if(rakeDimension == SOLID)
			pRake = new SolidRake();
	}
	const vector<double>&usrExts = vFlow->getDataMgr()->GetExtents((size_t)timestep);
	if(bRandom){
		
		//Setup for biased distribution:
		//First calc min/max of field in rake.
		float fieldMin = 0., fieldMax=0.;
		vector<string>varnames;
		
		varnames.push_back(varx);
		varnames.push_back(vary);
		varnames.push_back(varz);
		if( distribBias != 0.){
			bool rc = vFlow->getFieldMagBounds(&fieldMin, &fieldMax, varnames, 
				true, numRefinements, timeStep);
			if (!rc) {
				return false;
			}
		}
		//Then set up the FieldData
		FieldData* fData = 0;
		
		fData = vFlow->setupFieldData(varnames, true, numRefinements, timeStep, false);
		if (!fData) {delete pRake; return false;}
		
		bool rc = pRake->GenSeedBiased(usrExts,distribBias,fieldMin,fieldMax, fData,numSeeds,  rakeLocalMin,rakeLocalMax, pSeeds, randomSeed, stride);
		if (!rc) {
			MyBase::SetErrMsg(VAPOR_ERROR_SEEDS,
				"Unable to generate requested number of random seed points.\nEnsure the rake overlaps valid field data\nor choose a smaller bias.");
			delete fData; delete pRake; return false;
		}
		if (fData) {fData->releaseData(vFlow->getDataMgr());
			delete fData;
		}
	}
	else
		pRake->GenSeedRegular(usrExts,numSeeds, rakeLocalMin, rakeLocalMax, pSeeds, stride);
	delete pRake;
	return true;
}
예제 #2
0
파일: Rake.C 프로젝트: ChengLiOSU/OSUFlow
//////////////////////////////////////////////////////////////////////////
// generate seeds from rake
//////////////////////////////////////////////////////////////////////////
void SeedGenerator::GetSeeds(VECTOR3* pSeeds, const bool bRandom)
{
	Rake* pRake;

	if(rakeDimension == LINE)
		pRake = new LineRake();
	else if(rakeDimension == PLANE)
		pRake = new PlaneRake();
	else if(rakeDimension == SOLID)
		pRake = new SolidRake();
	
	if(bRandom)
		pRake->GenSeedRandom(numSeeds, rakeMin, rakeMax, pSeeds);
	else
		pRake->GenSeedRegular(numSeeds, rakeMin, rakeMax, pSeeds);
	delete pRake;
}