QPMS
/
qpms
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

Test the translation operator array function 

Former-commit-id: b2ccef5adc99d27e24e730b817ce40ae0c98d4c2
This commit is contained in:
Marek Nečada 2018-05-16 13:01:14 +03:00
parent f1bb8a7c35
commit 12f9191a09
1 changed files with 9 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
qpms/tests/test_vswf_translations_array.c
View File

@ -109,13 +109,18 @@ int test_sphwave_translation(const qpms_trans_calculator *c, qpms_bessel_t wavet
if(QPMS_SUCCESS != qpms_vswf_fill(NULL, M1, N1, lMax, w1s, wavetype, c->normalisation)) if(QPMS_SUCCESS != qpms_vswf_fill(NULL, M1, N1, lMax, w1s, wavetype, c->normalisation))
abort(); // original wave set abort(); // original wave set
complex double A_whole[nelem][nelem], B_whole[nelem][nelem];
if (qpms_trans_calculator_get_AB_arrays(c,(complex double *) A_whole, (complex double *) B_whole,
1, nelem, ss, (w1s.r > ss.r), wavetype)) abort();
for(qpms_y_t y1 = 0; y1 < nelem; ++y1) { //index of the wave originating in o1 that will be reconstructed in o2 for(qpms_y_t y1 = 0; y1 < nelem; ++y1) { //index of the wave originating in o1 that will be reconstructed in o2
qpms_m_t m1; qpms_m_t m1;
qpms_l_t l1; qpms_l_t l1;
qpms_y2mn_p(y1, &m1, &l1); qpms_y2mn_p(y1, &m1, &l1);
printf("*** wave l = %d, m = %d ***\n", l1, m1); printf("*** wave l = %d, m = %d ***\n", l1, m1);
complex double A[nelem], B[nelem]; //complex double A[nelem], B[nelem];
/* /*
for(qpms_y_t y2 = 0; y2 < nelem; ++y2){ for(qpms_y_t y2 = 0; y2 < nelem; ++y2){
qpms_m_t m2; qpms_l_t l2; qpms_m_t m2; qpms_l_t l2;
@ -124,7 +129,9 @@ int test_sphwave_translation(const qpms_trans_calculator *c, qpms_bessel_t wavet
abort(); abort();
} }
*/ */
// array function instead HERE: // array function instead above
complex double *A = A_whole[y1]; // CHECKME right index ??
complex double *B = B_whole[y1];
printf("M = "); printf("M = ");