diff --git a/bitpackv.c b/bitpackv.c
new file mode 100644
index 0000000..34a9e91
--- /dev/null
+++ b/bitpackv.c
@@ -0,0 +1,105 @@
+/**
+    Copyright (C) powturbo 2013-2015
+    GPL v2 License
+  
+    This program is free software; you can redistribute it and/or modify
+    it under the terms of the GNU General Public License as published by
+    the Free Software Foundation; either version 2 of the License, or
+    (at your option) any later version.
+
+    This program is distributed in the hope that it will be useful,
+    but WITHOUT ANY WARRANTY; without even the implied warranty of
+    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+    GNU General Public License for more details.
+
+    You should have received a copy of the GNU General Public License along
+    with this program; if not, write to the Free Software Foundation, Inc.,
+    51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
+
+    - homepage : https://sites.google.com/site/powturbo/
+    - github   : https://github.com/powturbo
+    - twitter  : https://twitter.com/powturbo
+    - email    : powturbo [_AT_] gmail [_DOT_] com
+**/
+//    bitpackv.c - "Integer Compression" SIMD bit packing 
+  #ifndef VSTI
+#include <emmintrin.h>
+#include "bitpack.h"
+#include "bitutil.h"
+
+#define PAD8(__x) (((__x)+8-1)/8)
+
+#define VSTI(ip, i, iv, parm)
+#define IPP(ip, i, iv) _mm_loadu_si128(ip++)
+#include __FILE__ 
+  
+unsigned char *bitpackv32(unsigned       *__restrict in, unsigned n, unsigned char *__restrict out, unsigned b) { unsigned char *pout = out+PAD8(n*b); BITPACKV32(in, n, b, out, 0); return pout; }
+unsigned char *bitpackv16(unsigned short *__restrict in, unsigned n, unsigned char *__restrict out, unsigned b) { unsigned char *pout = out+PAD8(n*b); BITPACKV32(in, n, b, out, 0); return pout; }
+#undef VSTI 
+#undef IPP
+
+//------------------------------------------------------------------------------------------------------------------------------
+#define VSTI(__ip, __i, __iv, __sv) v = _mm_loadu_si128(__ip++); __iv = DELTAV(v,__sv); __sv = v
+#define IPP(ip, i, __iv) __iv
+#include __FILE__ 
+
+unsigned char *bitdpackv32(unsigned       *__restrict in, unsigned n, unsigned char *__restrict out, unsigned start, unsigned b) { unsigned char *pout = out+PAD8(n*b);
+  __m128i v,sv = _mm_set1_epi32(start);
+  BITPACKV32(in, n, b, out, sv); return pout; 
+}
+#undef VSTI
+
+//------------------------------------------------------------------------------------------------------------------------------
+#define VSTI(__ip, __i, __iv, __sv) v = _mm_loadu_si128(__ip++); __iv = _mm_sub_epi32(DELTAV(v,__sv),cv); __sv = v
+
+unsigned char *bitd1packv32(unsigned       *__restrict in, unsigned n, unsigned char *__restrict out, unsigned start, unsigned b) { unsigned char *pout = out+PAD8(n*b);
+  __m128i v, sv = _mm_set1_epi32(start), cv = _mm_set1_epi32(1);
+  BITPACKV32(in, n, b, out, sv); return pout; 
+}
+  #else
+#include <strings.h>
+#include <emmintrin.h>
+
+#define OPPE(__op)
+#define IPPE(__op)
+
+#include "bitpackv32_.h"
+ 
+#define BITPACKV32(__pip, __n, __nbits, __pop, __parm) { __m128i *__ip=(__m128i *)__pip,*__op=(__m128i *)__pop; typeof(__ip[0]) *__ipe=(__ip)+(__n);\
+  switch(__nbits) {\
+    case  0:  												    		    break;\
+    case  1:{  BITPACKV32_1( __ip, __op, __parm); /*while(__ip < __ipe);*/} break;\
+    case  2:{  BITPACKV32_2( __ip, __op, __parm); /*while(__ip < __ipe)*/;} break;\
+    case  3:{  BITPACKV32_3( __ip, __op, __parm); /*while(__ip < __ipe);*/} break;\
+    case  4:{  BITPACKV32_4( __ip, __op, __parm); /*while(__ip < __ipe);*/} break;\
+    case  5:{  BITPACKV32_5( __ip, __op, __parm); /*while(__ip < __ipe);*/} break;\
+    case  6:{  BITPACKV32_6( __ip, __op, __parm); /*while(__ip < __ipe);*/} break;\
+    case  7:{  BITPACKV32_7( __ip, __op, __parm); /*while(__ip < __ipe);*/} break;\
+    case  8:{  BITPACKV32_8( __ip, __op, __parm); /*while(__ip < __ipe);*/} break;\
+    case  9:{  BITPACKV32_9( __ip, __op, __parm); /*while(__ip < __ipe);*/} break;\
+    case 10:{  BITPACKV32_10(__ip, __op, __parm); /*while(__ip < __ipe);*/} break;\
+    case 11:{  BITPACKV32_11(__ip, __op, __parm); /*while(__ip < __ipe);*/} break;\
+    case 12:{  BITPACKV32_12(__ip, __op, __parm); /*while(__ip < __ipe);*/} break;\
+    case 13:{  BITPACKV32_13(__ip, __op, __parm); /*while(__ip < __ipe);*/} break;\
+    case 14:{  BITPACKV32_14(__ip, __op, __parm); /*while(__ip < __ipe);*/} break;\
+    case 15:{  BITPACKV32_15(__ip, __op, __parm); /*while(__ip < __ipe);*/} break;\
+    case 16:{  BITPACKV32_16(__ip, __op, __parm); /*while(__ip < __ipe);*/} break;\
+    case 17:{  BITPACKV32_17(__ip, __op, __parm); /*while(__ip < __ipe);*/} break;\
+    case 18:{  BITPACKV32_18(__ip, __op, __parm); /*while(__ip < __ipe);*/} break;\
+    case 19:{  BITPACKV32_19(__ip, __op, __parm); /*while(__ip < __ipe);*/} break;\
+    case 20:{  BITPACKV32_20(__ip, __op, __parm); /*while(__ip < __ipe);*/} break;\
+    case 21:{  BITPACKV32_21(__ip, __op, __parm); /*while(__ip < __ipe);*/} break;\
+    case 22:{  BITPACKV32_22(__ip, __op, __parm); /*while(__ip < __ipe);*/} break;\
+    case 23:{  BITPACKV32_23(__ip, __op, __parm); /*while(__ip < __ipe);*/} break;\
+    case 24:{  BITPACKV32_24(__ip, __op, __parm); /*while(__ip < __ipe);*/} break;\
+    case 25:{  BITPACKV32_25(__ip, __op, __parm); /*while(__ip < __ipe);*/} break;\
+    case 26:{  BITPACKV32_26(__ip, __op, __parm); /*while(__ip < __ipe);*/} break;\
+    case 27:{  BITPACKV32_27(__ip, __op, __parm); /*while(__ip < __ipe);*/} break;\
+    case 28:{  BITPACKV32_28(__ip, __op, __parm); /*while(__ip < __ipe);*/} break;\
+    case 29:{  BITPACKV32_29(__ip, __op, __parm); /*while(__ip < __ipe);*/} break;\
+    case 30:{  BITPACKV32_30(__ip, __op, __parm); /*while(__ip < __ipe);*/} break;\
+    case 31:{  BITPACKV32_31(__ip, __op, __parm); /*while(__ip < __ipe);*/} break;\
+    case 32:{  BITPACKV32_32(__ip, __op, __parm); /*while(__ip < __ipe);*/} break;\
+  }\
+}
+  #endif