solve L problem with fixed G (least square)
(G)
| 809 | return nx.solve(K0, xst) |
| 810 | |
| 811 | def solve_L_bias(G): |
| 812 | """solve L problem with fixed G (least square)""" |
| 813 | xst = ns * nx.dot(G, xt) |
| 814 | return nx.solve(K0, nx.dot(K1.T, xst)) |
| 815 | |
| 816 | def solve_G(L, G0): |
| 817 | """Update G with CG algorithm""" |