From aee4d2db2a7d5a81427badfb801f7645237e900e Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Marek=20Ne=C4=8Dada?= Date: Wed, 13 Mar 2019 11:24:09 +0200 Subject: [PATCH] sss3 dump packed matrices Former-commit-id: 5b2f8a4ebe9829d63546ae91e24c25ec5966fe72 --- tests/sss3.c | 11 ++++++++++- 1 file changed, 10 insertions(+), 1 deletion(-) diff --git a/tests/sss3.c b/tests/sss3.c index 23cedcd..deac10a 100644 --- a/tests/sss3.c +++ b/tests/sss3.c @@ -155,9 +155,18 @@ int main() } complex double *S_packed[ss->sym->nirreps]; - for (qpms_iri_t iri = 0; iri < ss->sym->nirreps; ++iri) + for (qpms_iri_t iri = 0; iri < ss->sym->nirreps; ++iri) { S_packed[iri] = qpms_scatsys_irrep_pack_matrix(NULL, S_full, ss, iri); + fprintf(stderr, "--- Packed matrix for irrep %d (%s):\n", (int) iri, ss->sym->irreps[iri].name); + for (size_t row = 0; row < ss->saecv_sizes[iri]; ++row) { + for (size_t col = 0; col < ss->saecv_sizes[iri]; ++col) { + complex double elem = S_packed[iri][row * ss->saecv_sizes[iri] + col]; + fprintf(stderr, "%+.3f+%.3fj ", creal(elem), cimag(elem)); + } + fputc('\n', stderr); + } + } complex double *S_recfull = qpms_scatsys_irrep_unpack_matrix(NULL, S_packed[0], ss, 0, false);