def METHOD_NAME(self):
    timemesh = self.timemesh 
    N = len(timemesh)
    D = self.D
    for i in range(N):
        t = timemesh[i]
        b = self.get_right_vector(t)
        #self.uh1[:] =  spsolve(D, b)
        self.uh1[:] = self.ml.METHOD_NAME(b, tol=1e-12, accel='cg').reshape((-1,))
        self.current = i
        if self.current%10 == 0:
            self.show_soultion()
        self.uh0[:] = self.uh1[:]