Initial commit

Author: Kugelstadt

AVX_math.h

@@ -0,0 +1,72 @@
+cmake_minimum_required(VERSION 3.1)
+
+project(FastCorotDemo)
+
+set(CMAKE_MODULE_PATH "${CMAKE_SOURCE_DIR}/cmake/Modules" ${CMAKE_MODULE_PATH})
+set(PROJECT_PATH ${PROJECT_SOURCE_DIR})
+include_directories(${PROJECT_SOURCE_DIR})
+
+set(CMAKE_RUNTIME_OUTPUT_DIRECTORY "${CMAKE_SOURCE_DIR}/bin")
+set(CMAKE_RUNTIME_OUTPUT_DIRECTORY_DEBUG "${CMAKE_SOURCE_DIR}/bin")
+set(CMAKE_RUNTIME_OUTPUT_DIRECTORY_RELEASE "${CMAKE_SOURCE_DIR}/bin")
+set(CMAKE_RUNTIME_OUTPUT_DIRECTORY_RELWITHDEBINFO "${CMAKE_SOURCE_DIR}/bin")
+set(CMAKE_RUNTIME_OUTPUT_DIRECTORY_MINSIZEREL "${CMAKE_SOURCE_DIR}/bin")
+set(LIBRARY_OUTPUT_PATH ${CMAKE_SOURCE_DIR}/lib)
+set(CMAKE_DEBUG_POSTFIX "_d")
+set(CMAKE_RELWITHDEBINFO_POSTFIX "_rd")
+set(CMAKE_MINSIZEREL_POSTFIX "_ms")
+
+set(SOURCES
+	main.cpp
+	FastCorotFEM.cpp
+	TetModel.cpp
+	utilities/TetGenLoader.cpp
+	utilities/MiniGL.cpp
+)
+
+set(HEADERS
+	AVX_math.h
+	FastCorotFEM.h
+	Common.h
+	utilities/TetGenLoader.h
+	utilities/MiniGL.h
+	utilities/Timing.h
+	TetModel.h
+)
+
+if ("${CMAKE_CXX_COMPILER_ID}" STREQUAL "MSVC")
+	add_compile_options(/arch:AVX)
+elseif ("${CMAKE_CXX_COMPILER_ID}" STREQUAL "GNU")
+	add_compile_options("-mavx")
+endif()
+
+set(EIGEN3_INCLUDE_DIR "${CMAKE_SOURCE_DIR}/extern/eigen")
+find_package( Eigen3 REQUIRED )
+
+# Set include directories.
+include_directories(
+    ${EIGEN3_INCLUDE_DIR}
+)
+
+# OpenGL & GLEW library
+find_package(OpenGL)
+
+# executable
+add_executable(FastCorotDemo
+	${HEADERS}
+    ${SOURCES}
+	)
+	
+if (WIN32)
+  subdirs(extern/freeglut)
+  include_directories(${PROJECT_PATH}/extern/freeglut/include)
+  
+  set(GLUT_LIBRARIES freeglut)
+  add_dependencies(FastCorotDemo freeglut)
+else()
+  find_package(GLUT REQUIRED)
+endif()	
+	
+target_link_libraries(FastCorotDemo ${OPENGL_LIBRARIES} ${GLUT_LIBRARIES})
+
+add_definitions(-DDATA_PATH="${PROJECT_PATH}/meshes/")

CMakeLists.txt

@@ -0,0 +1,108 @@
+#ifndef COMMON_H
+#define COMMON_H
+
+#include <Eigen/Dense>
+#include "float.h"
+#define _USE_MATH_DEFINES
+#include <cmath> 
+
+
+#define USE_DOUBLE
+#define MIN_PARALLEL_SIZE 64
+
+#ifdef USE_DOUBLE
+typedef double Real;
+
+#define REAL_MAX DBL_MAX
+#define REAL_MIN DBL_MIN
+#else
+typedef float Real;
+
+#define REAL_MAX FLT_MAX
+#define REAL_MIN FLT_MIN
+#endif
+
+using Vector2r = Eigen::Matrix<Real, 2, 1>;
+using Vector3r = Eigen::Matrix<Real, 3, 1>;
+using Vector4r = Eigen::Matrix<Real, 4, 1>;
+using Vector6r = Eigen::Matrix<Real, 6, 1>;
+using Vector9r = Eigen::Matrix<Real, 9, 1>;
+using VectorXr = Eigen::Matrix<Real, Eigen::Dynamic, 1>;
+using Matrix2r = Eigen::Matrix<Real, 2, 2>;
+using Matrix3r = Eigen::Matrix<Real, 3, 3>;
+using Matrix4r = Eigen::Matrix<Real, 4, 4>;
+using Matrix6r = Eigen::Matrix<Real, 6, 6>;
+using Matrix9r = Eigen::Matrix<Real, 9, 9>;
+using MatrixXr = Eigen::Matrix<Real, Eigen::Dynamic, Eigen::Dynamic>;
+using AlignedBox2r = Eigen::AlignedBox<Real, 2>;
+using AlignedBox3r = Eigen::AlignedBox<Real, 3>;
+using AngleAxisr = Eigen::AngleAxis<Real>;
+using Quaternionr = Eigen::Quaternion<Real>;
+
+//allocators to be used in STL collections containing Eigen structures
+using Alloc_Vector2r = Eigen::aligned_allocator<Vector2r>;
+using Alloc_Vector3r = Eigen::aligned_allocator<Vector3r>;
+using Alloc_Vector4r = Eigen::aligned_allocator<Vector4r>;
+using Alloc_Matrix2r = Eigen::aligned_allocator<Matrix2r>;
+using Alloc_Matrix3r = Eigen::aligned_allocator<Matrix3r>;
+using Alloc_Matrix4r = Eigen::aligned_allocator<Matrix4r>;
+using Alloc_AlignedBox2r = Eigen::aligned_allocator<AlignedBox2r>;
+using Alloc_AlignedBox3r = Eigen::aligned_allocator<AlignedBox3r>;
+using Alloc_AngleAxisr = Eigen::aligned_allocator<AngleAxisr>;
+using Alloc_Quaternionr = Eigen::aligned_allocator<Quaternionr>;
+
+#if EIGEN_ALIGN
+	#define PDB_MAKE_ALIGNED_OPERATOR_NEW EIGEN_MAKE_ALIGNED_OPERATOR_NEW
+	#define REPORT_MEMORY_LEAKS
+
+#if defined(WIN32) || defined(_WIN32) || defined(WIN64)	   
+#ifdef _DEBUG
+	// Enable memory leak detection for Eigen new
+	#undef PDB_MAKE_ALIGNED_OPERATOR_NEW
+	#define PDB_MAKE_ALIGNED_OPERATOR_NEW	EIGEN_MAKE_ALIGNED_OPERATOR_NEW \
+		void *operator new(size_t size, int const block_use, char const*  file_name, int const line_number) { \
+		\
+			return _aligned_malloc_dbg(size, 16, file_name, line_number); \
+		} \
+		void *operator new[](size_t size, int const block_use, char const*  file_name, int const line_number) { \
+			return operator new(size, block_use, file_name, line_number); \
+		}\
+		void operator delete(void* block, int const block_use, char const*  file_name, int const line_number) { \
+		\
+			return _aligned_free_dbg(block); \
+		} \
+		void operator delete[](void* block, int const block_use, char const*  file_name, int const line_number) { \
+			return operator delete(block, block_use, file_name, line_number); \
+		}	
+#endif
+#endif
+#else
+	#define PDB_MAKE_ALIGNED_OPERATOR_NEW
+
+#if defined(WIN32) || defined(_WIN32) || defined(WIN64)	   
+	// Enable memory leak detection
+#ifdef _DEBUG
+	#define _CRTDBG_MAP_ALLOC 
+	#include <stdlib.h> 
+	#include <crtdbg.h> 
+	#define DEBUG_NEW new(_NORMAL_BLOCK, __FILE__, __LINE__) 	
+	#define REPORT_MEMORY_LEAKS _CrtSetDbgFlag ( _CRTDBG_ALLOC_MEM_DF | _CRTDBG_LEAK_CHECK_DF );
+#else
+	#define REPORT_MEMORY_LEAKS
+#endif
+#else
+	#define REPORT_MEMORY_LEAKS
+#endif
+
+#endif
+
+
+#endif
+
+#if defined(WIN32) || defined(_WIN32) || defined(WIN64)	   
+#define FORCE_INLINE __forceinline
+#else
+#define FORCE_INLINE __attribute__((always_inline))
+#endif
+
+