path: root/src/render/model.h

#include "../util/typedef.h"
#include "../util/string.h"
#include "../util/vector.h"
#include "../util/file.h"

#include <assert.h>

struct MODEL_DATA {
  VEC3 position;
  VEC2 uv;
  VEC3 normal;
};

struct MODEL {
  LIST<MODEL_DATA> vertices;
  LIST<U32> indices;
};

enum LINE_TYPES : I32 {
  comment      = ( '#' | ( ' ' << 8 ) ),
  transparency = ( 'd' | ( ' ' << 8 ) ),
  face         = ( 'f' | ( ' ' << 8 ) ),
  group        = ( 'g' | ( ' ' << 8 ) ),
  illum        = ( 'i' | ( 'l' << 8 ) ),
  ambient_clr  = ( 'K' | ( 'a' << 8 ) ),
  diffuse_clr  = ( 'K' | ( 'd' << 8 ) ),
  emissive_clr = ( 'K' | ( 'e' << 8 ) ),
  specular_clr = ( 'K' | ( 's' << 8 ) ),
  texname      = ( 'm' | ( 'a' << 8 ) ),
  mtlname      = ( 'm' | ( 't' << 8 ) ),
  newmtl       = ( 'n' | ( 'e' << 8 ) ),
  opticdensity = ( 'N' | ( 'i' << 8 ) ),
  specular_exp = ( 'N' | ( 's' << 8 ) ),
  object       = ( 'o' | ( ' ' << 8 ) ),
  shading      = ( 's' | ( ' ' << 8 ) ),
  vert         = ( 'v' | ( ' ' << 8 ) ),
  vertnorm     = ( 'v' | ( 'n' << 8 ) ),
  vertuv       = ( 'v' | ( 't' << 8 ) ),
  usemtl       = ( 'u' | ( 's' << 8 ) )
};

static inline char* skip_space( char* s ) {
  while( *s == ' ' || *s == '\t' ) ++s;
  return s;
}

static inline VEC3 read_vec3( char* s ) {
  VEC3 temp{};
  s = skip_space( s );
  temp.x = strtof( s, &s );
  s = skip_space( s );
  temp.y = strtof( s, &s );
  s = skip_space( s );
  temp.z = strtof( s, &s );
  return temp;
}

static inline LINE_TYPES get_first_two( char* s ) {
  I32 c0 = *s;
  U16 c1 = *(s + 1);
  return (LINE_TYPES)( c0 | c1 << 8 );
}

struct MTLDATA {
  STR name;
  F32 specular_exp;
  VEC3 ambient_clr;
  VEC3 diffuse_clr;
  VEC3 specular_clr;
  VEC3 emissive_clr;
  F32 opticdensity;
  F32 transparency;
  I32 illum;
  STR map_Kd;
};

struct MTLLIB {
  STR name;
  LIST<MTLDATA> data;
};

inline I32 line_type_offset( LINE_TYPES first_two ) {
  switch( first_two ) {
  case face:
  case vert:
  case group:
  case object:
  case shading:
  case transparency: return 2;
  case ambient_clr:
  case diffuse_clr:
  case emissive_clr:
  case specular_clr:
  case opticdensity:
  case specular_exp:
  case vertuv:
  case vertnorm: return 3;
  case illum: return 6;
  case mtlname:
  case newmtl:
  case usemtl:
  case texname: return 7;
  case comment: return 2;
  }
}

inline void mtl_parse_line( char* file, U64 idx, MTLLIB* mtl ) {
  LINE_TYPES first_two = get_first_two( file + idx );
  char* p = file + idx + line_type_offset( first_two );
  U64 mtl_idx = mtl->data.size - 1;

  switch( first_two ) {
  case newmtl: {
    char* end = p;
    while( *end && *end != '\n' ) ++end;
    U32 len = end - p;
    STR temp;
    temp.reserve( len + 1 );
    memcpy( temp.data, p, len );
    temp.data[len] = 0;
    mtl->data.resize( mtl->data.size + 1 );
    mtl->data[mtl->data.size - 1].name = temp.data;
  } break;
  case transparency: mtl->data[mtl_idx].transparency = strtof( p, &p ); break;
  case opticdensity: mtl->data[mtl_idx].opticdensity = strtof( p, &p ); break;
  case specular_exp: mtl->data[mtl_idx].specular_exp = strtof( p, &p ); break;
  case illum: mtl->data[mtl_idx].illum = (I32)strtol( p, &p, 10 ); break;
  case ambient_clr: mtl->data[mtl_idx].ambient_clr = read_vec3( p ); break;
  case diffuse_clr: mtl->data[mtl_idx].diffuse_clr = read_vec3( p ); break;
  case emissive_clr: mtl->data[mtl_idx].emissive_clr = read_vec3( p ); break;
  case specular_clr: mtl->data[mtl_idx].specular_clr = read_vec3( p ); break;
  case texname: {
    char* end = p;
    while( *end && *end != '\n' ) ++end;
    U32 len = end - p;
    STR temp;
    temp.reserve( len + 1 );
    memcpy( temp.data, p, len );
    temp.data[len] = 0;
    mtl->data[mtl_idx].map_Kd = temp.data;
  } break;
  case comment: break;
  default: {
    dlog(
      "mtl_from_file() : unhandled char combo %c%c\n",
      first_two & 0xFF, ( first_two >> 8 ) & 0xFF
    );
  } break;
  }
}

inline MTLLIB mtl_from_file( const char* path, const char* name ) {
  U64 mtl_size;
  STR full_path( path );
  full_path += name;
  char* f_mtl = (char*)file_read( full_path.data, &mtl_size );
  if( !f_mtl ) {
    dlog( "mtl_from_file() : failed to read file %s\n", full_path.data );
    return {};
  }

  MTLLIB mtl;
  mtl.name = name;
  for( U64 i = 0; i < mtl_size - 1 && f_mtl[i]; ++i ) {
    if( f_mtl[i] == '\n' )
      continue;
     mtl_parse_line( f_mtl, i, &mtl );
     while( f_mtl[i] && f_mtl[i] != '\n' ) ++i;
  }

   free( (void*)f_mtl );
  return mtl;
}

enum FACE_TYPES : U8 {
  V     = 0,
  VVT   = 1,
  VVN   = 2,
  VVTVN = 3
};

constexpr U32 STR_LEN = 0x40;

struct OBJFACE {
  FACE_TYPES type;
  LIST<U32> indices;
};

struct OBJDATA {
  LIST<VEC2> uvs;
  LIST<VEC3> normals;
  LIST<VEC3> vertices;
  LIST<MTLLIB> mtllib;
  struct UseMtl {
    LIST<U32> face_start;
    LIST<STR> mtl_name;
  } usemtl;
  LIST<OBJFACE> faces;
};

static inline U32 face_stride( FACE_TYPES t ) {
  switch( t ) {
  case V:     return 1;
  case VVT:   return 2;
  case VVN:   return 2;
  case VVTVN: return 3;
  }
  return 0;
}

inline OBJFACE triangulate_face( OBJFACE ref, U32 idx, U32 stride, FACE_TYPES type ) {
  OBJFACE tri{};
  tri.type = type;
  for( U32 ind = 0; ind < stride; ++ind )
    tri.indices.push( ref.indices.data[ind] );
  for( U32 ind = 0; ind < stride; ++ind )
    tri.indices.push( ref.indices.data[idx * stride + ind] );
  for( U32 ind = 0; ind < stride; ++ind )
    tri.indices.push( ref.indices.data[( idx + 1 ) * stride + ind] );
  return tri;
}

inline void obj_parse_f_line( OBJDATA* obj, char* p ) {
  FACE_TYPES type = (FACE_TYPES)0;
  bool type_set = 0;
  OBJFACE temp{};
  while( *p && *p != '\n' ) {
    bool has_vt = 0;
    bool has_vn = 0;
    U32 v = strtoul( p, &p, 10 );
    temp.indices.push( --v );
    if( *p == '/' ) {
      ++p;
      if( *p != '/' ) {
        U32 vt = strtoul( p, &p, 10 );
        temp.indices.push( --vt );
        has_vt = 1;
      }
      if( *p == '/' ) {
        ++p;
        U32 vn = strtoul( p, &p, 10 );
        temp.indices.push( --vn );
        has_vn = 1;
      }
    }
    if( !type_set ) {
      type_set = 1;
      if( has_vt && has_vn ) type = VVTVN;
      else if( has_vt )      type = VVT;
      else if( has_vn )      type = VVN;
      else                   type = V;
    }
    while( *p && *p != ' ' && *p != '\n' ) ++p;
    p = skip_space( p );
  }
  temp.type = type;
  U32 stride = face_stride( type );
  U32 count = temp.indices.size / stride;
  if ( count <= 3 )
    obj->faces.push( temp );
  else {
    for( U32 idx = 1; idx < count - 1; ++idx )
      obj->faces.push( triangulate_face( temp, idx, stride, type ) );
  }
}

inline void obj_parse_line( const char* path, char* file, U64 idx, OBJDATA* obj ) {
  LINE_TYPES first_two = get_first_two( file + idx );
  char* p = file + idx + line_type_offset( first_two );

  switch( first_two ) {
  case face: obj_parse_f_line( obj, p ); break;
  case vert:
  case vertnorm: {
    if( first_two == vert ) obj->vertices.push( read_vec3( p ) );
    else obj->normals.push( read_vec3( p ) );
  } break;
  case vertuv: {
    VEC2 temp{};
    p = skip_space( p );
    temp.x = strtof( p, &p );
    p = skip_space( p );
    temp.y = strtof( p, &p );
    obj->uvs.push( temp );
  } break;
  case mtlname:
  case usemtl: {
    bool is_mtllib = first_two == mtlname;
    char* end = p;
    while( *end && *end != '\n' ) ++end;
    U32 len = end - p;
    STR temp;
    temp.reserve( len + 1 );
    memcpy( temp.data, p, len );
    temp.data[len] = 0;
    if( !is_mtllib ) {
      if( obj->usemtl.mtl_name.size )
        obj->usemtl.face_start.push( obj->faces.size );
      obj->usemtl.mtl_name.push( temp.data );
    } else
      obj->mtllib.push( mtl_from_file( path, temp.data ) );
  } break;
  case comment:
  case group:
  case object:
  case shading: break;
  default: {
    dlog(
      "obj_from_file() : unhandled char combo %c%c\n",
      first_two & 0xFF, ( first_two >> 8 ) & 0xFF
    );
  } break;
  }
}

inline OBJDATA obj_from_file( const char* path, const char* name ) {
  char* f_obj = nullptr;
  STR obj_path( path );
  U64 obj_size = 0;
  obj_path += name;
  obj_path += ".obj";
  f_obj = (char*)file_read( obj_path.data, &obj_size );
  if( !f_obj ) {
    dlog( "obj_from_file() : failed to read obj file %s\n", obj_path.data );
    return {};
  }
  OBJDATA obj;
  for( U64 i = 0; i < obj_size - 1 && f_obj[i]; ++i ) {
    if( f_obj[i] == '\n' )
      continue;
    obj_parse_line( path, f_obj, i, &obj );
    while( f_obj[i] && f_obj[i] != '\n' ) ++i;
  }
  obj.usemtl.face_start.push( obj.faces.size );
  if ( !obj.usemtl.mtl_name.size )
    obj.usemtl.mtl_name.push( "default" );
  free( (void*)f_obj );
  return obj;
}

inline MODEL obj_to_model( OBJDATA obj ) {
  struct KEY { U32 v, vt, vn; };
  struct ENTRY { KEY key; U32 idx; };
  MODEL model;
  LIST<ENTRY> table;
  table.reserve( obj.faces.size * 3 );
  for( U32 f = 0; f < obj.faces.size; ++f ) {
    OBJFACE& face = obj.faces.data[f];
    U32 stride = face_stride( face.type );
    U32 count = face.indices.size / stride;
    for( U32 v = 0; v < count; ++v ) {
      KEY key = { 0xFFFFFFFF, 0xFFFFFFFF, 0xFFFFFFFF };
      U32 base = v * stride;
      key.v = face.indices.data[base + 0];
      if( face.type == VVT || face.type == VVTVN )
        key.vt = face.indices.data[base + 1];
      if( face.type == VVN )
        key.vn = face.indices.data[base + 1];
      if( face.type == VVTVN )
        key.vn = face.indices.data[base + 2];
      U32 idx = 0;
      bool found = 0;
      // this should be fine until ~1m triangles, right?
      for( U32 i = 0; i < table.size; ++i ) {
        ENTRY& e = table.data[i];
        if( e.key.vn != key.vn
         || e.key.vt != key.vt
         || e.key.v != key.v
        ) continue;
        idx = e.idx;
        found = 1;
        break;
      }
      if( !found ) {
        MODEL_DATA data{};
        data.position = obj.vertices.data[key.v];
        if( face.type == VVTVN || face.type == VVT )
          data.uv = obj.uvs.data[key.vt];
        if( face.type == VVTVN || face.type == VVN )
          data.normal = obj.normals.data[key.vn];
        idx = model.vertices.size;
        model.vertices.push( data );
        table.push({ key, idx });
      }
      model.indices.push( idx );
    }
  }
  return model;
}

inline MODEL model_from_obj( const I8* path, const I8* name ) {
  OBJDATA obj = obj_from_file( path, name );
  return obj_to_model( obj );
}