add LB class: compute LB operator and solve eigensystem