``broyden''

 

The modified Broyden algorithm discussed in class is:

1. Trial step --
   • Form trial step
2. Prediction step --
   • Form matrix

     

   • Form vector

     

   • Solve

     

   • Form final prediction

      

Important note: This formulation of the algorithm uses Rho and g in two different ways: to store the trial step and residual (which we called Rho' and g' in lecture) and to store the final prediction step. This becomes important for coding only at the stage where we form the final Rho . At the final step, you should first multiply the contents of Rho by b , and then loop over all of the preceding Rho , multiplying by b and accumulating the final result into Rho .

Beginning with your code ``test_it'', make a new program ``broyden'' with the modifications below.

• Matrix G and vector u
  • Declarations:
  • Residual vector :
  • Formation:
• Solution of G b = u
• Formation of final Rho[it+1]