#include "sdl_funcs.h"
#include "objects.h"
#include "raytrace.h"
-
+#include "vector.h"
+#include "stl.h"
SDL_Event event;
+rgb colorC(const color c) {
+ rgb colors;
+ char * arr = (char*)&colors;
+ c.rgb(arr);
+ return colors;
+}
+
+double alias[2] = {.25,.5};
+int threshold[2] = {30,10}; // corners sides
+inline bool checkT(rgb ca, rgb cb,int t){
+ return ((abs(((int)ca.r)-((int)cb.r)) +abs(((int)ca.g)-((int)cb.g)) + abs(((int)ca.b)-((int)cb.b)) > t*3));
+}
+vec3 rotate(vec3 v, const vec3 k)
+{
+ double cos_theta = cos(k.x);
+ double sin_theta = sin(k.x);
+ v = (v * cos_theta) + (vec3(1,0,0).cross(v) * sin_theta) + (vec3(1,0,0) * vec3(1,0,0).dot(v)) * (1 - cos_theta);
-void render(const std::vector<Sphere> &spheres, const std::vector<Light> &lights,const Cam &cam) {
- const float fov = cam.fov;
- std::vector<vec3> framebuffer(SCREEN_WIDTH*SCREEN_HEIGHT);
+ cos_theta = cos(k.y);
+ sin_theta = sin(k.y);
+ v = (v * cos_theta) + (vec3(0,1,0).cross(v) * sin_theta) + (vec3(0,1,0) * vec3(0,1,0).dot(v)) * (1 - cos_theta);
+
+ cos_theta = cos(k.z);
+ sin_theta = sin(k.z);
+ v = (v * cos_theta) + (vec3(0,0,1).cross(v) * sin_theta) + (vec3(0,0,1) * vec3(0,0,1).dot(v)) * (1 - cos_theta);
+
+ return v;
+}
+
+void render(const Objects &spheres, const Lights &lights,const Cam &cam) {
+ const float fov = cam.fov;
+ #if ANTIALIAS
+ rgb frame1[SCREEN_HEIGHT+2][SCREEN_WIDTH+2];
+ #endif
+ frame framebuffer(SCREEN_HEIGHT,SCREEN_WIDTH);
#pragma omp parallel for
for (size_t j = 0; j<SCREEN_HEIGHT; j++) { // actual rendering loop
for (size_t i = 0; i<SCREEN_WIDTH; i++) {
- float dir_x = (i + 0.5) - SCREEN_WIDTH/2.;
- float dir_y = -(j + 0.5) + SCREEN_HEIGHT/2.; // this flips the image at the same time
- float dir_z = -SCREEN_HEIGHT/(2.*tan(fov/2.));
- framebuffer[i+j*SCREEN_WIDTH] = cast_ray(cam.pos, ((vec3{dir_x, dir_y, dir_z}+cam.dir)).normalize(), spheres, lights);
+ double dir_x = (i + 0.5) - SCREEN_WIDTH/2.;
+ double dir_y = -(j + 0.5) + SCREEN_HEIGHT/2.; // this flips the image at the same time
+ double dir_z = -SCREEN_HEIGHT/(2.*tan(fov/2.));
+ framebuffer[j][i] = colorC(cast_ray(cam.pos, rotate(vec3{dir_x, dir_y, dir_z},cam.dir).normalize(), spheres, lights));
+
+ #if ANTIALIAS
+ frame1[j][i] = framebuffer[j][i];
+ #endif
}
}
- size_t index = 0;
- char * texture_pixels = sdl_pixels_lock();
-
- for (vec3 &c : framebuffer) {
- float max = std::max(c[0], std::max(c[1], c[2]));
- if (max>1) c = c*(1./max);
- texture_pixels[index++] = (char)(255 * c[0]);
- texture_pixels[index++] = (char)(255 * c[1]);
- texture_pixels[index++] = (char)(255 * c[2]);
+ #if ANTIALIAS
+ #pragma omp parallel for
+ for (size_t j = 0; j<SCREEN_HEIGHT; j++) { // actual rendering loop
+ for (size_t i = 0; i<SCREEN_WIDTH; i++) {
+ double dir_x = (i + 0.5) - SCREEN_WIDTH/2.;
+ double dir_y = -(j + 0.5) + SCREEN_HEIGHT/2.; // this flips the image at the same time
+ double dir_z = -SCREEN_HEIGHT/(2.*tan(fov/2.));
+ if(checkT(frame1[j][i],frame1[j+1][i+1],threshold[0])|| \
+ checkT(frame1[j][i],frame1[j-1][i+1],threshold[0])|| \
+ checkT(frame1[j][i],frame1[j-1][i-1],threshold[0])|| \
+ checkT(frame1[j][i],frame1[j+1][i-1],threshold[0])|| \
+ checkT(frame1[j][i],frame1[j][i+1],threshold[1])|| \
+ checkT(frame1[j][i],frame1[j][i-1],threshold[1])|| \
+ checkT(frame1[j][i],frame1[j-1][i],threshold[1])|| \
+ checkT(frame1[j][i],frame1[j+1][i],threshold[1])){
+ //color newC(framebuffer [j][i].r/512.0,framebuffer [j][i].g/512.0,framebuffer [j][i].b/512.0);
+ color newC(0,0,0);
+ // corners
+ newC += cast_ray(cam.pos, rotate(vec3{dir_x+0.25, dir_y+0.25, dir_z},cam.dir).normalize(), spheres, lights)*alias[0];
+ newC += cast_ray(cam.pos, rotate(vec3{dir_x-0.25, dir_y+0.25, dir_z},cam.dir).normalize(), spheres, lights)*alias[0];
+ newC += cast_ray(cam.pos, rotate(vec3{dir_x-0.25, dir_y-0.25, dir_z},cam.dir).normalize(), spheres, lights)*alias[0];
+ newC += cast_ray(cam.pos, rotate(vec3{dir_x+0.25, dir_y-0.25, dir_z},cam.dir).normalize(), spheres, lights)*alias[0];
+ // sides
+ newC += cast_ray(cam.pos, rotate(vec3{dir_x, dir_y+0.25, dir_z},cam.dir).normalize(), spheres, lights)*alias[1];
+ newC += cast_ray(cam.pos, rotate(vec3{dir_x, dir_y-0.25, dir_z},cam.dir).normalize(), spheres, lights)*alias[1];
+ newC += cast_ray(cam.pos, rotate(vec3{dir_x-0.25, dir_y, dir_z},cam.dir).normalize(), spheres, lights)*alias[1];
+ newC += cast_ray(cam.pos, rotate(vec3{dir_x+0.25, dir_y, dir_z},cam.dir).normalize(), spheres, lights)*alias[1];
+ newC /= 3; // .25*4 + .5*4 + .5
+ framebuffer[j][i] = colorC(newC);
+ }
+ }
}
- sdl_pixels_unlock();
+ #endif
+}
-
-}
void signal_hand(int signum) {
std::cout << "Caught signal " << signum << std::endl;
+ #pragma omp flush // fix seg fault from open mp accesing sdl doesnt work though
+ sdl_pixels_unlock();
sdl_close(0);
}
-
-int main() {
- int winsizeX, winsizeY;
- signal(SIGINT, signal_hand);
- if( !sdl_init() )
- {
- printf( "Failed to initialize!\n" );
- exit(-1);
- }
-
- //Free resources and close SDL
+void create_objects(Objects& objects ,Lights& lights){
const Material ivory = {1.0, {0.6, 0.3, 0.1, 0.0}, {0.4, 0.4, 0.3}, 50.};
const Material glass = {1.5, {0.0, 0.5, 0.1, 0.8}, {0.6, 0.7, 0.8}, 125.};
const Material red_rubber = {1.0, {0.9, 0.1, 0.0, 0.0}, {0.3, 0.1, 0.1}, 10.};
- const Material mirror = {1.0, {0.0, 10.0, 0.8, 0.0}, {1.0, 1.0, 1.0}, 1425.};
+ const Material mirror = {1.0, {0.0, 10.0, 0.8, 0.0}, {5.0, 1.0, 1.0}, 1425.};
+ const Material wood = {1.0, {1.2, 0.1, 0.0, 0.0}, {0.2, 0.1, 0.02}, 1.};
- std::vector<Sphere> spheres = {
- Sphere{vec3{-3, 0, -16}, 2, ivory},
+ objects ={{
+ /*Sphere{vec3{-3, 0, -16}, 2, ivory},
Sphere{vec3{-1.0, -1.5, -12}, 2, glass},
Sphere{vec3{ 1.5, -0.5, -18}, 3, red_rubber},
- Sphere{vec3{ 7, 5, -18}, 4, mirror}
+ Sphere{vec3{ 7, 5, -18}, 4, mirror},
+ Sphere{vec3{-3, 10, -17}, 2, wood}*/
+ },{
+ //Triangle{vec3{-3, 0, -16},vec3{-3, 10, -17},vec3{ 7, 5, -18},red_rubber}
+ }
};
- std::vector<Light> lights = {
- {{-20, 20, 20}, 1.5},
- {{ 30, 50, -25}, 1.8},
- {{ 30, 20, 30}, 1.7}
+ lights = {{
+ {{-20, 20, 20}, {1.5,1,1}},
+ {{ 30, 50, -25}, {1,1.8,1}},
+ {{ 30, 20, 30}, {1,1,1.7}}
+ }
};
+}
+
+int main(int argc, char*argv[]) {
+ signal(SIGINT, signal_hand);
+ if( !sdl_init() )
+ {
+ printf( "Failed to initialize!\n" );
+ exit(-1);
+ }
+
Cam cam = {
- {1,1,0},
- {-1,-1,0},
+ {-1.8,12.2,16.8},
+ {-0.52,-0.30,0.},
M_PI_2
};
+ std::vector<STL_Triangle> triangles;
+ if(argc == 2)
+ triangles = parsestl(std::string(argv[1]));
+ int winsizeX, winsizeY;
+ Objects objects;
+ Lights lights;
+ create_objects(objects,lights);
+ int xM = 0,yM = 0;
-
+ const Material default_mat = {1.0, {0.9, 0.1, 0.0, 0.0}, {0.3, 0.1, 0.1}, 10.};
+ for(STL_Triangle t: triangles){
+ objects.triangle.push_back(Triangle{vec3{t.a[0],t.a[2],t.a[1]},vec3{t.c[0],t.c[2],t.c[1]},vec3{t.b[0],t.b[2],t.b[1]},default_mat});
+ }
while(1){
- int x,y;
SDL_PumpEvents();
+ const Uint8 *keystates = SDL_GetKeyboardState(NULL);
while( SDL_PollEvent( &event ) )
{
if( event.type == SDL_QUIT )
sdl_close(0);
if(event.type == SDL_MOUSEMOTION){
- SDL_GetMouseState(&x,&y);
+ //SDL_GetMouseState(&,&yM);
+ xM += event.motion.xrel*MOUSE_SENSITIVITY;
+ yM += event.motion.yrel*MOUSE_SENSITIVITY;
SDL_GetWindowSize(sdl_getwindow(), &winsizeX,&winsizeY);
- cam.dir = vec3{((float)x)*((float)SCREEN_WIDTH/(float)winsizeX),((float)-y)*((float)SCREEN_HEIGHT/(float)winsizeY),0};
- }
- if( event.type == SDL_KEYDOWN )
- {
- //Select surfaces based on key press
- switch( event.key.keysym.sym )
- {
- case SDLK_w:
- cam.pos.z -= 1;
- break;
- case SDLK_s:
- cam.pos.z += 1;
- break;
- case SDLK_a:
- cam.pos.x -= 1;
- break;
- case SDLK_d:
- cam.pos.x += 1;
- break;
- case SDLK_SPACE:
- cam.pos.y += 1;
- break;
- case SDLK_LSHIFT:
- cam.pos.y -= 1;
- break;
- default:
- break;
- }
+ cam.dir = vec3{-(yM*M_PI/winsizeX)*1.,-(xM*M_PI/winsizeY)*1,0};
}
}
- render(spheres, lights,cam);
+ if(keystates[SDL_SCANCODE_ESCAPE])
+ signal_hand(0);
+ if(keystates[SDL_SCANCODE_W])
+ cam.pos -= rotate(vec3{0,0,MOVEMENT_SPEED},vec3{-(yM*M_PI/winsizeX)*1.,-(xM*M_PI/winsizeY)*1.,0});
+ if(keystates[SDL_SCANCODE_S])
+ cam.pos += rotate(vec3{0,0,MOVEMENT_SPEED},vec3{-(yM*M_PI/winsizeX)*1.,-(xM*M_PI/winsizeY)*1.,0});
+ if(keystates[SDL_SCANCODE_A])
+ cam.pos -= rotate(vec3{MOVEMENT_SPEED,0,0},vec3{-(yM*M_PI/winsizeX)*1.,-(xM*M_PI/winsizeY)*1.,0});
+ if(keystates[SDL_SCANCODE_D])
+ cam.pos += rotate(vec3{MOVEMENT_SPEED,0,0},vec3{-(yM*M_PI/winsizeX)*1.,-(xM*M_PI/winsizeY)*1.,0});
+ if(keystates[SDL_SCANCODE_SPACE])
+ cam.pos.y += MOVEMENT_SPEED;
+ if(keystates[SDL_SCANCODE_LSHIFT])
+ cam.pos.y -= MOVEMENT_SPEED;
+ render(objects, lights,cam);
+ //SDL_Log("{%lf,%lf,%lf},{%lf,%lf,%lf}",cam.pos.x,cam.pos.y,cam.pos.z,cam.dir.x,cam.dir.y,cam.dir.z);
sdl_frame();
}
sdl_free();
projects / raytracer.git / blobdiff