Skip to content
Snippets Groups Projects
  1. Jan 05, 2024
  2. Jan 03, 2024
  3. Jan 02, 2024
  4. Dec 04, 2023
  5. Nov 24, 2023
  6. Nov 23, 2023
  7. Nov 21, 2023
  8. Oct 30, 2023
  9. Oct 27, 2023
  10. Oct 26, 2023
  11. Oct 05, 2023
  12. Oct 04, 2023
    • Lisa Julia Nebel's avatar
      Divide the regularization parameter through the scaling parameter in the PN... · 3ea924df
      Lisa Julia Nebel authored
      Divide the regularization parameter through the scaling parameter in the PN solver instead of multiplying
      
      In the PN solver, we add identity*regularization/scaling to the stiffness matrix.
      This has two effects:
      - make the system matrix positive definite
      - penalize large corrections to ensure an energy decrease in each step
      
      The regularization parameter is adapted if either the system matrix is not positive definite
      if there is an energy increase. This works similar to the trust-region radius in the TR solver.
      The correction must be inside the trust-region, this ensures there is an energy decrease
      in each step.
      For each component of the correction, the trust-region radius as well as the regularization parameter
      is multiplied with a scaling factor, because the correction for the displacement works depends
      on the grid size where the correction for the rotation does not.
      
      To align the effect of the scaling parameter, we need to divide through the scaling parameter in
      the PN solver and multiply with it it in the TR solver.
      3ea924df
  13. Jun 20, 2023
  14. Jun 16, 2023
  15. May 24, 2023
  16. May 22, 2023
Loading