Context Navigation

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

Revision 2584, 21.3 KB checked in by mattausch, 16 years ago (diff)

Line
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)
177	importance = 1.0f/(pContributingRays + 1e-5);
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];
363	static Vector3 origs[packetSize];
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 );
376	origs[i] = termination;
377	// mlrtaStoreRayASEye4(&termination.x, &dirs[i].x, i);
378	}
379
380	// mlrtaTraverseGroupASEye4(&box.Min().x, &box.Max().x, hit_triangles, dist);
381	assert(preprocessor->mRayCaster);
382	preprocessor->mRayCaster->CastRaysPacket4(box.Min(), box.Max(), origs,
383	dirs, hit_triangles, dist);
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];
469	static Vector3 origs[packetSize];
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 );
489	// mlrtaStoreRayASEye4(&origin.x, &dirs[i].x, i);
490	origs[i] = origin;
491	}
492
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);
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: