Context Navigation

source: GTP/trunk/Lib/Vis/Preprocessing/src/ReverseGvs.cpp @ 2686

Revision 2686, 21.3 KB checked in by mattausch, 17 years ago (diff)
fixed several problems

Rev	Line
[2546]	1	#include "SamplingStrategy.h"
	2	#include "Ray.h"
	3	#include "Intersectable.h"
	4	#include "Preprocessor.h"
	5	#include "ViewCellsManager.h"
	6	#include "AxisAlignedBox3.h"
	7	#include "RssTree.h"
	8	#include "Vector2.h"
	9	#include "RndGauss.h"
	10	#include "Exporter.h"
	11	#include "Environment.h"
	12	#include "RayCaster.h"
	13	#include "ReverseGvs.h"
	14
	15
	16	#ifdef GTP_INTERNAL
	17	#include "ArchModeler2MLRT.hxx"
	18	#endif
	19
	20	namespace GtpVisibilityPreprocessor {
	21
	22	#define MUTATION_USE_CDF 0
	23	#define EVALUATE_MUTATION_STATS 1
	24	#define SORT_RAY_ENTRIES 1
	25
	26	// use avg ray contribution as importance
	27	// if 0 the importance is evaluated from the succ of mutations
	28	#define USE_AVG_CONTRIBUTION 1
	29
	30	ReverseGvs::RayEntry &
	31	ReverseGvs::GetEntry(const int index)
	32	{
	33	#if SORT_RAY_ENTRIES
	34	return mRays[index];
	35	#else
	36	return mRays[(mBufferStart+index)%mRays.size()];
	37	#endif
	38	}
	39
	40	void
	41	ReverseGvs::Update(VssRayContainer &vssRays)
	42	{
	43	// for (int i=0; i < mRays.size(); i++)
	44	// cout<<mRays[i].mMutations<<" ";
	45	// cout<<endl;
	46	cerr<<"Mutation update..."<<endl;
	47	cerr<<"rays = "<<mRays.size()<<endl;
	48	if ((int)mRays.size()) {
	49	cerr<<"Oversampling factors = "<<
	50	GetEntry(0).mMutations<<" "<<
	51	GetEntry(1).mMutations<<" "<<
	52	GetEntry(2).mMutations<<" "<<
	53	GetEntry(3).mMutations<<" "<<
	54	GetEntry(4).mMutations<<" "<<
	55	GetEntry(5).mMutations<<" ... "<<
	56	GetEntry((int)mRays.size()-6).mMutations<<" "<<
	57	GetEntry((int)mRays.size()-5).mMutations<<" "<<
	58	GetEntry((int)mRays.size()-4).mMutations<<" "<<
	59	GetEntry((int)mRays.size()-3).mMutations<<" "<<
	60	GetEntry((int)mRays.size()-2).mMutations<<" "<<
	61	GetEntry((int)mRays.size()-1).mMutations<<endl;
	62	}
	63	int contributingRays = 0;
	64
	65	int mutationRays = 0;
	66	int dummyNcMutations = 0;
	67	int dummyCMutations = 0;
	68
	69	int reverseCandidates = 0;
	70
	71	#if 0
	72	sort(mRays.begin(), mRays.end());
	73	// reset the start of the buffer
	74	mBufferStart = 0;
	75	#endif
	76
	77	for (int i=0; i < vssRays.size(); i++) {
	78	if (vssRays[i]->mPvsContribution) {
	79	// reset the counter of unsuccsseful mutation for a generating ray (if it exists)
	80	if (vssRays[i]->mDistribution == MUTATION_BASED_DISTRIBUTION &&
	81	vssRays[i]->mGeneratorId != -1
	82	) {
	83	mRays[vssRays[i]->mGeneratorId].mUnsuccessfulMutations = 0;
	84	#if EVALUATE_MUTATION_STATS
	85	mutationRays++;
	86
	87	Intersectable *newObject = vssRays[i]->mTerminationObject;
	88
	89	Intersectable *oldObject = mRays[vssRays[i]->mGeneratorId].mRay->mTerminationObject;
	90
	91	if (oldObject == newObject)
	92	dummyCMutations++;
	93	#endif
	94	}
	95	contributingRays++;
	96	if (mRays.size() < mMaxRays) {
	97	VssRay newRay = new VssRay(vssRays[i]);
	98	// add this ray
	99	newRay->Ref();
	100	mRays.push_back(RayEntry(newRay));
	101	} else {
	102	// unref the old ray
	103	mRays[mBufferStart].mRay = vssRays[i];
	104	mRays[mBufferStart].mMutations = 0;
	105	mRays[mBufferStart].mUnsuccessfulMutations = 0;
	106	mRays[mBufferStart].ResetReverseMutation();
	107	// mRays[mBufferStart] = RayEntry(newRay);
	108	mBufferStart++;
	109	if (mBufferStart >= mMaxRays)
	110	mBufferStart = 0;
	111	}
	112	} else {
	113	if (vssRays[i]->mDistribution == MUTATION_BASED_DISTRIBUTION &&
	114	vssRays[i]->mGeneratorId != -1
	115	) {
	116	// check whether not to store a new backward mutation candidate
	117	VssRay *oldRay = mRays[vssRays[i]->mGeneratorId].mRay;
	118	VssRay *newRay = vssRays[i];
	119
	120
	121	Intersectable *oldObject = oldRay->mTerminationObject;
	122
	123
	124	if (!mRays[newRay->mGeneratorId].HasReverseMutation()) {
	125	if (DotProd(oldRay->GetDir(), newRay->GetDir()) > 0.0f) {
	126	float oldDist = Magnitude(oldRay->mTermination - newRay->mOrigin);
	127	float newDist = Magnitude(newRay->mTermination - newRay->mOrigin);
	128
	129	if (newDist < oldDist - oldObject->GetBox().Radius()*mReverseSamplesDistance) {
	130	Vector3 origin, termination;
	131	if (ComputeReverseMutation(oldRay, newRay, origin, termination)) {
	132	mRays[newRay->mGeneratorId].SetReverseMutation(origin, termination);
	133	}
	134
	135	reverseCandidates++;
	136	//mReverseCandidates
	137	}
	138	}
	139	}
	140	#if EVALUATE_MUTATION_STATS
	141	mutationRays++;
	142
	143	Intersectable *newObject = vssRays[i]->mTerminationObject;
	144
	145
	146	if (oldObject == newObject)
	147	dummyNcMutations++;
	148	#endif
	149	}
	150	}
	151	}
	152
	153	if (mutationRays) {
	154	cout<<"Mutated rays:"<<mutationRays<<endl;
	155	cout<<"Dummy mutations ratio:"<<100.0f*(dummyCMutations + dummyNcMutations)/
	156	(float)mutationRays<<"%"<<endl;
	157	cout<<"Dummy NC mutations ratio:"<<100.0f*dummyNcMutations/(float)mutationRays<<"%"<<endl;
	158	cout<<"Dummy C mutations ratio:"<<100.0f*dummyCMutations/(float)mutationRays<<"%"<<endl;
	159	cout<<"Reverse candidates:"<<reverseCandidates<<endl;
	160	}
	161
	162	float pContributingRays = contributingRays/(float)vssRays.size();
	163
	164	cout<<"Percentage of contributing rays:"<<pContributingRays<<endl;
	165
	166	if (mUseUnsuccCountImportance) {
	167	// use unsucc mutation samples as feedback on importance
	168	for (int i=0; i < mRays.size(); i++) {
	169	const float minImportance = 0.1f;
	170	const int minImportanceSamples = 20;
	171	mRays[i].mImportance = minImportance +
	172	(1-minImportance)exp(-3.0fmRays[i].mUnsuccessfulMutations/minImportanceSamples);
	173	}
	174	} else {
	175	float importance = 1.0f;
	176	if (mUsePassImportance)
[2686]	177	importance = 1.0f/(pContributingRays + 1e-5f);
[2546]	178	// set this values for last contributingRays
	179	int index = mBufferStart - 1;
	180
	181	for (int i=0; i < contributingRays; i++, index--) {
	182	if (index < 0)
	183	index = (int)mRays.size()-1;
	184	mRays[index].mImportance = importance;
	185	}
	186	}
	187
	188	#if SORT_RAY_ENTRIES
	189	long t1 = GetTime();
	190	sort(mRays.begin(), mRays.end());
	191	// reset the start of the buffer
	192	mBufferStart = 0;
	193	mLastIndex = (int)mRays.size();
	194	cout<<"Mutation candidates sorted in "<<TimeDiff(t1, GetTime())<<" ms."<<endl;
	195	#endif
	196
	197	#if MUTATION_USE_CDF
	198	// compute cdf
	199	mRays[0].mCdf = mRays[0].mImportance/(mRays[0].mMutations+1);
	200	for (int i=1; i < mRays.size(); i++)
	201	mRays[i].mCdf = mRays[i-1].mCdf + mRays[i].mImportance/(mRays[i].mMutations+1);
	202
	203	float scale = 1.0f/mRays[i-1].mCdf;
	204	for (i=0; i < mRays.size(); i++) {
	205	mRays[i].mCdf *= scale;
	206	}
	207	#endif
	208
	209	cout<<"Importance = "<<
	210	GetEntry(0).mImportance<<" "<<
	211	GetEntry((int)mRays.size()-1).mImportance<<endl;
	212
	213	cout<<"Sampling factor = "<<
	214	GetEntry(0).GetSamplingFactor()<<" "<<
	215	GetEntry((int)mRays.size()-1).GetSamplingFactor()<<endl;
	216
	217	cerr<<"Mutation update done."<<endl;
	218	}
	219
	220
	221	Vector3
	222	ReverseGvs::ComputeOriginMutation(const VssRay &ray,
	223	const Vector3 &U,
	224	const Vector3 &V,
	225	const Vector2 vr2,
	226	const float radius
	227	)
	228	{
	229	#if 0
	230	Vector3 v;
	231	if (d.DrivingAxis() == 0)
	232	v = Vector3(0, r[0]-0.5f, r[1]-0.5f);
	233	else
	234	if (d.DrivingAxis() == 1)
	235	v = Vector3(r[0]-0.5f, 0, r[1]-0.5f);
	236	else
	237	v = Vector3(r[0]-0.5f, r[1]-0.5f, 0);
	238	return v(2radius);
	239	#endif
	240	#if 0
	241	return (U(r[0] - 0.5f) + V(r[1] - 0.5f))(2radius);
	242	#endif
	243
	244
	245	// Output random variable
	246	Vector2 gaussvec2;
	247
	248	// Here we apply transform to gaussian, so 2D bivariate
	249	// normal distribution
	250	// float sigma = ComputeSigmaFromRadius(radius);
	251	float sigma = radius;
	252	GaussianOn2D(vr2,
	253	sigma, // input
	254	gaussvec2); // output
	255
	256
	257	// Here we tranform the point correctly to 3D space using base
	258	// vectors of the 3D space defined by the direction
	259	Vector3 shift = gaussvec2.xx * U + gaussvec2.yy * V;
	260
	261	// cout<<shift<<endl;
	262	return shift;
	263	}
	264
	265	Vector3
	266	ReverseGvs::ComputeTerminationMutation(const VssRay &ray,
	267	const Vector3 &U,
	268	const Vector3 &V,
	269	const Vector2 vr2,
	270	const float radius
	271	)
	272	{
	273	#if 0
	274	Vector3 v;
	275	// mutate the termination
	276	if (d.DrivingAxis() == 0)
	277	v = Vector3(0, r[2]-0.5f, r[3]-0.5f);
	278	else
	279	if (d.DrivingAxis() == 1)
	280	v = Vector3(r[2]-0.5f, 0, r[3]-0.5f);
	281	else
	282	v = Vector3(r[2]-0.5f, r[3]-0.5f, 0);
	283
	284	// Vector3 nv;
	285
	286	// if (Magnitude(v) > Limits::Small)
	287	// nv = Normalize(v);
	288	// else
	289	// nv = v;
	290
	291	// v = nvsize + vsize;
	292
	293	return v(4.0fradius);
	294	#endif
	295	#if 0
	296	return (U(vr2.xx - 0.5f) + V(vr2.yy - 0.5f))(4.0fradius);
	297	#endif
	298	Vector2 gaussvec2;
	299	#if 1
	300	float sigma = radius;
	301	GaussianOn2D(vr2,
	302	sigma, // input
	303	gaussvec2); // output
	304	Vector3 shift = gaussvec2.xx * U + gaussvec2.yy * V;
	305	// cout<<shift<<endl;
	306	return shift;
	307	#endif
	308	#if 0
	309	// Here we estimate standard deviation (sigma) from radius
	310	float sigma = 1.1f*ComputeSigmaFromRadius(radius);
	311	Vector3 vr3(vr2.xx, vr2.yy, RandomValue(0,1));
	312	PolarGaussianOnDisk(vr3,
	313	sigma,
	314	radius, // input
	315	gaussvec2); // output
	316
	317	// Here we tranform the point correctly to 3D space using base
	318	// vectors of the 3D space defined by the direction
	319	Vector3 shift = gaussvec2.xx * U + gaussvec2.yy * V;
	320
	321	// cout<<shift<<endl;
	322	return shift;
	323	#endif
	324	}
	325
	326	bool
	327	ReverseGvs::ComputeReverseMutation(
	328	const VssRay &oldRay,
	329	const VssRay &newRay,
	330	Vector3 &origin,
	331	Vector3 &termination
	332	)
	333	{
	334	#ifdef GTP_INTERNAL
	335
	336	// first reconstruct the termination point
	337	Vector3 oldDir = Normalize(oldRay.GetDir());
	338	Plane3 oldPlane(oldDir, oldRay.mTermination);
	339
	340	termination = oldPlane.FindIntersection(newRay.mOrigin,
	341	newRay.mTermination);
	342
	343	// now find the new origin of the ray by casting ray backward from the termination and termining
	344	// silhouette point with respect to the occluding object (object containing the newRay termination)
	345
	346	Plane3 newPlane(oldDir, newRay.mTermination);
	347
	348	Vector3 oldPivot = newPlane.FindIntersection(oldRay.mOrigin,
	349	oldRay.mTermination);
	350
	351	Vector3 newPivot = newRay.mTermination;
	352	Vector3 line = 2.0f*(oldPivot - newPivot);
	353
	354	Intersectable *occluder = newRay.mTerminationObject;
	355
	356	AxisAlignedBox3 box = occluder->GetBox();
	357	box.Scale(2.0f);
	358
	359	const int packetSize = 4;
	360	static int hit_triangles[packetSize];
	361	static float dist[packetSize];
	362	static Vector3 dirs[packetSize];
[2582]	363	static Vector3 origs[packetSize];
[2546]	364	static Vector3 shifts[packetSize];
	365	// now find the silhouette along the line
	366	int i;
	367	float left = 0.0f;
	368	float right = 1.0f;
	369
	370	// cast rays to find silhouette ray
	371	for (int j=0; j < mSilhouetteSearchSteps; j++) {
	372	for (i=0; i < packetSize; i++) {
	373	float r = left + (i+1)*(right-left)/(packetSize+1);
	374	shifts[i] = r*line;
	375	dirs[i] = Normalize(newPivot + shifts[i] - termination );
[2582]	376	origs[i] = termination;
	377	// mlrtaStoreRayASEye4(&termination.x, &dirs[i].x, i);
[2546]	378	}
	379
[2582]	380	// mlrtaTraverseGroupASEye4(&box.Min().x, &box.Max().x, hit_triangles, dist);
	381	assert(preprocessor->mRayCaster);
[2584]	382	preprocessor->mRayCaster->CastRaysPacket4(box.Min(), box.Max(), origs,
[2582]	383	dirs, hit_triangles, dist);
[2546]	384
	385	for (i=0; i < packetSize; i++) {
	386	if (hit_triangles[i] == -1) {
	387	// break on first passing ray
	388	break;
	389	}
	390	}
	391	float rr = left + (i+1)*(right-left)/(packetSize+1);
	392	float rl = left + i*(right-left)/(packetSize+1);
	393	left = rl;
	394	right = rr;
	395	}
	396
	397	float t = right;
	398	if (right==1.0f)
	399	return false;
	400
	401	if (i == packetSize)
	402	origin = newPivot + right*line;
	403	else
	404	origin = newPivot + shifts[i];
	405
	406	if (0) {
	407
	408	static VssRayContainer rRays;
	409	static int counter = 0;
	410	char filename[256];
	411
	412	if (counter < 50) {
	413	sprintf(filename, "reverse_rays_%03d.x3d", counter++);
	414
	415	VssRay tRay(origin, termination, NULL, NULL);
	416	rRays.push_back((VssRay *)&oldRay);
	417	rRays.push_back((VssRay *)&newRay);
	418	rRays.push_back(&tRay);
	419
	420	Exporter *exporter = NULL;
	421	exporter = Exporter::GetExporter(filename);
	422
	423	exporter->SetFilled();
	424
	425	Intersectable *occludee =
	426	oldRay.mTerminationObject;
	427
	428	exporter->SetForcedMaterial(RgbColor(0,0,1));
	429	exporter->ExportIntersectable(occluder);
	430	exporter->SetForcedMaterial(RgbColor(0,1,0));
	431	exporter->ExportIntersectable(occludee);
	432	exporter->ResetForcedMaterial();
	433
	434	exporter->SetWireframe();
	435
	436
	437	exporter->ExportRays(rRays, RgbColor(1, 0, 0));
	438	delete exporter;
	439	rRays.clear();
	440	}
	441	}
	442
	443	#else
	444	cerr << "warning: reverse mutation not supported!" << endl;
	445	#endif
	446
	447	return true;
	448
	449	// now the origin and termination is swapped compred to the generator ray
	450	// swap(origin, termination);???
	451	// -> perhaps not neccessary as the reverse mutation wil only be used once!
	452	}
	453
	454	Vector3
	455	ReverseGvs::ComputeSilhouetteTerminationMutation(const VssRay &ray,
	456	const Vector3 &origin,
	457	const AxisAlignedBox3 &box,
	458	const Vector3 &U,
	459	const Vector3 &V,
	460	const float radius
	461	)
	462	{
	463	#ifdef GTP_INTERNAL
	464
	465	const int packetSize = 4;
	466	static int hit_triangles[packetSize];
	467	static float dist[packetSize];
	468	static Vector3 dirs[packetSize];
[2582]	469	static Vector3 origs[packetSize];
[2546]	470	static Vector3 shifts[packetSize];
	471	// mutate the
	472	float alpha = RandomValue(0.0f, Real(2.0*M_PI));
	473	//float alpha = vr2.x2.0fM_PI;
	474
	475	// direction along which we will mutate the ray
	476	Vector3 line = sin(alpha)U + cos(alpha)V;
	477
	478	// cout<<line<<endl;
	479	// create 16 rays along the selected dir
	480	int i;
	481	float left = 0.0f;
	482	float right = radius;
	483	// cast rays to find silhouette ray
	484	for (int j=0; j < mSilhouetteSearchSteps; j++) {
	485	for (i=0; i < packetSize; i++) {
	486	float r = left + (i+1)*(right-left)/(packetSize+1);
	487	shifts[i] = r*line;
	488	dirs[i] = Normalize(ray.mTermination + shifts[i] - origin );
[2582]	489	// mlrtaStoreRayASEye4(&origin.x, &dirs[i].x, i);
	490	origs[i] = origin;
[2546]	491	}
	492
[2582]	493	// mlrtaTraverseGroupASEye4(&box.Min().x, &box.Max().x, hit_triangles, dist);
	494	assert(preprocessor->mRayCaster);
	495	preprocessor->mRayCaster->CastRaysPacket4(box.Min(), box.Max(), origs,
	496	dirs, hit_triangles, dist);
[2546]	497
	498	for (i=0; i < packetSize; i++) {
	499	if (hit_triangles[i] == -1) {
	500	// if (hit_triangles[i] == -1 \|\| !box.IsInside(origin + dist[i]*dirs[i])) {
	501	// break on first passing ray
	502	break;
	503	}
	504	}
	505	float rr = left + (i+1)*(right-left)/(packetSize+1);
	506	float rl = left + i*(right-left)/(packetSize+1);
	507	left = rl;
	508	right = rr;
	509	}
	510
	511	if (i == packetSize) {
	512	// cerr<<"Warning: hit the same box here should never happen!"<<endl;
	513	// shift the ray even a bit more
	514	//cout<<"W"<<i<<endl;
	515	// return (RandomValue(1.0f, 1.5f)radius)line;
	516	return right*line;
	517	}
	518
	519	// cout<<i<<endl;
	520	return shifts[i];
	521
	522	#else
	523	cerr << "warning: shiluette mutation not supported" << endl;
	524	return Vector3(0, 0, 0);
	525	#endif
	526
	527	}
	528
	529
	530	bool
	531	ReverseGvs::GenerateSample(SimpleRay &sray)
	532	{
	533
	534	if (mRays.size() == 0) {
	535	float rr[5];
	536	// use direction based distribution
	537	Vector3 origin, direction;
	538
	539	for (int i=0; i < 5; i++)
	540	rr[i] = RandomValue(0,1);
	541
	542	mPreprocessor.mViewCellsManager->GetViewPoint(origin,
	543	Vector3(rr[0], rr[1], rr[2]));
	544
	545
	546	direction = UniformRandomVector(rr[3], rr[4]);
	547
	548	const float pdf = 1.0f;
	549	sray = SimpleRay(origin, direction, MUTATION_BASED_DISTRIBUTION, pdf);
	550	sray.mGeneratorId = -1;
	551
	552	return true;
	553	}
	554
	555	int index;
	556
	557	#if !MUTATION_USE_CDF
	558	#if SORT_RAY_ENTRIES
	559	index = mLastIndex - 1;
	560	if (index < 0 \|\| index >= mRays.size()-1) {
	561	index = (int)mRays.size() - 1;
	562	} else
	563	if (
	564	mRays[index].GetSamplingFactor() >= mRays[mLastIndex].GetSamplingFactor()) {
	565	// make another round
	566
	567	// cout<<"R2"<<endl;
	568	// cout<<mLastIndex<<endl;
	569	// cout<<index<<endl;
	570	index = (int)mRays.size() - 1;
	571	}
	572	#else
	573	// get tail of the buffer
	574	index = (mLastIndex+1)%mRays.size();
	575	if (mRays[index].GetSamplingFactor() >
	576	mRays[mLastIndex].GetSamplingFactor()) {
	577	// search back for index where this is valid
	578	index = (mLastIndex - 1 + mRays.size())%mRays.size();
	579	for (int i=0; i < mRays.size(); i++) {
	580
	581	// if (mRays[index].mMutations > mRays[mLastIndex].mMutations)
	582	// break;
	583	if (mRays[index].GetSamplingFactor() >
	584	mRays[mLastIndex].GetSamplingFactor() )
	585	break;
	586	index = (index - 1 + mRays.size())%mRays.size();
	587	}
	588	// go one step back
	589	index = (index+1)%mRays.size();
	590	}
	591	#endif
	592	#else
	593	static HaltonSequence iHalton;
	594	iHalton.GetNext(1, rr);
	595	//rr[0] = RandomValue(0,1);
	596	// use binary search to find index with this cdf
	597	int l=0, r=(int)mRays.size()-1;
	598	while(l<r) {
	599	int i = (l+r)/2;
	600	if (rr[0] < mRays[i].mCdf )
	601	r = i;
	602	else
	603	l = i+1;
	604	}
	605	index = l;
	606	// if (rr[0] >= mRays[r].mCdf)
	607	// index = r;
	608	// else
	609	// index = l;
	610
	611
	612	#endif
	613	// cout<<index<<" "<<rr[0]<<" "<<mRays[index].mCdf<<" "<<mRays[(index+1)%mRays.size()].mCdf<<endl;
	614
	615	mLastIndex = index;
	616	// Debug<<index<<" "<<mRays[index].GetSamplingFactor()<<endl;
	617
	618	if (mRays[index].HasReverseMutation()) {
	619	//cout<<"R "<<mRays[index].mutatedOrigin<<" "<<mRays[index].mutatedTermination<<endl;
	620	sray = SimpleRay(mRays[index].mutatedOrigin,
	621	Normalize(mRays[index].mutatedTermination - mRays[index].mutatedOrigin),
	622	MUTATION_BASED_DISTRIBUTION,
	623	1.0f);
	624
	625
	626	sray.mGeneratorId = index;
	627	mRays[index].ResetReverseMutation();
	628	mRays[index].mMutations++;
	629	mRays[index].mUnsuccessfulMutations++;
	630
	631	return true;
	632	}
	633
	634	return GenerateMutation(index, sray);
	635	}
	636
	637
	638
	639
	640
	641	bool
	642	ReverseGvs::GenerateMutationCandidate(const int index,
	643	SimpleRay &sray,
	644	Intersectable *object,
	645	const AxisAlignedBox3 &box
	646	)
	647	{
	648	float rr[4];
	649
	650	VssRay *ray = mRays[index].mRay;
	651
	652	// rr[0] = RandomValue(0.0f,0.99999f);
	653	// rr[1] = RandomValue(0.0f,0.99999f);
	654	// rr[2] = RandomValue(0.0f,0.99999f);
	655	// rr[3] = RandomValue(0.0f,0.99999f);
	656
	657	// mutate the origin
	658	Vector3 d = ray->GetDir();
	659
	660	float objectRadius = box.Radius();
	661	// cout<<objectRadius<<endl;
	662	if (objectRadius < Limits::Small)
	663	return false;
	664
	665	// Compute right handed coordinate system from direction
	666	Vector3 U, V;
	667	Vector3 nd = Normalize(d);
	668	nd.RightHandedBase(U, V);
	669
	670	Vector3 origin = ray->mOrigin;
	671	Vector3 termination = ray->mTermination; //box.Center(); //ray->mTermination; //box.Center();
	672
	673
	674	// use probabilitistic approach to decide for the type of mutation
	675	float a = RandomValue(0.0f,1.0f);
	676	bool bidirectional = true;
	677
	678	if (mUseSilhouetteSamples && a < mSilhouetteProb) {
	679
	680	termination += ComputeSilhouetteTerminationMutation(*ray,
	681	origin,
	682	box,
	683	U, V,
	684	2.0f*objectRadius);
	685	bidirectional = false;
	686	} else {
	687	mRays[index].mHalton.GetNext(4, rr);
	688
	689	// fuzzy random mutation
	690	origin += ComputeOriginMutation(*ray, U, V,
	691	Vector2(rr[0], rr[1]),
	692	mMutationRadiusOrigin*objectRadius);
	693
	694	termination += ComputeTerminationMutation(*ray, U, V,
	695	Vector2(rr[2], rr[3]),
	696	mMutationRadiusTermination*objectRadius);
	697	}
	698
	699	Vector3 direction = termination - origin;
	700
	701	if (Magnitude(direction) < Limits::Small)
	702	return false;
	703
	704	// shift the origin a little bit
	705	origin += direction*0.5f;
	706
	707	direction.Normalize();
	708
	709	// $$ jb the pdf is yet not correct for all sampling methods!
	710	const float pdf = 1.0f;
	711
	712	sray = SimpleRay(origin, direction, MUTATION_BASED_DISTRIBUTION, pdf);
	713	sray.mGeneratorId = index;
	714	sray.SetBidirectional(bidirectional);
	715
	716	return true;
	717	}
	718
	719	bool
	720	ReverseGvs::GenerateMutation(const int index, SimpleRay &sray)
	721	{
	722	VssRay *ray = mRays[index].mRay;
	723
	724	Intersectable *object = ray->mTerminationObject;
	725
	726	AxisAlignedBox3 box = object->GetBox();
	727
	728	if (GenerateMutationCandidate(index, sray, object, box)) {
	729	mRays[index].mMutations++;
	730	mRays[index].mUnsuccessfulMutations++;
	731
	732	return true;
	733	}
	734	return false;
	735	}
	736
	737
	738
	739	ReverseGvs::ReverseGvs(Preprocessor &preprocessor
	740	) :
	741	SamplingStrategy(preprocessor)
	742	{
	743	mType = MUTATION_BASED_DISTRIBUTION;
	744	mBufferStart = 0;
	745	mLastIndex = 0;
	746
	747	Environment::GetSingleton()->GetIntValue("Mutation.bufferSize",
	748	mMaxRays);
	749
	750	mRays.reserve(mMaxRays);
	751
	752	Environment::GetSingleton()->GetFloatValue("Mutation.radiusOrigin",
	753	mMutationRadiusOrigin);
	754
	755	Environment::GetSingleton()->GetFloatValue("Mutation.radiusTermination",
	756	mMutationRadiusTermination);
	757
	758	bool useEnhancedFeatures;
	759
	760	#ifdef GTP_INTERNAL
	761	int rayCaster;
	762	Environment::GetSingleton()->GetIntValue("Preprocessor.rayCastMethod",
	763	rayCaster);
	764	useEnhancedFeatures = (rayCaster ==
	765	RayCaster::INTEL_RAYCASTER);
	766	#else
	767	useEnhancedFeatures = false;
	768	#endif
	769
	770	if (useEnhancedFeatures) {
	771	Environment::GetSingleton()->GetBoolValue("Mutation.useReverseSamples",
	772	mUseReverseSamples);
	773
	774	Environment::GetSingleton()->GetFloatValue("Mutation.reverseSamplesDistance",
	775
	776	mReverseSamplesDistance);
	777
	778	Environment::GetSingleton()->GetFloatValue("Mutation.silhouetteProb",
	779	mSilhouetteProb);
	780
	781	} else {
	782	mUseReverseSamples = false;
	783	mReverseSamplesDistance = 1e20f;
	784	mSilhouetteProb = 0.0f;
	785	}
	786
	787	Environment::GetSingleton()->GetBoolValue("Mutation.useSilhouetteSamples",
	788	mUseSilhouetteSamples);
	789
	790	Environment::GetSingleton()->GetIntValue("Mutation.silhouetteSearchSteps",
	791	mSilhouetteSearchSteps);
	792
	793
	794	Environment::GetSingleton()->GetBoolValue("Mutation.usePassImportance",
	795	mUsePassImportance);
	796
	797
	798	Environment::GetSingleton()->GetBoolValue("Mutation.useUnsuccCountImportance",
	799	mUseUnsuccCountImportance);
	800
	801
	802
	803	}
	804
	805
	806
	807	}

Note: See TracBrowser for help on using the repository browser.

Download in other formats: