#include <stdio.h>
#include "Defs.h"
#include "math.h"
#include "neuralshooter.h"
#include <stdlib.h>

#define W_XS 10
#define W_YS 42

static REAL hW[W_XS][W_YS];
static REAL lW[W_XS];
static REAL bhW[W_XS];
static REAL blW = 0;

void readWeightFromFile(const char * filename){
	FILE * file;
	file = fopen(filename, "r");
	if(file != NULL){
		INT i,j;
		for(i=0; i<W_YS; i++){
		//	printf("linha : %d\n", i);
			for(j=0; j<W_XS; j++){
				fscanf(file, "%f\n", &hW[j][i]);
		//		printf("c[%d]= %f | ", j, hW[j][i]);
			}
		//	puts("\n");
		}
		for(i=0; i<W_XS; i++)
			fscanf(file, "%f\n", &bhW[i]);
		for(i=0; i<W_XS; i++)
			fscanf(file, "%f\n", &lW[i]);
		fscanf(file, "%f\n", &blW);
	}
	else
		fclose(file);
}

REAL tansig(REAL x){
	REAL y, up, down;
	REAL e = exp(x*-2);
	up = 2;
	down = 1 + e;
	y = (up/down)-1;
	return y;
}

REAL logsig(REAL x){
	REAL y, down;
	REAL e = exp(-x);
	y = 1 / ( 1 + e );
	return y;
}

REAL sim(REAL pattern[], INT size){
	REAL netO =0;
	INT i,j;
	REAL net, y;
	for(i=0; i<W_XS; i++){
		net = 0;
		for(j=0; j<W_YS; j++){
			net += hW[i][j] * pattern[j]; 
		}
		net += bhW[i];
		y = tansig(net);
		netO += y * lW[i];
	}
	netO += blW;
	return logsig(netO);
}

void readLog(REAL pattern[], const char * filename){
	FILE * f;
	REAL x;
	INT i,j;
	f = fopen(filename, "r");
	if(f != NULL){
		for(i=0; i<W_YS; i++){
			fscanf(f,"%f ", &x);
			pattern[i] = x/128.0;
			printf("%f\n", pattern[i]);
		}
	}
}
/*
int main(void){
	INT i,j;
	REAL pattern[42];
	readWeightFromFile("logW.txt");

//	for(j=0;j<42;j++)
//		pattern[j] = 0.0;
	readLog(pattern,"log.txt");
	printf("result %f \n", sim(pattern, 42));
	// teste das funcoes 
//	REAL x, y;
//	for(i=0; i< 10; i++){
//		x = ((2.0*rand())/RAND_MAX)-1;
//		y = tansig(x);
//		printf("tansig %f: %f \n",x,y); 
//	}
	
}*/