c - FFTW real-to-real transform strided array -

-

i have matrix stored in row-major order. trying compute dct of ub-matrix using fftw, , nonsense. in following paragraphs describe problem , solution, , in understanding why not work.

given i , l, compute dct of sub-matrix consists of rows k k mod l == i. example, assume l = 3 , i = 2. in following matrix, sub-matrix wish transform marked in red (2 mod 3 = 2, 5 mod 3 = 2, 8 mod 3 = 2).

matrix

the source , destination arrays have same layout, , transformed matrix should stored in same location in destination array.

void transform(double* src, double* dest, size_t rows, size_t cols, size_t l, size_t i) {     int rank = 2;     fftw_iodim64 dims[] = {         { rows / l, l, l },         { cols, rows, rows } };     fftw_r2r_kind kind = fftw_redft10;      fftw_plan plan = fftw_plan_guru64_r2r(rank, dims, 0, null, src + l, dest + l, &kind, fftw_estimate | fftw_unaligned | fftw_preserve_input);     fftw_execute(plan);     fftw_destroy_plan(pla); }

update

i tested simple case when i=l=1. in case nonsense. tested 3x4 matrix 1 of dct basis vectors:

a(i,j) = cos((i + 0.5)*2*pi/3) * cos((j + 0.5)*3*pi/4)

i expected result elements (close to) 0 except one. resulting matrix looks this:

0             -2.22045e-016  1.33227e-015  2.22045e-016 2.22045e-016  -2.77556e-016  9.99201e-016  5.55112e-017 -8.88178e-016 -1.62359       7.83938       1.62359

seems quite weird.

update 2

i tested simple matrix (0,0) element 1 , rest zero. in case, too, i=l=1 (the sub-matrix whole matrix). here result:

      2       2       2       0 1.73205 1.73205 1.73205       0       1       1       1       0

kind supposed array of size rank. see e.g. http://www.fftw.org/doc/guru-real_002dto_002dreal-transforms.html#guru-real_002dto_002dreal-transforms.


Comments

Post a Comment

Popular posts from this blog

Qt4: how to send QString inside a struct via QSharedMemory -

-
i have struct struct control_data{ int column_number; qstring cell; }; i need send thread of qsharememory. read can't because qstring contains pointers inside. other ways? you have serialize struct byte array. can convert qstring const char* this: mystring.tostdstring().c_str(); but serializing qstring should work. the first step serialize struct qdatastream using qt, example here . then once struct can read , written can pass shared memory. a complete example of using qsharedmemory can found here . here relevant code: // first, test whether shared memory segment attached process. // if so, detach if (sharedmem.isattached()) { sharedmem.detach(); } ... qbuffer buffer; buffer.open( qbuffer::readwrite ); qdatastream out( &buffer ); out << youstruct; int size = buffer.size(); // size of int + size of qstring in bytes if ( !sharedmem.create( size ) ) { return; } // write shared memory sharedmem.lock(); char *to = (char*...
Read more

python - jinja2: TemplateSyntaxError: expected token ',', got 'string' -

-
i new jinja2 , having issue using python regular expression (re). in following code bold lines have error string in them. {% block content %} <div class="container"> {% l in lines %} {% if re.search(r"error", l) %} {# <<< throws error #} <b> {{ l }} </b> {% else %} {{ l }} <hr> {% endif %} {% endfor %} </div> {% endblock %} the re.search above throws following error: jinja2.exceptions.templatesyntaxerror templatesyntaxerror: expected token ',', got 'string' raw python code not supported in jinja2 template syntax. {% if re.search(r"error", l) %} replace line with {% if "error" in l %} can fix problem. if logical condition more complicated, should consider defining own custom filters (which can call python code) or complicated things in view layer. go check global namespace .
Read more

node.js - NodeJS remote terminal to Dropbear OpenWRT-Server -

-
i wandering whether there valid solution have remote terminal openwrt-box out of nodejs-app? connecting terminal works: ssh -i ~/.mykeys/id_rsa root@192.168.178.39 busybox v1.23.2 (2015-04-22 23:25:48 utc) built-in shell (ash) root@openwrt:~# the only interactive ssh solution nodejs doesn't interactive part described in readme.md following: var client = require('ssh2').client; var conn = new client(); conn.on('ready', function() { console.log('client :: ready'); conn.shell(function(err, stream) { if (err) throw err; stream.on('close', function() { console.log('stream :: close'); conn.end(); }).on('data', function(data) { console.log('stdout: ' + data); }).stderr.on('data', function(data) { console.log('stderr: ' + data); }); stream.end('ls -l\nexit\n'); }); }).connect({ host: '192.168.100.100', port: 22, username...
Read more