|
|
@ -114,16 +114,6 @@ int main(void){ |
|
|
|
} |
|
|
|
} |
|
|
|
printf("};\n"); |
|
|
|
printf("};\n"); |
|
|
|
{ |
|
|
|
{ |
|
|
|
int level=2; |
|
|
|
|
|
|
|
int w= width >> (s.spatial_decomposition_count-level); |
|
|
|
|
|
|
|
//int h= height >> (s.spatial_decomposition_count-level);
|
|
|
|
|
|
|
|
int stride= width << (s.spatial_decomposition_count-level); |
|
|
|
|
|
|
|
DWTELEM *buf= buffer[0]; |
|
|
|
|
|
|
|
int64_t error=0; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
buf+=w; |
|
|
|
|
|
|
|
buf+=stride>>1; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
memset(buffer[0], 0, sizeof(int)*width*height); |
|
|
|
memset(buffer[0], 0, sizeof(int)*width*height); |
|
|
|
for(y=0; y<height; y++){ |
|
|
|
for(y=0; y<height; y++){ |
|
|
|
for(x=0; x<width; x++){ |
|
|
|
for(x=0; x<width; x++){ |
|
|
@ -135,7 +125,6 @@ int main(void){ |
|
|
|
for(y=0; y<height; y++){ |
|
|
|
for(y=0; y<height; y++){ |
|
|
|
for(x=0; x<width; x++){ |
|
|
|
for(x=0; x<width; x++){ |
|
|
|
int64_t d= buffer[0][x + y*width]; |
|
|
|
int64_t d= buffer[0][x + y*width]; |
|
|
|
error += d*d; |
|
|
|
|
|
|
|
if(FFABS(width/2-x)<9 && FFABS(height/2-y)<9) printf("%8"PRId64" ", d); |
|
|
|
if(FFABS(width/2-x)<9 && FFABS(height/2-y)<9) printf("%8"PRId64" ", d); |
|
|
|
} |
|
|
|
} |
|
|
|
if(FFABS(height/2-y)<9) printf("\n"); |
|
|
|
if(FFABS(height/2-y)<9) printf("\n"); |
|
|
|