コード例 #1
0
int main()
{
RBM *a = new RBM;

//a->input(stdin);
// for xor
a->build(2,2);  
int x[2],y[2],s[2],sz[2];
x[0] = -1;
x[1] =  1;
y[0] =  1;
y[1] = -1;
s[0] =  1;
s[1] =  1;
sz[0] = -1;
sz[1] = -1;
for( int i=0; i< 100; i++)
	{
a->train((float)1.001,2,x);
a->train((float)1.001,2,y);
printf( "xor %f %f same %f %f\n",a->energy(2,x),a->energy(2,y),a->energy(2,s),a->energy(2,sz));

a->dump(stdout);
	}

s[1] = 0;
a->fillin(2,s);
printf("%d %d\n",s[0],s[1]);
}
コード例 #2
0
int main()
{
RBM *a = new RBM;
const int code_width = 5;
const int window_width = 40;

FILE *ip; 
ip = fopen("one.10.rbm","r");
a->input(ip);
fclose(ip);
ip = fopen("pdbaanr","r");

fasta *d = new fasta(ip);

// static issue WTF man
static int encoded[window_width*code_width];
static int test[window_width*code_width];



while(d->get_next())
  while( d->encode(window_width,window_width*code_width,encoded))
  {
	int same;

	for( int ia=0;ia<window_width; ia++)
	{
	same = 0;
	for( int i=0; i< window_width*code_width; i++)
		test[i] = encoded[i];
	for( int i=ia*code_width; i< ia*code_width+code_width; i++)
		test[i] = 0;
//	for( int i=ia*code_width; i< ia*code_width+code_width; i++)
//		printf("%d %d\n", encoded[i],test[i]);

	a->fillin(window_width*code_width, test);
//	for( int i=ia*code_width; i< ia*code_width+code_width; i++)
//		printf("%d %d\n", encoded[i],test[i]);
	for( int i=ia*code_width; i< ia*code_width+code_width; i++)
		if( test[i] == encoded[i]) same ++;
	printf("%d ", same);
   
	
	}
	printf("\n");
//     a->train((float)1.001,window_width*code_width,encoded);
  }
}