From da3889d4706b982da33dd1416ee3ee6961adb602 Mon Sep 17 00:00:00 2001
From: Alex <shadow8t4@tamu.edu>
Date: Sun, 27 Nov 2016 01:12:12 -0600
Subject: [PATCH] Template .obj now duplicates on modular scale with modular
 spacing between layers. Will clean up comments later.

---
 source/main.cpp | 108 ++++++++++++++++++++++++++++++++++++++++++++----
 1 file changed, 101 insertions(+), 7 deletions(-)

diff --git a/source/main.cpp b/source/main.cpp
index 0dbe6bc..09456fd 100644
--- a/source/main.cpp
+++ b/source/main.cpp
@@ -7,6 +7,9 @@
 #include <vector>
 #include "../include/CompFab.h"
 #include "../include/Mesh.h"
+#include <math.h>
+
+#define PI 3.14159265
 
 void findLW(Mesh &m, double &l, double &w)
 {
@@ -23,38 +26,129 @@ void findLW(Mesh &m, double &l, double &w)
         }
         if(m.v[i].m_y < minw)
         {
-            minw = m.v[i].m_x;
+            minw = m.v[i].m_y;
         }
         if(m.v[i].m_y > maxw)
         {
-            maxw = m.v[i].m_x;
+            maxw = m.v[i].m_y;
         }
     }
     l = maxl - minl;
     w = maxw - minw;
 }
 
+// Calculate translation matrices and add them to the vector of translation matrices provided.
+std::vector<CompFab::Vec3> createVec3d(int layers, double spacing, double length, double width)
+{
+    std::vector<CompFab::Vec3> *output = new std::vector<CompFab::Vec3>();
+    
+    CompFab::Vec3 *temp = new CompFab::Vec3(-(length + spacing), -(width + spacing), 0);
+    CompFab::Vec3 *trans = new CompFab::Vec3(0, spacing, 0);
+    
+    for(int cl = 1; cl < layers; cl++)
+    {
+        *temp = CompFab::Vec3(-(length + spacing)*cl, -(width + spacing)*cl, 0);
+        
+        for(int c = 0; c < cl*8; c++)
+        {
+            /*
+            double cosine = cos(((2*PI)/(cl*8))*c);
+            double sine = sin(((2*PI)/(cl*8))*c);
+            //
+            double xcoord = 1;
+            double angle = tan((2*PI)/(cl*8))*c;
+            double ycoord = 1;
+            if(cosine != 0)
+            {
+                ycoord = sine/cosine;
+            }
+            if(sine != 0)
+            {
+                xcoord = cosine/sine;
+            }
+            //
+            double x = cosine;//(sqrt(1 - (sine*sine)/2));
+            if(x > 0)
+            {
+                x = floor(x);
+            }
+            else
+            {
+                x = ceil(x);
+            }
+            double y = sine;//(sqrt(1-(cosine*cosine)/2));
+            if(y > 0)
+            {
+                y = floor(x);
+            }
+            else
+            {
+                y = ceil(x);
+            }
+            */
+            *trans = CompFab::Vec3(spacing*cos(floor((c/(2*cl)))*(0.5*PI)) + length*cos(floor((c/(2*cl)))*(0.5*PI)), spacing*sin(floor((c/(2*cl)))*(0.5*PI)) + width*sin(floor((c/(2*cl)))*(0.5*PI)), 0);
+            *temp = *temp + *trans;
+            
+            output->push_back(*temp);
+        }
+    }
+    
+    /*
+    temp->push_back(CompFab::Vec3(length,0,0));
+    temp->push_back(CompFab::Vec3(length,width,0));
+    temp->push_back(CompFab::Vec3(0,width,0));
+    temp->push_back(CompFab::Vec3(-length,width,0));
+    temp->push_back(CompFab::Vec3(-length,0,0));
+    temp->push_back(CompFab::Vec3(-length,-width,0));
+    temp->push_back(CompFab::Vec3(0,-width,0));
+    temp->push_back(CompFab::Vec3(length,-width,0));
+    */
+    
+    return *output;
+}
+
 int main(int argc, char **argv)
 {
-    unsigned int num = 16; //number of voxels (e.g. 16x16x16)
     
+    // Error checking.
     if(argc < 3)
     {
         std::cout << "Usage: [executable] [template].obj output.obj [optional: -d for debugging output]" << std::endl;
         std::exit(1);
     }
+    
+    // Modularize this later.
+    double spacing = 1.0;
+    
+    // Create Mesh object from file, output to manipulate from template Mesh.
     Mesh *test = new Mesh(argv[1], false);
     Mesh *output = new Mesh(test->v, test->t);
+    
+    int layers = 10;
     double l = 0, w = 0;
     double *length = &l, *width = &w;
+    
+    // Find dimensions for the mesh. Assumes the mesh is facing upright.
     findLW(*test, *length, *width);
-    for(int i = 0; i < test->v.size(); i++)
+    
+    std::vector<CompFab::Vec3> d = createVec3d(layers, spacing, *length, *width);
+    
+    // Duplicating template, will later be replaced with a much more robust procedural generation function.
+    for(int i = 0; i < d.size(); i++)
     {
-        output->v.push_back(*new CompFab::Vec3(test->v[i].m_x + *length + 1, test->v[i].m_y, test->v[i].m_z));
+        for(int j = 0; j < test->v.size(); j++)
+        {
+            output->v.push_back(CompFab::Vec3(test->v[j] + d[i]));
+        }
     }
-    for(int k = 0; k < test->t.size(); k++)
+    
+    // Copying needed triangle data.
+    for(int n = 1; n < (2*layers - 1)*(2*layers - 1); n++)
     {
-        output->t.push_back(*new CompFab::Vec3i(test->t[k].m_x + test->v.size(), test->t[k].m_y + test->v.size(), test->t[k].m_z + test->v.size()));
+        for(int k = 0; k < test->t.size(); k++)
+        {
+            output->t.push_back(CompFab::Vec3i(test->t[k].m_x + test->v.size()*n, test->t[k].m_y + test->v.size()*n, test->t[k].m_z + test->v.size()*n));
+        }
     }
     
     // Debugging