From d769c253e8a6ddc67cf8424da2017d309e93f11b Mon Sep 17 00:00:00 2001
From: amy <ameliasquires@disroot.org>
Date: Fri, 14 Apr 2023 13:22:32 +0000
Subject: init

---
 src/point.c | 212 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 212 insertions(+)
 create mode 100644 src/point.c

(limited to 'src/point.c')

diff --git a/src/point.c b/src/point.c
new file mode 100644
index 0000000..d8e4c56
--- /dev/null
+++ b/src/point.c
@@ -0,0 +1,212 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <math.h>
+#include "glfww.h"
+#include <unistd.h>
+typedef struct {
+  double x;
+  double y;
+  double z;
+} cord;
+typedef struct {
+  int roll; //x axis
+  int pitch; //y axis
+  int yaw; //z axis
+} rot;
+typedef struct {
+  float opacity;
+} opts;
+typedef struct {
+  cord loc;
+  rot rot;
+  int uid;
+  int id;
+  opts opts;
+} point;
+
+double binominal(int n, int k){
+  if(n==k)
+    return 1.0;
+  double v = 1.0; 
+  for(int i = 1; i<=k; i++){
+    v=v*((float)(n+1-i)/i);
+  } 
+  return v;
+}
+int ma = 4;
+void mul_matr(double a[][ma],double b[][ma], double res[][ma]){
+  for(int i = 0; i<ma; i++){
+    for(int j = 0; j<ma;j++){
+      res[i][j] = 0;
+      for(int k = 0; k<ma;k++){
+        res[i][j]+=a[i][k]*b[k][j];
+        
+      }
+    }
+  }
+}
+void perspective_proj(GLFWwindow* b,cord* c, int len,double ctx,double cty,double ctz,double cx, double cy, double cz){
+  //double cx = -100;
+  //double cy = 0;
+  //double cz = 0;
+  //double ctx = 0;
+  //double cty = 0;
+  //double ctz = 0;
+  //cz=100-cz;
+  double coy = cos(cty);
+  double siz = sin(ctz);
+  double coz = cos(ctz);
+  double six = sin(ctx);
+  double siy = sin(cty);
+  double cox = cos(ctx);
+  double fov = .1;
+  double ex = cx;
+  double ey = cy;
+  //glfwGetFramebufferSize(b,&w,&h);
+  refresh_size(b);
+  glColor3f(1.0f,0.0f,0.0f);
+  double ez = 1/tan(fov/2);
+  glBegin(GL_POINTS);
+  for(int i = 0; i!=len; i++){
+    double ax = c[i].x;
+    double ay = c[i].y;
+    double az = c[i].z;
+    
+    double eyz = (coz*(ay-cy)-siz*ax-cx);
+    double yzm = (coy*(az-cz) + siy*(siz*(ay-cy) + coz*(ax-cx)));
+    double dx = coy * (siz*(ay-cy) + coz*(ax-cx)) - (siy*(az-cz));
+    double dy = six * yzm + cox*eyz;
+    double dz = cox * yzm - six*eyz; 
+    
+    //printf("%f %f %f, %f %f\n",dx,dy,ez,ez/dz*dx+dx,dy*ez);
+    double bx = ez/dz*dx+dx;
+    double by = ez/dz*dy+dy;
+    //printf("%f %f | %f %f %f\n",bx,by,ctx,cty,ctz);
+    //int aaa = round((150-dz)/2);
+    if(dz>=0)
+    glfw_circle_partial(b,nearbyint(bx),nearbyint(by),/*(0>=aaa?1:*/1/*)*/); 
+    //glfw_pixel_partial(b,round(bx),round(by)); 
+  }
+  glEnd();
+}
+cord* basier3d(double*xx,double*yy,double*zz,int n){
+  cord* pa = malloc(sizeof(*pa)*(3000));
+  //double xx[5] = {5.0,5.0,50.0,5.0,5.0};
+  //double yy[5] = {5.0,5.0,5.0,5.0,10.0};
+  //double zz[5] = {10.0,5.0,5.0,5.0,5.0};
+  //int n = 5-1;
+  n-=1;
+  for(int iy = 0; iy<=100;iy+=1){
+    double t = 0.01*iy;
+    double bcx = 0;
+    double bcy = 0;
+    double bcz = 0;
+    for(int i = 0; i <=n;i++){
+      double pp = binominal(n,i) * pow((1 - t),(n - i)) * pow(t,i);
+      bcx +=  pp * xx[i];
+      bcy +=  pp * yy[i];
+      bcz +=  pp * zz[i];
+      //if(ii==3){
+      //  printf("pp: %f bi : %f p1 : %f p2 : %f| %f, %f, %f\n",pp,binominal(n,i),bcx,bcy,bcz,pow((1-t),(n-i)),pow(t,i));
+      //}
+    }
+    //int ii = floor(100*(double)(t/1.0));
+    //printf("%i\n",iy);
+    pa[iy].x = bcx;
+    pa[iy].y = bcy;
+    pa[iy].z = bcz;
+    //printf("(%f,%f,%f)\n",bcx,bcy,bcz);
+  }
+  //bw b = sdl_init();
+  /*for(int i = 0; i!=100; i ++){
+    SDL_SetRenderDrawColor(b.r,0,255-255*(pa[i].z/100),0,255);
+    SDL_RenderDrawPoint(b.r,pa[i].x,pa[i].y);
+    sdl_circle(b,round(pa[i].x),round(pa[i].y),round(100-pa[i].z));
+    //printf("%i : (%f,%f,%f)\n",i,pa[i].x,pa[i].y,(100-pa[i].z));
+  }*/
+  //perspective_proj(b,pa,100);
+  //sdl_loop(b);
+  return pa;
+}
+void join_cords(cord* a, cord* b, int a_len, int b_len){
+  a = realloc(a,sizeof(*a)*(a_len+b_len+2));
+  for(int i = 0; i<=a_len+b_len; i++){
+    a[a_len+i] = b[i]; 
+  }
+}
+int main(){
+  double xx[5] = {5.0,5.0,5.0};
+  double yy[5] = {5.0,100.0,200.0};
+  double zz[5] = {95.0,95.0,95.0};
+  cord* a = basier3d(xx,yy,zz,3);
+ 
+  double xx2[5] = {200.0,200.0,200.0};
+  double yy2[5] = {5.0,100.0,200.0};
+  double zz2[5] = {95.0,95.0,95.0};
+  cord* b = basier3d(xx2,yy2,zz2,3);
+  
+  double xx3[5] = {5.0,100.0,200.0};
+  double yy3[5] = {5.0,5.0,5.0};
+  double zz3[5] = {95.0,95.0,95.0};
+  cord* c = basier3d(xx3,yy3,zz3,3);
+  
+  double xx4[5] = {5.0,100.0,200.0};
+  double yy4[5] = {200.0,200.0,200.0};
+  double zz4[5] = {95.0,95.0,95.0};
+  cord* d = basier3d(xx4,yy4,zz4,3); 
+  
+  double xx5[5] = {5.0,100.0,200.0};
+  double yy5[5] = {200.0,200.0,200.0};
+  double zz5[5] = {150.0,150.0,150.0};
+  cord* e = basier3d(xx5,yy5,zz5,3); 
+  
+  double xx6[5] = {5.0,100.0,200.0};
+  double yy6[5] = {5.0,5.0,5.0};
+  double zz6[5] = {150.0,150.0,150.0};
+  cord* f = basier3d(xx6,yy6,zz6,3);
+ 
+  double xx7[5] = {200.0,200.0,200.0};
+  double yy7[5] = {5.0,100.0,200.0};
+  double zz7[5] = {150.0,150.0,150.0};
+  cord* g = basier3d(xx7,yy7,zz7,3);
+  
+  double xx8[5] = {5.0,5.0,5.0};
+  double yy8[5] = {5.0,100.0,200.0};
+  double zz8[5] = {150.0,150.0,150.0};
+  cord* h = basier3d(xx8,yy8,zz8,3);
+  join_cords(a,b,100,100);
+  join_cords(a,c,200,100);
+  join_cords(a,d,300,100);
+  join_cords(a,e,400,100);
+  join_cords(a,f,500,100);
+  join_cords(a,g,600,100);
+  join_cords(a,h,700,100);
+  GLFWwindow* w = glfw_init();
+  for(double rr = 0.01;rr<=20;rr+=0.01){
+    double p1 = 0; 
+    double p2 = rr;
+    double p3 = 0;
+    double p4 = -100;
+    double p5 = 0;
+    double p6 = rr*100*0;
+    perspective_proj(w,a,800,p1,p2,p3,p4,p5,p6);
+    //perspective_proj(w,b,100,p1,p2,p3,p4,p5,p6);
+    //perspective_proj(w,c,100,p1,p2,p3,p4,p5,p6);
+    //perspective_proj(w,d,100,p1,p2,p3,p4,p5,p6);
+    //perspective_proj(w,e,100,p1,p2,p3,p4,p5,p6);
+    //perspective_proj(w,f,100,p1,p2,p3,p4,p5,p6);
+    //perspective_proj(w,g,100,p1,p2,p3,p4,p5,p6);
+    //perspective_proj(w,h,100,p1,p2,p3,p4,p5,p6);
+    glfw_load(w);
+    usleep(500);
+    glfw_clear(w);
+    if(glfwWindowShouldClose(w))break;
+  }
+  free(a);
+  free(b);
+  free(c);
+  free(d);
+  glfwDestroyWindow(w);
+  win_clean();
+  return 0;
+}
-- 
cgit v1.2.3