From 12f9191a09c1d6d87cc1981b160d8b2b3373d4fe Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Marek=20Ne=C4=8Dada?= Date: Wed, 16 May 2018 13:01:14 +0300 Subject: [PATCH] Test the translation operator array function Former-commit-id: b2ccef5adc99d27e24e730b817ce40ae0c98d4c2 --- qpms/tests/test_vswf_translations_array.c | 11 +++++++++-- 1 file changed, 9 insertions(+), 2 deletions(-) diff --git a/qpms/tests/test_vswf_translations_array.c b/qpms/tests/test_vswf_translations_array.c index d2b8767..1d8c047 100644 --- a/qpms/tests/test_vswf_translations_array.c +++ b/qpms/tests/test_vswf_translations_array.c @@ -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)) 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 qpms_m_t m1; qpms_l_t l1; qpms_y2mn_p(y1, &m1, &l1); 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){ 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(); } */ - // 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 = ");