From e7688a05a194e3c8baf3c11fbf09b7f5e8640a77 Mon Sep 17 00:00:00 2001
From: Joseph Redmon <pjreddie@gmail.com>
Date: Wed, 06 May 2015 21:08:16 +0000
Subject: [PATCH] no idea

---
 src/cost_layer.c |   82 +++++++++++++++++++++++++++++++----------
 1 files changed, 62 insertions(+), 20 deletions(-)

diff --git a/src/cost_layer.c b/src/cost_layer.c
index 1df0ed4..1f36232 100644
--- a/src/cost_layer.c
+++ b/src/cost_layer.c
@@ -1,53 +1,95 @@
 #include "cost_layer.h"
 #include "utils.h"
-#include "mini_blas.h"
+#include "cuda.h"
+#include "blas.h"
 #include <math.h>
+#include <string.h>
 #include <stdlib.h>
 #include <stdio.h>
 
-cost_layer *make_cost_layer(int batch, int inputs)
+COST_TYPE get_cost_type(char *s)
+{
+    if (strcmp(s, "sse")==0) return SSE;
+    if (strcmp(s, "masked")==0) return MASKED;
+    fprintf(stderr, "Couldn't find activation function %s, going with SSE\n", s);
+    return SSE;
+}
+
+char *get_cost_string(COST_TYPE a)
+{
+    switch(a){
+        case SSE:
+            return "sse";
+        case MASKED:
+            return "masked";
+    }
+    return "sse";
+}
+
+cost_layer *make_cost_layer(int batch, int inputs, COST_TYPE type)
 {
     fprintf(stderr, "Cost Layer: %d inputs\n", inputs);
     cost_layer *layer = calloc(1, sizeof(cost_layer));
     layer->batch = batch;
     layer->inputs = inputs;
+    layer->type = type;
     layer->delta = calloc(inputs*batch, sizeof(float));
     layer->output = calloc(1, sizeof(float));
     #ifdef GPU
-    layer->delta_cl = cl_make_array(layer->delta, inputs*batch);
+    layer->delta_gpu = cuda_make_array(layer->delta, inputs*batch);
     #endif
     return layer;
 }
 
-void forward_cost_layer(cost_layer layer, float *input, float *truth)
+void forward_cost_layer(cost_layer layer, network_state state)
 {
-    if (!truth) return;
-    copy_cpu(layer.batch*layer.inputs, truth, 1, layer.delta, 1);
-    axpy_cpu(layer.batch*layer.inputs, -1, input, 1, layer.delta, 1);
+    if (!state.truth) return;
+    if(layer.type == MASKED){
+        int i;
+        for(i = 0; i < layer.batch*layer.inputs; ++i){
+            if(state.truth[i] == 0) state.input[i] = 0;
+        }
+    }
+    copy_cpu(layer.batch*layer.inputs, state.truth, 1, layer.delta, 1);
+    axpy_cpu(layer.batch*layer.inputs, -1, state.input, 1, layer.delta, 1);
     *(layer.output) = dot_cpu(layer.batch*layer.inputs, layer.delta, 1, layer.delta, 1);
+    //printf("cost: %f\n", *layer.output);
 }
 
-void backward_cost_layer(const cost_layer layer, float *input, float *delta)
+void backward_cost_layer(const cost_layer layer, network_state state)
 {
-    copy_cpu(layer.batch*layer.inputs, layer.delta, 1, delta, 1);
+    copy_cpu(layer.batch*layer.inputs, layer.delta, 1, state.delta, 1);
 }
 
 #ifdef GPU
-void forward_cost_layer_gpu(cost_layer layer, cl_mem input, cl_mem truth)
+
+void pull_cost_layer(cost_layer layer)
 {
-    if (!truth) return;
-
-    copy_ongpu(layer.batch*layer.inputs, truth, 1, layer.delta_cl, 1);
-    axpy_ongpu(layer.batch*layer.inputs, -1, input, 1, layer.delta_cl, 1);
-
-    cl_read_array(layer.delta_cl, layer.delta, layer.batch*layer.inputs);
-    *(layer.output) = dot_cpu(layer.batch*layer.inputs, layer.delta, 1, layer.delta, 1);
-    //printf("%f\n", *layer.output);
+    cuda_pull_array(layer.delta_gpu, layer.delta, layer.batch*layer.inputs);
 }
 
-void backward_cost_layer_gpu(const cost_layer layer, cl_mem input, cl_mem delta)
+void push_cost_layer(cost_layer layer)
 {
-    copy_ongpu(layer.batch*layer.inputs, layer.delta_cl, 1, delta, 1);
+    cuda_push_array(layer.delta_gpu, layer.delta, layer.batch*layer.inputs);
+}
+
+void forward_cost_layer_gpu(cost_layer layer, network_state state)
+{
+    if (!state.truth) return;
+    if (layer.type == MASKED) {
+        mask_ongpu(layer.batch*layer.inputs, state.input, state.truth);
+    }
+    
+    copy_ongpu(layer.batch*layer.inputs, state.truth, 1, layer.delta_gpu, 1);
+    axpy_ongpu(layer.batch*layer.inputs, -1, state.input, 1, layer.delta_gpu, 1);
+
+    cuda_pull_array(layer.delta_gpu, layer.delta, layer.batch*layer.inputs);
+    *(layer.output) = dot_cpu(layer.batch*layer.inputs, layer.delta, 1, layer.delta, 1);
+}
+
+void backward_cost_layer_gpu(const cost_layer layer, network_state state)
+{
+    copy_ongpu(layer.batch*layer.inputs, layer.delta_gpu, 1, state.delta, 1);
 }
 #endif
 

--
Gitblit v1.10.0