From a87abe8de58427f842b371c0169b4c41dc1cb6f1 Mon Sep 17 00:00:00 2001
From: Alexey <AlexeyAB@users.noreply.github.com>
Date: Sun, 23 Apr 2017 22:52:45 +0000
Subject: [PATCH] Update Readme.md
---
src/blas.c | 32 +++++++++++++++++++++++++++-----
1 files changed, 27 insertions(+), 5 deletions(-)
diff --git a/src/blas.c b/src/blas.c
index d6ab88b..31bd86b 100644
--- a/src/blas.c
+++ b/src/blas.c
@@ -5,8 +5,30 @@
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
+void reorg_cpu(float *x, int w, int h, int c, int batch, int stride, int forward, float *out)
+{
+ int b,i,j,k;
+ int out_c = c/(stride*stride);
-void reorg(float *x, int size, int layers, int batch, int forward)
+ for(b = 0; b < batch; ++b){
+ for(k = 0; k < c; ++k){
+ for(j = 0; j < h; ++j){
+ for(i = 0; i < w; ++i){
+ int in_index = i + w*(j + h*(k + c*b));
+ int c2 = k % out_c;
+ int offset = k / out_c;
+ int w2 = i*stride + offset % stride;
+ int h2 = j*stride + offset / stride;
+ int out_index = w2 + w*stride*(h2 + h*stride*(c2 + out_c*b));
+ if(forward) out[out_index] = x[in_index];
+ else out[in_index] = x[out_index];
+ }
+ }
+ }
+ }
+}
+
+void flatten(float *x, int size, int layers, int batch, int forward)
{
float *swap = calloc(size*layers*batch, sizeof(float));
int i,c,b;
@@ -189,12 +211,12 @@
if(input[i] > largest) largest = input[i];
}
for(i = 0; i < n; ++i){
- sum += exp(input[i]/temp-largest/temp);
+ float e = exp(input[i]/temp - largest/temp);
+ sum += e;
+ output[i] = e;
}
- if(sum) sum = largest/temp+log(sum);
- else sum = largest-100;
for(i = 0; i < n; ++i){
- output[i] = exp(input[i]/temp-sum);
+ output[i] /= sum;
}
}
--
Gitblit v1.10.0