Method F

cpp/demo/hyperelasticity/main.cpp:159–173

Compute F at current point x

  /// Compute F at current point x
  void F(const Vec x)
  {
    VecGhostUpdateBegin(x, INSERT_VALUES, SCATTER_FORWARD);
    VecGhostUpdateEnd(x, INSERT_VALUES, SCATTER_FORWARD);

    // Assemble b and update ghosts
    std::span b(_b_vec.array());
    std::ranges::fill(b, 0);
    fem::assemble_vector(b, _l);
    VecGhostUpdateBegin(_b, ADD_VALUES, SCATTER_REVERSE);
    VecGhostUpdateEnd(_b, ADD_VALUES, SCATTER_REVERSE);

    // Set bcs
    fem::petsc::set_bc(_b, _bcs, x, -1);
  }

  /// Compute J = F' at current point x
  void J(const Vec, Mat A)