Compute Q-matrix
q.matrix3(sig, sig.a, sig.b, act.a, act.b, ss)
sig
no descr
sig.a
sig.b
act.a
act.b
ss
n.stadler