lineintersection

epgpu.cu

@@ -144,6 +144,18 @@ __device__ inline int instabil_ep(const vec3& S, const vec3& C, const vec3& D) {
     return cnt;
 }
 
+__device__ inline vec3 lineintersect(const vec3& S, const vec3& A, const vec3& B) {
+    return A + (dot(S - A, B - A) * (B - A)) / dot(B - A, B - A); 
+}
+
+__device__ inline int nyereg_ep(const vec3& S, const vec3& C, const vec3& D) {
+    int cnt = 0;
+    vec3 A(0, 0, 0), B(0, 0, 1.0);
+
+
+    return cnt;
+}
+
 __device__ void ABC_oldal(int v, int w, const vec3& C, const vec3& D, char* egysulyi_mtx) {
     int pos = blockDim.x * blockIdx.x + threadIdx.x;