UpdateThetaSQUAREM_BBM function

Update theta SQUAREM