Added method to calculate Wilson loops in ToricCode_1_q class