Inspired by the simple two-step threshold dynamics algorithm which iteratively does convolution and thresholding to simulate the motion of grain boundaries, we developed an algorithm to approach the minimal compliance problem in topology optimization with proved unconditional energy stability. However, the algorithm alone usually locates a local minimum which could be improved further. We discuss some possible remedies to get a better result that is comparable to the results obtained by existing optimization methods.