[1991] | 1 | #ifndef _Mutation_H__
|
---|
| 2 | #define _Mutation_H__
|
---|
| 3 |
|
---|
| 4 | #include <vector>
|
---|
| 5 | using namespace std;
|
---|
| 6 |
|
---|
| 7 | #include "common.h"
|
---|
| 8 | #include "Halton.h"
|
---|
| 9 |
|
---|
| 10 | namespace GtpVisibilityPreprocessor {
|
---|
| 11 |
|
---|
| 12 | class Vector2;
|
---|
| 13 | class Vector3;
|
---|
| 14 | class VssRay;
|
---|
| 15 | class Preprocessor;
|
---|
| 16 | struct SimpleRay;
|
---|
| 17 | class SimpleRayContainer;
|
---|
| 18 |
|
---|
| 19 | struct VssRayContainer;
|
---|
| 20 |
|
---|
| 21 | class MutationBasedDistribution: public SamplingStrategy
|
---|
| 22 | {
|
---|
| 23 | public:
|
---|
| 24 | MutationBasedDistribution(Preprocessor &preprocessor);
|
---|
| 25 | virtual void Update(VssRayContainer &vssRays);
|
---|
| 26 |
|
---|
| 27 | virtual bool RequiresRays() { return true; }
|
---|
| 28 |
|
---|
| 29 | private:
|
---|
| 30 | virtual bool GenerateSample(SimpleRay &ray);
|
---|
| 31 |
|
---|
| 32 | bool
|
---|
| 33 | GenerateMutation(const int index, SimpleRay &ray);
|
---|
| 34 |
|
---|
| 35 | bool
|
---|
| 36 | GenerateMutationCandidate(const int index,
|
---|
| 37 | SimpleRay &ray,
|
---|
| 38 | Intersectable *object,
|
---|
| 39 | const AxisAlignedBox3 &box
|
---|
| 40 | );
|
---|
| 41 |
|
---|
| 42 |
|
---|
| 43 |
|
---|
| 44 | struct RayEntry {
|
---|
| 45 | // halton sequence for generatin gmutations of this ray
|
---|
| 46 | VssRay *mRay;
|
---|
[2001] | 47 | short mMutations;
|
---|
| 48 | short mUnsuccessfulMutations;
|
---|
[1991] | 49 | HaltonSequence mHalton;
|
---|
| 50 | float mImportance;
|
---|
| 51 | float mCdf;
|
---|
| 52 |
|
---|
[2001] | 53 | Vector3 mutatedOrigin;
|
---|
| 54 | Vector3 mutatedTermination;
|
---|
| 55 |
|
---|
[1991] | 56 | float GetSamplingFactor() const { return mMutations/mImportance; }
|
---|
[2011] | 57 |
|
---|
| 58 | friend bool operator<(const RayEntry &a, const RayEntry &b) {
|
---|
| 59 | return a.GetSamplingFactor() > b.GetSamplingFactor();
|
---|
| 60 | }
|
---|
| 61 |
|
---|
[1991] | 62 | RayEntry() {}
|
---|
| 63 | RayEntry(VssRay *r):mRay(r),
|
---|
| 64 | mMutations(0),
|
---|
| 65 | mUnsuccessfulMutations(0),
|
---|
| 66 | mHalton(),
|
---|
[2001] | 67 | mImportance(1.0f)
|
---|
| 68 | {
|
---|
| 69 | ResetReverseMutation();
|
---|
| 70 | }
|
---|
| 71 |
|
---|
| 72 | void ResetReverseMutation() {
|
---|
| 73 | mutatedOrigin = mutatedTermination = Vector3(0,0,0);
|
---|
| 74 | }
|
---|
| 75 | bool HasReverseMutation() const {
|
---|
| 76 | return !(mutatedOrigin == mutatedTermination);
|
---|
| 77 | }
|
---|
| 78 |
|
---|
| 79 | void SetReverseMutation(const Vector3 &a, const Vector3 &b) {
|
---|
| 80 | mutatedOrigin = a;
|
---|
| 81 | mutatedTermination = b;
|
---|
| 82 | }
|
---|
| 83 |
|
---|
[1991] | 84 | };
|
---|
[2001] | 85 |
|
---|
[1991] | 86 |
|
---|
| 87 | Vector3
|
---|
| 88 | ComputeOriginMutation(const VssRay &ray,
|
---|
| 89 | const Vector3 &U,
|
---|
| 90 | const Vector3 &V,
|
---|
| 91 | const Vector2 vr2,
|
---|
| 92 | const float radius
|
---|
| 93 | );
|
---|
| 94 |
|
---|
| 95 | Vector3
|
---|
| 96 | ComputeTerminationMutation(const VssRay &ray,
|
---|
| 97 | const Vector3 &U,
|
---|
| 98 | const Vector3 &V,
|
---|
| 99 | const Vector2 vr2,
|
---|
| 100 | const float radius
|
---|
| 101 | );
|
---|
| 102 |
|
---|
[1997] | 103 |
|
---|
| 104 | Vector3
|
---|
| 105 | ComputeSilhouetteTerminationMutation(const VssRay &ray,
|
---|
| 106 | const Vector3 &origin,
|
---|
| 107 | const AxisAlignedBox3 &box,
|
---|
| 108 | const Vector3 &U,
|
---|
| 109 | const Vector3 &V,
|
---|
| 110 | const float radius
|
---|
| 111 | );
|
---|
| 112 |
|
---|
[2001] | 113 |
|
---|
| 114 | bool
|
---|
| 115 | ComputeReverseMutation(
|
---|
| 116 | const VssRay &oldRay,
|
---|
| 117 | const VssRay &newRay,
|
---|
| 118 | Vector3 &origin,
|
---|
| 119 | Vector3 &termination
|
---|
| 120 | );
|
---|
| 121 |
|
---|
[2011] | 122 | RayEntry &GetEntry(const int index);
|
---|
| 123 |
|
---|
[1991] | 124 |
|
---|
| 125 | vector<RayEntry> mRays;
|
---|
| 126 | int mBufferStart;
|
---|
| 127 | int mLastIndex;
|
---|
[2060] | 128 |
|
---|
| 129 |
|
---|
| 130 | int mMaxRays;
|
---|
| 131 | float mMutationRadiusOrigin;
|
---|
| 132 | float mMutationRadiusTermination;
|
---|
| 133 | bool mUseReverseSamples;
|
---|
| 134 | float mReverseSamplesDistance;
|
---|
| 135 | bool mUseSilhouetteSamples;
|
---|
| 136 | int mSilhouetteSearchSteps;
|
---|
| 137 | float mSilhouetteProb;
|
---|
| 138 | bool mUsePassImportance;
|
---|
| 139 | bool mUseUnsuccCountImportance;
|
---|
| 140 |
|
---|
[1991] | 141 | };
|
---|
| 142 |
|
---|
| 143 |
|
---|
| 144 | }
|
---|
| 145 |
|
---|
| 146 | #endif
|
---|
| 147 |
|
---|