// Watermark embedding process described in Cox et al.
void ss_embed_to_signal(complex *freq_buffer, int f_buffer_len)
{ //{{{
// extract V, represented here by a set of indices from the
// frequency domain.
int *embed_indices;
// avoid sending dc offset at potential index
int noise_len = extract_sequence_indices(freq_buffer, f_buffer_len-1,
&embed_indices);
if(embed_iteration == embed_debug){
//printf("embed indices:\n");
//print_i_array(embed_indices, 20);
}
// noise_seq corresponds to W
double noise_seq[noise_len];
generate_noise(noise_seq, noise_len);
embed_to_noise(noise_seq, noise_len);
// embed W into V
for(int i = 0; i < noise_len; i++){
int freq_i = embed_indices[i];
watermark_elt(noise_seq[i], &(freq_buffer[freq_i]));
}
free(embed_indices);
} //}}}
void CpuPointSet::purturb_points (size_t group_size, float purturb_factor)
{
ASSERT_DIVIDES(group_size, size);
const size_t num_probes = size / group_size;
const size_t block_size = dim * group_size;
float sigma = purturb_factor * get_radius() / sqrt(dim);
LOG("purturbing point groups of size " << group_size << " by " << sigma);
Vector<int8_t> noise(block_size);
generate_noise(noise, sigma);
int8_t * restrict dx = noise;
for (size_t i_probe = 0; i_probe < num_probes; ++i_probe) {
uint8_t * restrict x = m_points + block_size * i_probe;
for (size_t i = 0; i < block_size; ++i) {
x[i] = bound_to(0, 255, int(x[i]) + int(dx[i]));
}
}
}
int main(int argc, char *argv[])
{
FLAC__uint32 test = 1;
unsigned channels;
int pattern01[] = { 1, -1, 0 };
int pattern02[] = { 1, 1, -1, 0 };
int pattern03[] = { 1, -1, -1, 0 };
int pattern04[] = { 1, -1, 1, -1, 0 };
int pattern05[] = { 1, -1, -1, 1, 0 };
int pattern06[] = { 1, -1, 1, 1, -1, 0 };
int pattern07[] = { 1, -1, -1, 1, -1, 0 };
(void)argc;
(void)argv;
is_big_endian_host = (*((FLAC__byte*)(&test)))? false : true;
#if !defined _MSC_VER && !defined __MINGW32__
struct timeval tv;
if(gettimeofday(&tv, 0) < 0) {
fprintf(stderr, "WARNING: couldn't seed RNG with time\n");
tv.tv_usec = 4321;
}
srandom(tv.tv_usec);
#else
srand(time(0));
#endif
if(!generate_01()) return 1;
if(!generate_02()) return 1;
if(!generate_03()) return 1;
if(!generate_04()) return 1;
if(!generate_fsd8("fsd8-01.raw", pattern01, 100)) return 1;
if(!generate_fsd8("fsd8-02.raw", pattern02, 100)) return 1;
if(!generate_fsd8("fsd8-03.raw", pattern03, 100)) return 1;
if(!generate_fsd8("fsd8-04.raw", pattern04, 100)) return 1;
if(!generate_fsd8("fsd8-05.raw", pattern05, 100)) return 1;
if(!generate_fsd8("fsd8-06.raw", pattern06, 100)) return 1;
if(!generate_fsd8("fsd8-07.raw", pattern07, 100)) return 1;
if(!generate_fsd16("fsd16-01.raw", pattern01, 100)) return 1;
if(!generate_fsd16("fsd16-02.raw", pattern02, 100)) return 1;
if(!generate_fsd16("fsd16-03.raw", pattern03, 100)) return 1;
if(!generate_fsd16("fsd16-04.raw", pattern04, 100)) return 1;
if(!generate_fsd16("fsd16-05.raw", pattern05, 100)) return 1;
if(!generate_fsd16("fsd16-06.raw", pattern06, 100)) return 1;
if(!generate_fsd16("fsd16-07.raw", pattern07, 100)) return 1;
if(!generate_fsd24("fsd24-01.raw", pattern01, 100)) return 1;
if(!generate_fsd24("fsd24-02.raw", pattern02, 100)) return 1;
if(!generate_fsd24("fsd24-03.raw", pattern03, 100)) return 1;
if(!generate_fsd24("fsd24-04.raw", pattern04, 100)) return 1;
if(!generate_fsd24("fsd24-05.raw", pattern05, 100)) return 1;
if(!generate_fsd24("fsd24-06.raw", pattern06, 100)) return 1;
if(!generate_fsd24("fsd24-07.raw", pattern07, 100)) return 1;
if(!generate_wbps16("wbps16-01.raw", 1000)) return 1;
if(!generate_sine8_1("sine8-00.raw", 48000.0, 200000, 441.0, 0.50, 441.0, 0.49)) return 1;
if(!generate_sine8_1("sine8-01.raw", 96000.0, 200000, 441.0, 0.61, 661.5, 0.37)) return 1;
if(!generate_sine8_1("sine8-02.raw", 44100.0, 200000, 441.0, 0.50, 882.0, 0.49)) return 1;
if(!generate_sine8_1("sine8-03.raw", 44100.0, 200000, 441.0, 0.50, 4410.0, 0.49)) return 1;
if(!generate_sine8_1("sine8-04.raw", 44100.0, 200000, 8820.0, 0.70, 4410.0, 0.29)) return 1;
if(!generate_sine8_2("sine8-10.raw", 48000.0, 200000, 441.0, 0.50, 441.0, 0.49, 1.0)) return 1;
if(!generate_sine8_2("sine8-11.raw", 48000.0, 200000, 441.0, 0.61, 661.5, 0.37, 1.0)) return 1;
if(!generate_sine8_2("sine8-12.raw", 96000.0, 200000, 441.0, 0.50, 882.0, 0.49, 1.0)) return 1;
if(!generate_sine8_2("sine8-13.raw", 44100.0, 200000, 441.0, 0.50, 4410.0, 0.49, 1.0)) return 1;
if(!generate_sine8_2("sine8-14.raw", 44100.0, 200000, 8820.0, 0.70, 4410.0, 0.29, 1.0)) return 1;
if(!generate_sine8_2("sine8-15.raw", 44100.0, 200000, 441.0, 0.50, 441.0, 0.49, 0.5)) return 1;
if(!generate_sine8_2("sine8-16.raw", 44100.0, 200000, 441.0, 0.61, 661.5, 0.37, 2.0)) return 1;
if(!generate_sine8_2("sine8-17.raw", 44100.0, 200000, 441.0, 0.50, 882.0, 0.49, 0.7)) return 1;
if(!generate_sine8_2("sine8-18.raw", 44100.0, 200000, 441.0, 0.50, 4410.0, 0.49, 1.3)) return 1;
if(!generate_sine8_2("sine8-19.raw", 44100.0, 200000, 8820.0, 0.70, 4410.0, 0.29, 0.1)) return 1;
if(!generate_sine16_1("sine16-00.raw", 48000.0, 200000, 441.0, 0.50, 441.0, 0.49)) return 1;
if(!generate_sine16_1("sine16-01.raw", 96000.0, 200000, 441.0, 0.61, 661.5, 0.37)) return 1;
if(!generate_sine16_1("sine16-02.raw", 44100.0, 200000, 441.0, 0.50, 882.0, 0.49)) return 1;
if(!generate_sine16_1("sine16-03.raw", 44100.0, 200000, 441.0, 0.50, 4410.0, 0.49)) return 1;
if(!generate_sine16_1("sine16-04.raw", 44100.0, 200000, 8820.0, 0.70, 4410.0, 0.29)) return 1;
if(!generate_sine16_2("sine16-10.raw", 48000.0, 200000, 441.0, 0.50, 441.0, 0.49, 1.0)) return 1;
if(!generate_sine16_2("sine16-11.raw", 48000.0, 200000, 441.0, 0.61, 661.5, 0.37, 1.0)) return 1;
if(!generate_sine16_2("sine16-12.raw", 96000.0, 200000, 441.0, 0.50, 882.0, 0.49, 1.0)) return 1;
if(!generate_sine16_2("sine16-13.raw", 44100.0, 200000, 441.0, 0.50, 4410.0, 0.49, 1.0)) return 1;
if(!generate_sine16_2("sine16-14.raw", 44100.0, 200000, 8820.0, 0.70, 4410.0, 0.29, 1.0)) return 1;
if(!generate_sine16_2("sine16-15.raw", 44100.0, 200000, 441.0, 0.50, 441.0, 0.49, 0.5)) return 1;
if(!generate_sine16_2("sine16-16.raw", 44100.0, 200000, 441.0, 0.61, 661.5, 0.37, 2.0)) return 1;
if(!generate_sine16_2("sine16-17.raw", 44100.0, 200000, 441.0, 0.50, 882.0, 0.49, 0.7)) return 1;
if(!generate_sine16_2("sine16-18.raw", 44100.0, 200000, 441.0, 0.50, 4410.0, 0.49, 1.3)) return 1;
if(!generate_sine16_2("sine16-19.raw", 44100.0, 200000, 8820.0, 0.70, 4410.0, 0.29, 0.1)) return 1;
if(!generate_sine24_1("sine24-00.raw", 48000.0, 200000, 441.0, 0.50, 441.0, 0.49)) return 1;
if(!generate_sine24_1("sine24-01.raw", 96000.0, 200000, 441.0, 0.61, 661.5, 0.37)) return 1;
if(!generate_sine24_1("sine24-02.raw", 44100.0, 200000, 441.0, 0.50, 882.0, 0.49)) return 1;
if(!generate_sine24_1("sine24-03.raw", 44100.0, 200000, 441.0, 0.50, 4410.0, 0.49)) return 1;
if(!generate_sine24_1("sine24-04.raw", 44100.0, 200000, 8820.0, 0.70, 4410.0, 0.29)) return 1;
if(!generate_sine24_2("sine24-10.raw", 48000.0, 200000, 441.0, 0.50, 441.0, 0.49, 1.0)) return 1;
if(!generate_sine24_2("sine24-11.raw", 48000.0, 200000, 441.0, 0.61, 661.5, 0.37, 1.0)) return 1;
if(!generate_sine24_2("sine24-12.raw", 96000.0, 200000, 441.0, 0.50, 882.0, 0.49, 1.0)) return 1;
if(!generate_sine24_2("sine24-13.raw", 44100.0, 200000, 441.0, 0.50, 4410.0, 0.49, 1.0)) return 1;
if(!generate_sine24_2("sine24-14.raw", 44100.0, 200000, 8820.0, 0.70, 4410.0, 0.29, 1.0)) return 1;
if(!generate_sine24_2("sine24-15.raw", 44100.0, 200000, 441.0, 0.50, 441.0, 0.49, 0.5)) return 1;
if(!generate_sine24_2("sine24-16.raw", 44100.0, 200000, 441.0, 0.61, 661.5, 0.37, 2.0)) return 1;
if(!generate_sine24_2("sine24-17.raw", 44100.0, 200000, 441.0, 0.50, 882.0, 0.49, 0.7)) return 1;
if(!generate_sine24_2("sine24-18.raw", 44100.0, 200000, 441.0, 0.50, 4410.0, 0.49, 1.3)) return 1;
if(!generate_sine24_2("sine24-19.raw", 44100.0, 200000, 8820.0, 0.70, 4410.0, 0.29, 0.1)) return 1;
/* WATCHOUT: the size of noise.raw is hardcoded into test/test_flac.sh */
if(!generate_noise("noise.raw", 65536 * 8 * 3)) return 1;
if(!generate_noise("noise8m32.raw", 32)) return 1;
if(!generate_wackywavs()) return 1;
for(channels = 1; channels <= 8; channels++) {
unsigned bits_per_sample;
for(bits_per_sample = 4; bits_per_sample <= 24; bits_per_sample++) {
static const unsigned nsamples[] = { 1, 111, 4777 } ;
unsigned samples;
for(samples = 0; samples < sizeof(nsamples)/sizeof(nsamples[0]); samples++) {
char fn[64];
sprintf(fn, "rt-%u-%u-%u.aiff", channels, bits_per_sample, nsamples[samples]);
if(!generate_aiff(fn, 44100, channels, bits_per_sample, nsamples[samples]))
return 1;
sprintf(fn, "rt-%u-%u-%u.wav", channels, bits_per_sample, nsamples[samples]);
if(!generate_wav(fn, 44100, channels, bits_per_sample, nsamples[samples], /*strict=*/true))
return 1;
if(bits_per_sample % 8 == 0) {
sprintf(fn, "rt-%u-%u-%u.raw", channels, bits_per_sample, nsamples[samples]);
if(!generate_noise(fn, channels * bits_per_sample/8 * nsamples[samples]))
return 1;
}
}
}
}
return 0;
}
int main(int argc, char *argv[])
{
FLAC__uint32 test = 1;
int pattern01[] = { 1, -1, 0 };
int pattern02[] = { 1, 1, -1, 0 };
int pattern03[] = { 1, -1, -1, 0 };
int pattern04[] = { 1, -1, 1, -1, 0 };
int pattern05[] = { 1, -1, -1, 1, 0 };
int pattern06[] = { 1, -1, 1, 1, -1, 0 };
int pattern07[] = { 1, -1, -1, 1, -1, 0 };
(void)argc;
(void)argv;
is_big_endian_host = (*((FLAC__byte*)(&test)))? false : true;
if(!generate_01()) return 1;
if(!generate_02()) return 1;
if(!generate_03()) return 1;
if(!generate_04()) return 1;
if(!generate_fsd8("fsd8-01.raw", pattern01, 100)) return 1;
if(!generate_fsd8("fsd8-02.raw", pattern02, 100)) return 1;
if(!generate_fsd8("fsd8-03.raw", pattern03, 100)) return 1;
if(!generate_fsd8("fsd8-04.raw", pattern04, 100)) return 1;
if(!generate_fsd8("fsd8-05.raw", pattern05, 100)) return 1;
if(!generate_fsd8("fsd8-06.raw", pattern06, 100)) return 1;
if(!generate_fsd8("fsd8-07.raw", pattern07, 100)) return 1;
if(!generate_fsd16("fsd16-01.raw", pattern01, 100)) return 1;
if(!generate_fsd16("fsd16-02.raw", pattern02, 100)) return 1;
if(!generate_fsd16("fsd16-03.raw", pattern03, 100)) return 1;
if(!generate_fsd16("fsd16-04.raw", pattern04, 100)) return 1;
if(!generate_fsd16("fsd16-05.raw", pattern05, 100)) return 1;
if(!generate_fsd16("fsd16-06.raw", pattern06, 100)) return 1;
if(!generate_fsd16("fsd16-07.raw", pattern07, 100)) return 1;
if(!generate_fsd24("fsd24-01.raw", pattern01, 100)) return 1;
if(!generate_fsd24("fsd24-02.raw", pattern02, 100)) return 1;
if(!generate_fsd24("fsd24-03.raw", pattern03, 100)) return 1;
if(!generate_fsd24("fsd24-04.raw", pattern04, 100)) return 1;
if(!generate_fsd24("fsd24-05.raw", pattern05, 100)) return 1;
if(!generate_fsd24("fsd24-06.raw", pattern06, 100)) return 1;
if(!generate_fsd24("fsd24-07.raw", pattern07, 100)) return 1;
if(!generate_wbps16("wbps16-01.raw", 1000)) return 1;
if(!generate_sine8_1("sine8-00.raw", 48000.0, 200000, 441.0, 0.50, 441.0, 0.49)) return 1;
if(!generate_sine8_1("sine8-01.raw", 96000.0, 200000, 441.0, 0.61, 661.5, 0.37)) return 1;
if(!generate_sine8_1("sine8-02.raw", 44100.0, 200000, 441.0, 0.50, 882.0, 0.49)) return 1;
if(!generate_sine8_1("sine8-03.raw", 44100.0, 200000, 441.0, 0.50, 4410.0, 0.49)) return 1;
if(!generate_sine8_1("sine8-04.raw", 44100.0, 200000, 8820.0, 0.70, 4410.0, 0.29)) return 1;
if(!generate_sine8_2("sine8-10.raw", 48000.0, 200000, 441.0, 0.50, 441.0, 0.49, 1.0)) return 1;
if(!generate_sine8_2("sine8-11.raw", 48000.0, 200000, 441.0, 0.61, 661.5, 0.37, 1.0)) return 1;
if(!generate_sine8_2("sine8-12.raw", 96000.0, 200000, 441.0, 0.50, 882.0, 0.49, 1.0)) return 1;
if(!generate_sine8_2("sine8-13.raw", 44100.0, 200000, 441.0, 0.50, 4410.0, 0.49, 1.0)) return 1;
if(!generate_sine8_2("sine8-14.raw", 44100.0, 200000, 8820.0, 0.70, 4410.0, 0.29, 1.0)) return 1;
if(!generate_sine8_2("sine8-15.raw", 44100.0, 200000, 441.0, 0.50, 441.0, 0.49, 0.5)) return 1;
if(!generate_sine8_2("sine8-16.raw", 44100.0, 200000, 441.0, 0.61, 661.5, 0.37, 2.0)) return 1;
if(!generate_sine8_2("sine8-17.raw", 44100.0, 200000, 441.0, 0.50, 882.0, 0.49, 0.7)) return 1;
if(!generate_sine8_2("sine8-18.raw", 44100.0, 200000, 441.0, 0.50, 4410.0, 0.49, 1.3)) return 1;
if(!generate_sine8_2("sine8-19.raw", 44100.0, 200000, 8820.0, 0.70, 4410.0, 0.29, 0.1)) return 1;
if(!generate_sine16_1("sine16-00.raw", 48000.0, 200000, 441.0, 0.50, 441.0, 0.49)) return 1;
if(!generate_sine16_1("sine16-01.raw", 96000.0, 200000, 441.0, 0.61, 661.5, 0.37)) return 1;
if(!generate_sine16_1("sine16-02.raw", 44100.0, 200000, 441.0, 0.50, 882.0, 0.49)) return 1;
if(!generate_sine16_1("sine16-03.raw", 44100.0, 200000, 441.0, 0.50, 4410.0, 0.49)) return 1;
if(!generate_sine16_1("sine16-04.raw", 44100.0, 200000, 8820.0, 0.70, 4410.0, 0.29)) return 1;
if(!generate_sine16_2("sine16-10.raw", 48000.0, 200000, 441.0, 0.50, 441.0, 0.49, 1.0)) return 1;
if(!generate_sine16_2("sine16-11.raw", 48000.0, 200000, 441.0, 0.61, 661.5, 0.37, 1.0)) return 1;
if(!generate_sine16_2("sine16-12.raw", 96000.0, 200000, 441.0, 0.50, 882.0, 0.49, 1.0)) return 1;
if(!generate_sine16_2("sine16-13.raw", 44100.0, 200000, 441.0, 0.50, 4410.0, 0.49, 1.0)) return 1;
if(!generate_sine16_2("sine16-14.raw", 44100.0, 200000, 8820.0, 0.70, 4410.0, 0.29, 1.0)) return 1;
if(!generate_sine16_2("sine16-15.raw", 44100.0, 200000, 441.0, 0.50, 441.0, 0.49, 0.5)) return 1;
if(!generate_sine16_2("sine16-16.raw", 44100.0, 200000, 441.0, 0.61, 661.5, 0.37, 2.0)) return 1;
if(!generate_sine16_2("sine16-17.raw", 44100.0, 200000, 441.0, 0.50, 882.0, 0.49, 0.7)) return 1;
if(!generate_sine16_2("sine16-18.raw", 44100.0, 200000, 441.0, 0.50, 4410.0, 0.49, 1.3)) return 1;
if(!generate_sine16_2("sine16-19.raw", 44100.0, 200000, 8820.0, 0.70, 4410.0, 0.29, 0.1)) return 1;
if(!generate_sine24_1("sine24-00.raw", 48000.0, 200000, 441.0, 0.50, 441.0, 0.49)) return 1;
if(!generate_sine24_1("sine24-01.raw", 96000.0, 200000, 441.0, 0.61, 661.5, 0.37)) return 1;
if(!generate_sine24_1("sine24-02.raw", 44100.0, 200000, 441.0, 0.50, 882.0, 0.49)) return 1;
if(!generate_sine24_1("sine24-03.raw", 44100.0, 200000, 441.0, 0.50, 4410.0, 0.49)) return 1;
if(!generate_sine24_1("sine24-04.raw", 44100.0, 200000, 8820.0, 0.70, 4410.0, 0.29)) return 1;
if(!generate_sine24_2("sine24-10.raw", 48000.0, 200000, 441.0, 0.50, 441.0, 0.49, 1.0)) return 1;
if(!generate_sine24_2("sine24-11.raw", 48000.0, 200000, 441.0, 0.61, 661.5, 0.37, 1.0)) return 1;
if(!generate_sine24_2("sine24-12.raw", 96000.0, 200000, 441.0, 0.50, 882.0, 0.49, 1.0)) return 1;
if(!generate_sine24_2("sine24-13.raw", 44100.0, 200000, 441.0, 0.50, 4410.0, 0.49, 1.0)) return 1;
if(!generate_sine24_2("sine24-14.raw", 44100.0, 200000, 8820.0, 0.70, 4410.0, 0.29, 1.0)) return 1;
if(!generate_sine24_2("sine24-15.raw", 44100.0, 200000, 441.0, 0.50, 441.0, 0.49, 0.5)) return 1;
if(!generate_sine24_2("sine24-16.raw", 44100.0, 200000, 441.0, 0.61, 661.5, 0.37, 2.0)) return 1;
if(!generate_sine24_2("sine24-17.raw", 44100.0, 200000, 441.0, 0.50, 882.0, 0.49, 0.7)) return 1;
if(!generate_sine24_2("sine24-18.raw", 44100.0, 200000, 441.0, 0.50, 4410.0, 0.49, 1.3)) return 1;
if(!generate_sine24_2("sine24-19.raw", 44100.0, 200000, 8820.0, 0.70, 4410.0, 0.29, 0.1)) return 1;
if(!generate_noise("noise.raw", 65536 * 8 * 3)) return 1;
if(!generate_wackywavs()) return 1;
return 0;
}
int
main (int argc, char *argv[])
{
SchroEncoder *encoder;
SchroParams params;
SchroVideoFormat video_format;
int filter;
SchroFrame *frame;
int transform_depth;
int i,j;
int k;
schro_init();
encoder = schro_encoder_new ();
frame = schro_frame_new_and_alloc (NULL, SCHRO_FRAME_FORMAT_S16_444, SIZE, SIZE);
schro_fft_generate_tables_f32 (costable, sintable, 2*SHIFT);
filter = 6;
transform_depth = 4;
schro_video_format_set_std_video_format (&video_format, 0);
video_format.width = SIZE;
video_format.height = SIZE;
memset (¶ms, 0, sizeof(params));
params.video_format = &video_format;
schro_params_init (¶ms, 0);
params.wavelet_filter_index = filter;
params.transform_depth = transform_depth;
schro_params_calculate_iwt_sizes (¶ms);
for(i=0;i<SIZE*SIZE;i++) power[i] = 0;
for(k=0;k<N_TRIALS;k++){
generate_noise (frame, transform_depth,
encoder->subband_weights[filter][transform_depth-1]);
schro_frame_inverse_iwt_transform (frame, ¶ms, tmp);
for(j=0;j<SIZE;j++){
int16_t *line;
line = OFFSET(frame->components[0].data,
frame->components[0].stride * j);
for(i=0;i<SIZE;i++){
sr[j*SIZE+i] = line[i];
si[j*SIZE+i] = 0;
}
}
schro_fft_fwd_f32 (dr, di, sr, si, costable, sintable, 2*SHIFT);
for(i=0;i<SIZE*SIZE;i++) {
power[i] += (dr[i]*dr[i]+di[i]*di[i])*(1.0/(SIZE*SIZE));
}
}
//power[0] *= 4.0/SIZE;
//power[0] = 0;
for(j=0;j<SIZE/2;j+=CHUNK_SIZE){
for(i=0;i<SIZE/2;i+=CHUNK_SIZE){
int ii,jj;
double sum = 0;
for(jj=0;jj<CHUNK_SIZE;jj++){
for(ii=0;ii<CHUNK_SIZE;ii++){
sum += power[(j+jj)*SIZE+(i+ii)];
}
}
sum /= N_TRIALS*CHUNK_SIZE*CHUNK_SIZE;
printf("%d %d %g\n", j, i, sqrt(sum)/AMPLITUDE);
}
printf("\n");
}
return fail;
}
int main() {
int k = 1, err = 10, iter = 0;
float ratio = 0;
int A[4][5] = {
{ 1, 1, 1, 1, 1},
{ 1, 0, 0, 0, 1},
{ 1, 1, 1, 1, 1},
{ 1, 0, 0, 0, 1} };
int C[4][5] = {
{ 1, 1, 1, 1, 1},
{ 1, 0, 0, 0, 0},
{ 1, 0, 0, 0, 0},
{ 1, 1, 1, 1, 1} };
int L[4][5];
int last_used_matrix[4][5] = {
{ -1, -1, -1, -1, -1},
{ -1, -1, -1, -1, -1},
{ -1, -1, -1, -1, -1},
{ -1, -1, -1, -1, -1} };
t_neuron d_neuron = { .potential = 0, .threshold = 0.5, .output = 0 };
t_synapse synapses[20];
// Création du réseau
for (int i = 0; i < 20; i++) {
synapses[i].source = create_sensor();
synapses[i].target = &d_neuron;
synapses[i].weight = 0;
}
// Phase d'apprentissage
do {
feed_neural_network(&synapses, A);
learning_phase(&synapses, EPSILON, 1, 0);
feed_neural_network(&synapses, C);
learning_phase(&synapses, EPSILON, 0, 1);
feed_neural_network(&synapses, A);
calculate_potential(&synapses, &d_neuron);
if (d_neuron.output != 1)
err++;
feed_neural_network(&synapses, C);
calculate_potential(&synapses, &d_neuron);
if (d_neuron.output != 0)
err++;
iter++;
ratio = (float) iter / (float) err;
//printf("iter: %d && err: %d\n", iter, err);
//printf("ratio:%f\n", ratio);
} while ( ratio < 0.8);
// Phase de test pour A
while (k < 20) {
int i = 0;
int success = 0;
while (i < 100000) {
while (true) {
// Generer copie d'une lettre
lcpy(&L, &A);
// Generer du bruit sur cette copie
generate_noise(&L, k);
//print_letter(L);
//printf("\n");
if (!matrices_are_equal(L, last_used_matrix)) {
lcpy(&last_used_matrix, &L);
break;
}
}
feed_neural_network(&synapses, L);
calculate_potential(&synapses, &d_neuron);
if (d_neuron.output == 1)
success++;
i++;
}
success /= 1000;
printf("%d%% de succes sur %d essais avec %d%% de bruit.\n", success, i, k*5);
k++;
}
// Phase de test pour C
k = 1;
while (k < 20) {
int i = 0;
int success = 0;
while (i < 100000) {
while (true) {
// Generer copie d'une lettre
lcpy(&L, &C);
// Generer du bruit sur cette copie
generate_noise(&L, k);
if (!matrices_are_equal(L, last_used_matrix)) {
lcpy(&last_used_matrix, &L);
break;
}
}
feed_neural_network(&synapses, L);
calculate_potential(&synapses, &d_neuron);
if (d_neuron.output == 0)
success++;
i++;
}
success /= 1000;
printf("%d%% de succes sur %d essais avec %d%% de bruit.\n", success, i, k*5);
k++;
}
}
