Switch compute from recursive to looped

[f0uriest]

Resolves #1336

[github-actions[bot]]

|             benchmark_name             |         dt(%)          |         dt(s)          |        t_new(s)        |        t_old(s)        | 
| -------------------------------------- | ---------------------- | ---------------------- | ---------------------- | ---------------------- |
 test_build_transform_fft_lowres         |     +0.46 +/- 2.25     | +2.41e-03 +/- 1.18e-02 |  5.25e-01 +/- 7.8e-03  |  5.23e-01 +/- 8.8e-03  |
 test_equilibrium_init_medres            |     -0.14 +/- 1.00     | -5.83e-03 +/- 4.11e-02 |  4.11e+00 +/- 3.6e-02  |  4.12e+00 +/- 2.1e-02  |
 test_equilibrium_init_highres           |     +0.26 +/- 1.12     | +1.42e-02 +/- 5.99e-02 |  5.39e+00 +/- 4.2e-02  |  5.37e+00 +/- 4.3e-02  |
 test_objective_compile_dshape_current   |     -0.13 +/- 4.32     | -5.05e-03 +/- 1.75e-01 |  4.03e+00 +/- 1.3e-01  |  4.04e+00 +/- 1.2e-01  |
 test_objective_compute_dshape_current   |     -1.45 +/- 1.63     | -7.56e-05 +/- 8.50e-05 |  5.14e-03 +/- 6.0e-05  |  5.22e-03 +/- 6.1e-05  |
 test_objective_jac_dshape_current       |     +0.49 +/- 7.25     | +2.06e-04 +/- 3.05e-03 |  4.23e-02 +/- 2.1e-03  |  4.21e-02 +/- 2.2e-03  |
 test_perturb_2                          |     +0.35 +/- 1.58     | +6.75e-02 +/- 3.05e-01 |  1.94e+01 +/- 2.8e-01  |  1.94e+01 +/- 1.1e-01  |
 test_proximal_freeb_jac                 |     -0.02 +/- 0.68     | -1.59e-03 +/- 4.99e-02 |  7.34e+00 +/- 4.4e-02  |  7.35e+00 +/- 2.3e-02  |
 test_solve_fixed_iter                   |     +0.16 +/- 1.75     | +5.07e-02 +/- 5.51e-01 |  3.16e+01 +/- 4.1e-01  |  3.16e+01 +/- 3.7e-01  |
 test_LinearConstraintProjection_build   |     +0.37 +/- 0.82     | +3.77e-02 +/- 8.28e-02 |  1.02e+01 +/- 6.3e-02  |  1.02e+01 +/- 5.4e-02  |
 test_build_transform_fft_midres         |     -0.64 +/- 1.38     | -3.89e-03 +/- 8.40e-03 |  6.04e-01 +/- 5.5e-03  |  6.08e-01 +/- 6.3e-03  |
 test_build_transform_fft_highres        |     -0.22 +/- 7.10     | -2.09e-03 +/- 6.87e-02 |  9.66e-01 +/- 6.8e-02  |  9.68e-01 +/- 1.1e-02  |
 test_equilibrium_init_lowres            |     -3.58 +/- 3.44     | -1.42e-01 +/- 1.36e-01 |  3.83e+00 +/- 2.2e-02  |  3.97e+00 +/- 1.3e-01  |
 test_objective_compile_atf              |     -1.47 +/- 2.01     | -1.23e-01 +/- 1.68e-01 |  8.24e+00 +/- 4.8e-02  |  8.37e+00 +/- 1.6e-01  |
 test_objective_compute_atf              |     +0.56 +/- 1.61     | +8.77e-05 +/- 2.54e-04 |  1.58e-02 +/- 2.0e-04  |  1.58e-02 +/- 1.6e-04  |
 test_objective_jac_atf                  |     -0.48 +/- 2.32     | -9.48e-03 +/- 4.56e-02 |  1.96e+00 +/- 2.6e-02  |  1.97e+00 +/- 3.7e-02  |
 test_perturb_1                          |     -0.25 +/- 2.39     | -3.73e-02 +/- 3.49e-01 |  1.46e+01 +/- 1.9e-01  |  1.46e+01 +/- 2.9e-01  |
 test_proximal_jac_atf                   |     -0.52 +/- 0.87     | -4.38e-02 +/- 7.30e-02 |  8.33e+00 +/- 5.8e-02  |  8.37e+00 +/- 4.4e-02  |
 test_proximal_freeb_compute             |     -0.44 +/- 1.56     | -8.75e-04 +/- 3.13e-03 |  2.00e-01 +/- 2.3e-03  |  2.01e-01 +/- 2.2e-03  |
 test_solve_fixed_iter_compiled          |     +0.31 +/- 3.58     | +6.32e-02 +/- 7.37e-01 |  2.07e+01 +/- 6.9e-01  |  2.06e+01 +/- 2.5e-01  |