https://gcc.gnu.org/bugzilla/show_bug.cgi?id=117874
Richard Biener <rguenth at gcc dot gnu.org> changed: What |Removed |Added ---------------------------------------------------------------------------- Last reconfirmed| |2024-12-02 Status|UNCONFIRMED |ASSIGNED Ever confirmed|0 |1 Assignee|unassigned at gcc dot gnu.org |rguenth at gcc dot gnu.org --- Comment #2 from Richard Biener <rguenth at gcc dot gnu.org> --- Testcase: typedef struct { double real; double imag; } complex; typedef struct { complex e[3][3]; } su3_matrix; void mult_su3_an( su3_matrix *a, su3_matrix *b, su3_matrix *c ) { int j; double a0r,a0i,a1r,a1i,a2r,a2i; double b0r,b0i,b1r,b1i,b2r,b2i; for(j=0;j<3;j++) { a0r=a->e[0][0].real; a0i=a->e[0][0].imag; b0r=b->e[0][j].real; b0i=b->e[0][j].imag; a1r=a->e[1][0].real; a1i=a->e[1][0].imag; b1r=b->e[1][j].real; b1i=b->e[1][j].imag; a2r=a->e[2][0].real; a2i=a->e[2][0].imag; b2r=b->e[2][j].real; b2i=b->e[2][j].imag; c->e[0][j].real = a0r*b0r + a0i*b0i + a1r*b1r + a1i*b1i + a2r*b2r + a2i*b2i; c->e[0][j].imag = a0r*b0i - a0i*b0r + a1r*b1i - a1i*b1r + a2r*b2i - a2i*b2r; a0r=a->e[0][1].real; a0i=a->e[0][1].imag; b0r=b->e[0][j].real; b0i=b->e[0][j].imag; a1r=a->e[1][1].real; a1i=a->e[1][1].imag; b1r=b->e[1][j].real; b1i=b->e[1][j].imag; a2r=a->e[2][1].real; a2i=a->e[2][1].imag; b2r=b->e[2][j].real; b2i=b->e[2][j].imag; c->e[1][j].real = a0r*b0r + a0i*b0i + a1r*b1r + a1i*b1i + a2r*b2r + a2i*b2i; c->e[1][j].imag = a0r*b0i - a0i*b0r + a1r*b1i - a1i*b1r + a2r*b2i - a2i*b2r; a0r=a->e[0][2].real; a0i=a->e[0][2].imag; b0r=b->e[0][j].real; b0i=b->e[0][j].imag; a1r=a->e[1][2].real; a1i=a->e[1][2].imag; b1r=b->e[1][j].real; b1i=b->e[1][j].imag; a2r=a->e[2][2].real; a2i=a->e[2][2].imag; b2r=b->e[2][j].real; b2i=b->e[2][j].imag; c->e[2][j].real = a0r*b0r + a0i*b0i + a1r*b1r + a1i*b1i + a2r*b2r + a2i*b2i; c->e[2][j].imag = a0r*b0i - a0i*b0r + a1r*b1i - a1i*b1r + a2r*b2i - a2i*b2r; } }