-
-#include <TH/TH.h>
-
/*
Example of FFI extension I started from:
*/
+#include <TH/TH.h>
+
#include "sequence_generator.h"
-int generate_sequence(long nb_sequences, THByteTensor *output) {
- long nb_images_per_sequence = 5;
- long depth = 3;
- long width = 64;
- long height = 64;
- long s;
+THByteTensor *generate_sequence(long nb_sequences, long nb_images_per_sequence, long image_width, long image_height) {
+ long nb_channels = 3;
unsigned char *a, *b;
- int c, k, i, j, st0, st1, st2, st3, st4;
+ long s, c, k, i, j, st0, st1, st2, st3, st4;
- THByteTensor_resize5d(output, nb_sequences, nb_images_per_sequence, depth, height, width);
+ THLongStorage *size = THLongStorage_newWithSize(5);
+ size->data[0] = nb_sequences;
+ size->data[1] = nb_images_per_sequence;
+ size->data[2] = nb_channels;
+ size->data[3] = image_height;
+ size->data[4] = image_width;
- st0 = THByteTensor_stride(output, 0);
- st1 = THByteTensor_stride(output, 1);
- st2 = THByteTensor_stride(output, 2);
- st3 = THByteTensor_stride(output, 3);
- st4 = THByteTensor_stride(output, 4);
+ THByteTensor *result = THByteTensor_newWithSize(size, NULL);
- a =
- THByteTensor_storage(output)->data + THByteTensor_storageOffset(output);
+ THLongStorage_free(size);
+
+ st0 = THByteTensor_stride(result, 0);
+ st1 = THByteTensor_stride(result, 1);
+ st2 = THByteTensor_stride(result, 2);
+ st3 = THByteTensor_stride(result, 3);
+ st4 = THByteTensor_stride(result, 4);
+
+ unsigned char tmp_buffer[nb_images_per_sequence * nb_channels * image_width * image_height];
for(s = 0; s < nb_sequences; s++) {
- unsigned char result[nb_images_per_sequence * depth * width * height];
- unsigned char *r = result;
- fl_generate_sequences(1, nb_images_per_sequence, width, height, result);
+ a =
+ THByteTensor_storage(result)->data + THByteTensor_storageOffset(result) + s * st0;
+
+ fl_generate_sequences(1, nb_images_per_sequence, image_width, image_height, tmp_buffer);
+
+ unsigned char *r = tmp_buffer;
for(k = 0; k < nb_images_per_sequence; k++) {
- for(c = 0; c < depth; c++) {
- for(i = 0; i < height; i++) {
+ for(c = 0; c < nb_channels; c++) {
+ for(i = 0; i < image_height; i++) {
b = a
- + s * st0 + k * st1 + c * st2 + i * st3;
- for(j = 0; j < width; j++) {
+ + k * st1 + c * st2 + i * st3;
+ for(j = 0; j < image_width; j++) {
*b = (unsigned char) (*r);
r++;
b += st4;
}
}
- return 1;
+ return result;
}