Successfully computed bandwidth of Q matrix only. It is 65000 on 64^3 (I expected low).
So run RCM on Q matrix only. Try boost::ublas to do this. Follow Evan.