File: MotionTransform.cpp

package info (click to toggle)
ospray 3.2.0-2
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid
  • size: 10,048 kB
  • sloc: cpp: 80,569; ansic: 951; sh: 805; makefile: 170; python: 69
file content (95 lines) | stat: -rw-r--r-- 3,258 bytes parent folder | download | duplicates (2)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
// Copyright 2021 Intel Corporation
// SPDX-License-Identifier: Apache-2.0

// ospray
#include "MotionTransform.h"

namespace ospray {

void MotionTransform::readParams(ISPCDeviceObject &obj)
{
  transforms = obj.getParamDataT<affine3f>("motion.transform");
  if (transforms) {
    motionBlur = transforms->size() > 1;
    quaternion = false;
    transform = (*transforms)[0];
  } else { // look for SRT
    auto shift = obj.getParamDataT<vec3f>("motion.pivot");
    auto scale = obj.getParamDataT<vec3f>("motion.scale");
    auto rotation = obj.getParamDataT<quatf>("motion.rotation");
    auto translation = obj.getParamDataT<vec3f>("motion.translation");
    quaternion = shift || scale || rotation || translation;
    if (quaternion) {
      // determine (minimum) size
      size_t size = std::numeric_limits<size_t>::max();
      if (shift)
        size = shift->size();
      if (scale)
        size = std::min(size, scale->size());
      if (rotation)
        size = std::min(size, rotation->size());
      if (translation)
        size = std::min(size, translation->size());
      motionBlur = size > 1;
      quaternions.resize(size);
      for (size_t i = 0; i < size; i++) {
        auto *qdecomp = &quaternions[i];
        rtcInitQuaternionDecomposition(qdecomp);
        if (shift) {
          const vec3f &s = (*shift)[i];
          rtcQuaternionDecompositionSetShift(qdecomp, s.x, s.y, s.z);
        }
        if (scale) {
          const vec3f &s = (*scale)[i];
          rtcQuaternionDecompositionSetScale(qdecomp, s.x, s.y, s.z);
        }
        if (rotation) {
          const quatf &q = (*rotation)[i];
          rtcQuaternionDecompositionSetQuaternion(qdecomp, q.r, q.i, q.j, q.k);
        }
        if (translation) {
          const vec3f &t = (*translation)[i];
          rtcQuaternionDecompositionSetTranslation(qdecomp, t.x, t.y, t.z);
        }
      }
    } else { // add single transform
      quaternions.clear();
      motionBlur = false;
      transform = obj.getParam<affine3f>(
          "transform", obj.getParam<affine3f>("xfm", affine3f(one)));
    }
  }
  time = range1f(0.0f, 1.0f);
  if (motionBlur) {
    time = obj.getParam<range1f>("time", range1f(0.0f, 1.0f));
    if (time.upper < time.lower)
      time.upper = time.lower;
  }
}

void MotionTransform::setEmbreeTransform(RTCGeometry embreeGeometry) const
{
  if (quaternion) {
    rtcSetGeometryTimeStepCount(embreeGeometry, quaternions.size());
    for (size_t i = 0; i < quaternions.size(); i++)
      rtcSetGeometryTransformQuaternion(embreeGeometry, i, &quaternions[i]);
    rtcSetGeometryTimeRange(embreeGeometry, time.lower, time.upper);
  } else {
    if (motionBlur) {
      rtcSetGeometryTimeStepCount(embreeGeometry, transforms->size());
      for (size_t i = 0; i < transforms->size(); i++)
        rtcSetGeometryTransform(embreeGeometry,
            i,
            RTC_FORMAT_FLOAT3X4_COLUMN_MAJOR,
            &(*transforms)[i]);
      rtcSetGeometryTimeRange(embreeGeometry, time.lower, time.upper);
    } else {
      rtcSetGeometryTimeStepCount(embreeGeometry, 1);
      rtcSetGeometryTransform(
          embreeGeometry, 0, RTC_FORMAT_FLOAT3X4_COLUMN_MAJOR, transform);
    }
  }
  rtcCommitGeometry(embreeGeometry);
}

} // namespace ospray