#include #include #include "cpx.hh" #include "su2.hh" using namespace std ; using namespace cpx4 ; using namespace su24 ; int main() { using atype = float ; using btype = cpx ; auto i = btype(0,1) ; auto sx = su2(0, 1, 1, 0) ; auto sy = su2(0, -i, i, 0) ; auto sum = sx+sy ; cout << ((i+sy)^sx) << endl ; // ... exactly ! raise an SU(2) to the power of an SU(2) // these objects you have in the tool-kit are smart // ... note the pencil_on_paper ease, i+sy are different // types, but are properly handled }