diff --git a/src/main.cpp b/src/main.cpp
index 88dfbe4..446191f 100644
--- a/src/main.cpp
+++ b/src/main.cpp
@@ -2,6 +2,7 @@
 #include "scene/Scene.h"
 #include "util/ModelLoader.h"
 #include "window/Window.h"
+#include "util/ModelLoader.h"
 
 int main()
 {
@@ -17,6 +18,7 @@ int main()
             .height = 1080,
             .numSamples = 10000,
         });
+
     while (true)
     {
         window.update(scene.getImage());
diff --git a/src/scene/AABB.h b/src/scene/AABB.h
index f97e003..aa7ef05 100644
--- a/src/scene/AABB.h
+++ b/src/scene/AABB.h
@@ -3,6 +3,7 @@
 #include <glm/glm.hpp>
 #include <algorithm>
 #include "util/Ray.h"
+#include <iostream>
 
 struct AABB
 {
@@ -45,9 +46,13 @@ struct AABB
     }
     bool intersects(Ray ray, float tmin, float tmax)
     {
+      std::cout << "Hello1" << std::endl;
         glm::vec3 invD = 1.0f / ray.direction;
+      std::cout << "Hello2: " << min.x << " " << min.y << " " << min.z << " " << max.x << " " << max.y << " " << max.z << " " << std::endl;
         glm::vec3 t0s = glm::vec3(min - ray.origin) * invD;
+      std::cout << "Hello3" << std::endl;
         glm::vec3 t1s = glm::vec3(max - ray.origin) * invD;
+      std::cout << "Hello4" << std::endl;
 
         glm::vec3 tsmaller = glm::min(t0s, t1s);
         glm::vec3 tbigger = glm::max(t0s, t1s);
diff --git a/src/scene/BVH.cpp b/src/scene/BVH.cpp
index e754f9c..d7d5dee 100644
--- a/src/scene/BVH.cpp
+++ b/src/scene/BVH.cpp
@@ -1,6 +1,7 @@
 #include "BVH.h"
 #include <algorithm>
 #include <ranges>
+#include <iostream>
 
 void BVH::addModel(PModel model, glm::mat4 transform)
 {
@@ -10,10 +11,10 @@ void BVH::addModel(PModel model, glm::mat4 transform)
 
 void BVH::addModels(std::vector<PModel> _models, glm::mat4 transform)
 {
-    for (int i = 0; i < _models.size(); ++i)
+    for (auto& _model : _models)
     {
-        _models[i]->boundingBox.transform(transform);
-        models.push_back(std::move(_models[i]));
+        _model->boundingBox.transform(transform);
+        models.push_back(std::move(_model));
     }
 }
 
@@ -60,13 +61,13 @@ std::optional<IntersectionInfo> BVH::traceRay(Ray ray)
 {
     auto results = generateIntersections(hierarchy, ray);
     float closestT = std::numeric_limits<float>::max();
-    IntersectionInfo info;
-    for (uint32_t i = 0; i < results.size(); ++i)
+    std::optional<IntersectionInfo> info = {};
+    for (auto& result : results)
     {
-        if (results[i].t < closestT)
+        if (result.t < closestT)
         {
-            closestT = results[i].t;
-            info = results[i];
+            closestT = result.t;
+            info = result;
         }
     }
     if (closestT < std::numeric_limits<float>::max())
diff --git a/src/scene/Scene.cpp b/src/scene/Scene.cpp
index 93d0efc..dc488fc 100644
--- a/src/scene/Scene.cpp
+++ b/src/scene/Scene.cpp
@@ -1,8 +1,16 @@
 #include "Scene.h"
+#include "util/ModelLoader.h"
 #include <iostream>
 #include <chrono>
 
-Scene::Scene() {}
+Scene::Scene()
+{
+  bvh.addModels(ModelLoader::loadModel("../cube.fbx"),
+               glm::mat4(glm::vec4(1.0f, 0.0f, 0.0f, 0.0f),
+                         glm::vec4(0.0f, 1.0f, 0.0f, 0.0f),
+                         glm::vec4(0.0f, 0.0f, 1.0f, 0.0f),
+                         glm::vec4(0.0f, 0.0f, 0.0f, 1.0f)));
+}
 
 Scene::~Scene() {}
 
@@ -35,10 +43,16 @@ void Scene::render(Camera cam, RenderParameter params)
                         {
                             for (int h = 0; h < params.height; ++h)
                             {
-                                // Ray r = Ray();
-                                // bvh.traceRay(r);
-                                accumulator[w + h * params.width] +=
-                                    glm::vec3(w / float(params.width * params.numSamples), h / float(params.height * params.numSamples), 0);
+                              Ray r = Ray{
+                                .origin = glm::vec3(0.0f),
+                                .direction = glm::vec3(0.0f, 0.0f, 1.0f)
+                              };
+                              bvh.traceRay(r);
+
+                              std::cout << "HELLO" << std::endl;
+
+                              accumulator[w + h * params.width] +=
+                                  glm::vec3(w / float(params.width * params.numSamples), h / float(params.height * params.numSamples), 0);
                             }
                         });
                 }
diff --git a/src/util/ModelLoader.cpp b/src/util/ModelLoader.cpp
index c2330f8..02a8954 100644
--- a/src/util/ModelLoader.cpp
+++ b/src/util/ModelLoader.cpp
@@ -28,8 +28,8 @@ std::vector<PModel> ModelLoader::loadModel(std::string_view filename)
             model->indices.push_back(face.mIndices[1]);
             model->indices.push_back(face.mIndices[2]);
 
-            auto e0 = model->positions[face.mIndices[1] - face.mIndices[0]];
-            auto e1 = model->positions[face.mIndices[2] - face.mIndices[0]];
+            auto e0 = model->positions[face.mIndices[1]] - model->positions[face.mIndices[0]];
+            auto e1 = model->positions[face.mIndices[2]] - model->positions[face.mIndices[0]];
             model->es.push_back(e0);
             model->es.push_back(e1);