instance_intersector1.cpp 5.9 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149
  1. // ======================================================================== //
  2. // Copyright 2009-2017 Intel Corporation //
  3. // //
  4. // Licensed under the Apache License, Version 2.0 (the "License"); //
  5. // you may not use this file except in compliance with the License. //
  6. // You may obtain a copy of the License at //
  7. // //
  8. // http://www.apache.org/licenses/LICENSE-2.0 //
  9. // //
  10. // Unless required by applicable law or agreed to in writing, software //
  11. // distributed under the License is distributed on an "AS IS" BASIS, //
  12. // WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. //
  13. // See the License for the specific language governing permissions and //
  14. // limitations under the License. //
  15. // ======================================================================== //
  16. #include "instance_intersector1.h"
  17. #include "../common/scene.h"
  18. namespace embree
  19. {
  20. namespace isa
  21. {
  22. void InstanceBoundsFunction(void* userPtr, const Instance* instance, size_t item, size_t itime, BBox3fa& bounds_o)
  23. {
  24. assert(itime < instance->numTimeSteps);
  25. unsigned num_time_segments = instance->numTimeSegments();
  26. if (num_time_segments == 0) {
  27. bounds_o = xfmBounds(instance->local2world[itime],instance->object->bounds.bounds());
  28. }
  29. else {
  30. const float ftime = float(itime) / float(num_time_segments);
  31. const BBox3fa obounds = instance->object->bounds.interpolate(ftime);
  32. bounds_o = xfmBounds(instance->local2world[itime],obounds);
  33. }
  34. }
  35. RTCBoundsFunc3 InstanceBoundsFunc() {
  36. return (RTCBoundsFunc3) InstanceBoundsFunction;
  37. }
  38. void FastInstanceIntersector1::intersect(const Instance* instance, Ray& ray, size_t item)
  39. {
  40. const AffineSpace3fa world2local =
  41. likely(instance->numTimeSteps == 1) ? instance->getWorld2Local() : instance->getWorld2Local(ray.time);
  42. const Vec3fa ray_org = ray.org;
  43. const Vec3fa ray_dir = ray.dir;
  44. const int ray_geomID = ray.geomID;
  45. const int ray_instID = ray.instID;
  46. ray.org = xfmPoint (world2local,ray_org);
  47. ray.dir = xfmVector(world2local,ray_dir);
  48. ray.geomID = RTC_INVALID_GEOMETRY_ID;
  49. ray.instID = instance->id;
  50. IntersectContext context(instance->object,nullptr);
  51. instance->object->intersect((RTCRay&)ray,&context);
  52. ray.org = ray_org;
  53. ray.dir = ray_dir;
  54. if (ray.geomID == RTC_INVALID_GEOMETRY_ID) {
  55. ray.geomID = ray_geomID;
  56. ray.instID = ray_instID;
  57. }
  58. }
  59. void FastInstanceIntersector1::occluded (const Instance* instance, Ray& ray, size_t item)
  60. {
  61. const AffineSpace3fa world2local =
  62. likely(instance->numTimeSteps == 1) ? instance->getWorld2Local() : instance->getWorld2Local(ray.time);
  63. const Vec3fa ray_org = ray.org;
  64. const Vec3fa ray_dir = ray.dir;
  65. ray.org = xfmPoint (world2local,ray_org);
  66. ray.dir = xfmVector(world2local,ray_dir);
  67. ray.instID = instance->id;
  68. IntersectContext context(instance->object,nullptr);
  69. instance->object->occluded((RTCRay&)ray,&context);
  70. ray.org = ray_org;
  71. ray.dir = ray_dir;
  72. }
  73. DEFINE_SET_INTERSECTOR1(InstanceIntersector1,FastInstanceIntersector1);
  74. void FastInstanceIntersector1M::intersect(const Instance* instance, RTCIntersectContext* context, Ray** rays, size_t M, size_t item)
  75. {
  76. assert(M<=MAX_INTERNAL_STREAM_SIZE);
  77. Ray lrays[MAX_INTERNAL_STREAM_SIZE];
  78. AffineSpace3fa world2local = instance->getWorld2Local();
  79. for (size_t i=0; i<M; i++)
  80. {
  81. if (unlikely(instance->numTimeSteps != 1))
  82. world2local = instance->getWorld2Local(rays[i]->time);
  83. lrays[i].org = xfmPoint (world2local,rays[i]->org);
  84. lrays[i].dir = xfmVector(world2local,rays[i]->dir);
  85. lrays[i].tnear = rays[i]->tnear;
  86. lrays[i].tfar = rays[i]->tfar;
  87. lrays[i].time = rays[i]->time;
  88. lrays[i].mask = rays[i]->mask;
  89. lrays[i].geomID = RTC_INVALID_GEOMETRY_ID;
  90. lrays[i].instID = instance->id;
  91. }
  92. rtcIntersect1M((RTCScene)instance->object,context,(RTCRay*)lrays,M,sizeof(Ray));
  93. for (size_t i=0; i<M; i++)
  94. {
  95. if (lrays[i].geomID == RTC_INVALID_GEOMETRY_ID) continue;
  96. rays[i]->instID = lrays[i].instID;
  97. rays[i]->geomID = lrays[i].geomID;
  98. rays[i]->primID = lrays[i].primID;
  99. rays[i]->u = lrays[i].u;
  100. rays[i]->v = lrays[i].v;
  101. rays[i]->tfar = lrays[i].tfar;
  102. rays[i]->Ng = lrays[i].Ng;
  103. }
  104. }
  105. void FastInstanceIntersector1M::occluded (const Instance* instance, RTCIntersectContext* context, Ray** rays, size_t M, size_t item)
  106. {
  107. assert(M<MAX_INTERNAL_STREAM_SIZE);
  108. Ray lrays[MAX_INTERNAL_STREAM_SIZE];
  109. AffineSpace3fa world2local = instance->getWorld2Local();
  110. for (size_t i=0; i<M; i++)
  111. {
  112. if (unlikely(instance->numTimeSteps != 1))
  113. world2local = instance->getWorld2Local(rays[i]->time);
  114. lrays[i].org = xfmPoint (world2local,rays[i]->org);
  115. lrays[i].dir = xfmVector(world2local,rays[i]->dir);
  116. lrays[i].tnear = rays[i]->tnear;
  117. lrays[i].tfar = rays[i]->tfar;
  118. lrays[i].time = rays[i]->time;
  119. lrays[i].mask = rays[i]->mask;
  120. lrays[i].geomID = RTC_INVALID_GEOMETRY_ID;
  121. lrays[i].instID = instance->id;
  122. }
  123. rtcOccluded1M((RTCScene)instance->object,context,(RTCRay*)lrays,M,sizeof(Ray));
  124. for (size_t i=0; i<M; i++)
  125. {
  126. if (lrays[i].geomID == RTC_INVALID_GEOMETRY_ID) continue;
  127. rays[i]->geomID = 0;
  128. }
  129. }
  130. DEFINE_SET_INTERSECTOR1M(InstanceIntersector1M,FastInstanceIntersector1M);
  131. }
  132. }