<!DOCTYPE html>
<html>
  <head>
    <meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
  </head>
  <body>
    <p>Thanks for your response. You are correct. I overlooked this
      step.</p>
    <p>Now I am trying to correct my "shell matrix approach" for
      ex1f.F90 of snes solver
(<a class="moz-txt-link-freetext" href="https://github.com/hyharry/small_petsc_test/blob/master/test_shell_jac/ex1f.F90">https://github.com/hyharry/small_petsc_test/blob/master/test_shell_jac/ex1f.F90</a>).
      I realized that I need to record the base vector X in the context
      of shell matrix and then use this info to carry MyMult. However,
      the context cannot be obtained through MatShellGetContext(). Here
      are the critical parts of my code.</p>
    <p>       INTERFACE MatCreateShell<br>
               SUBROUTINE
      MatCreateShell(comm,mloc,nloc,m,n,ctx,mat,ierr)<br>
                 USE solver_context<br>
                 MPI_Comm :: comm<br>
                 PetscInt :: mloc,nloc,m,n<br>
                 Vec :: ctx<br>
                 Mat :: mat<br>
                 PetscErrorCode :: ierr<br>
               END SUBROUTINE MatCreateShell<br>
             END INTERFACE MatCreateShell<br>
      <br>
             INTERFACE MatShellSetContext<br>
               SUBROUTINE MatShellSetContext(mat,ctx,ierr)<br>
                 USE solver_context<br>
                 Mat :: mat<br>
                 !TYPE(MatCtx) :: ctx<br>
                 Vec :: ctx<br>
                 PetscErrorCode :: ierr<br>
               END SUBROUTINE MatShellSetContext<br>
             END INTERFACE MatShellSetContext<br>
      <br>
             INTERFACE MatShellGetContext<br>
               SUBROUTINE MatShellGetContext(mat,ctx,ierr)<br>
                 USE solver_context<br>
                 Mat :: mat<br>
                 Vec, Pointer :: ctx<br>
                 PetscErrorCode :: ierr<br>
               END SUBROUTINE MatShellGetContext<br>
             END INTERFACE MatShellGetContext</p>
    <p>in my FormShellJacobian() I did</p>
    <p>subroutine FormJacobianShell(snes,X,jac,B,dummy,ierr)</p>
    <p>......<br>
      <br>
        call MatShellSetContext(jac,X,ierr)</p>
    <p>......</p>
    <p>Then in MyMult() I tried to recover this context by <br>
    </p>
    <p>call MatShellGetContext(J,x,ierr)</p>
    <p>call VecView(x,PETSC_VIEWER_STDOUT_WORLD,ierr)<br>
    </p>
    <p>Then the program failed with <br>
    </p>
    <p>[0]PETSC ERROR: --------------------- Error Message
      --------------------------------------------------------------<br>
      [0]PETSC ERROR: Null argument, when expecting valid pointer<br>
      [0]PETSC ERROR: Null Pointer: Parameter # 1<br>
    </p>
    <p>In MyMult, I actually defined x to be a pointer. So I am confused
      here. </p>
    <p>Best regards,</p>
    <p>Yi <br>
    </p>
    <div class="moz-cite-prefix">On 1/31/24 03:18, Barry Smith wrote:<br>
    </div>
    <blockquote type="cite"
      cite="mid:6C071B22-8D59-4DB8-BCF9-6F31300A648B@petsc.dev">
      <meta http-equiv="content-type" content="text/html; charset=UTF-8">
      <div><br>
      </div>
         It is not running an extra KSP iteration. This "extra" matmult
      is normal and occurs in many of the SNESLineSearchApply_*
      functions, for example, <a
href="https://petsc.org/release/src/snes/linesearch/impls/bt/linesearchbt.c.html#SNESLineSearchApply_BT"
        moz-do-not-send="true" class="moz-txt-link-freetext">https://petsc.org/release/src/snes/linesearch/impls/bt/linesearchbt.c.html#SNESLineSearchApply_BT</a> It
      is used to decide if the Newton step results in sufficient
      decrease of the function value.
      <div><br>
      </div>
      <div>  Barry<br>
        <div><br>
        </div>
        <div><br id="lineBreakAtBeginningOfMessage">
          <div><br>
            <blockquote type="cite">
              <div>On Jan 30, 2024, at 3:19 PM, Yi Hu
                <a class="moz-txt-link-rfc2396E" href="mailto:y.hu@mpie.de"><y.hu@mpie.de></a> wrote:</div>
              <br class="Apple-interchange-newline">
              <div>
                <meta http-equiv="Content-Type"
                  content="text/html; charset=UTF-8">
                <div>
                  <p>Hello Barry,</p>
                  <p>Thanks for your reply. The monitor options are
                    fine. I actually meant my modification of snes
                    tutorial ex1f.F90 does not work and has some
                    unexpected behavior. I basically wanted to test if I
                    can use a shell matrix as my jacobian (code is here
                    <a class="moz-txt-link-freetext"
href="https://github.com/hyharry/small_petsc_test/blob/master/test_shell_jac/ex1f.F90"
                      moz-do-not-send="true">https://github.com/hyharry/small_petsc_test/blob/master/test_shell_jac/ex1f.F90</a>).
                    After compile my modified version and run with these
                    monitor options, it gives me the following,<br>
                  </p>
                  <p> ( in rhs )<br>
                     ( leave rhs )<br>
                      0 SNES Function norm 6.041522986797e+00 <br>
                     ++++++++++++ in jac shell +++++++++++<br>
                        0 KSP Residual norm 6.041522986797e+00 <br>
                     === start mymult ===<br>
                     === done mymult ===<br>
                        1 KSP Residual norm 5.065392549852e-16 <br>
                      Linear solve converged due to CONVERGED_RTOL
                    iterations 1<br>
                     === start mymult ===<br>
                     === done mymult ===<br>
                     ( in rhs )<br>
                     ( leave rhs )<br>
                      1 SNES Function norm 3.512662245652e+00 <br>
                     ++++++++++++ in jac shell +++++++++++<br>
                        0 KSP Residual norm 3.512662245652e+00 <br>
                     === start mymult ===<br>
                     === done mymult ===<br>
                        1 KSP Residual norm 6.230314124713e-16 <br>
                      Linear solve converged due to CONVERGED_RTOL
                    iterations 1<br>
                     === start mymult ===<br>
                     === done mymult ===<br>
                     ( in rhs )<br>
                     ( leave rhs )<br>
                     ( in rhs )<br>
                     ( leave rhs )<br>
                      2 SNES Function norm 8.969285922373e-01 <br>
                     ++++++++++++ in jac shell +++++++++++<br>
                        0 KSP Residual norm 8.969285922373e-01 <br>
                     === start mymult ===<br>
                     === done mymult ===<br>
                        1 KSP Residual norm 0.000000000000e+00 <br>
                      Linear solve converged due to CONVERGED_ATOL
                    iterations 1<br>
                     === start mymult ===<br>
                     === done mymult ===<br>
                     ( in rhs )<br>
                     ( leave rhs )<br>
                     ( in rhs )<br>
                     ( leave rhs )<br>
                      3 SNES Function norm 4.863816734540e-01 <br>
                     ++++++++++++ in jac shell +++++++++++<br>
                        0 KSP Residual norm 4.863816734540e-01 <br>
                     === start mymult ===<br>
                     === done mymult ===<br>
                        1 KSP Residual norm 0.000000000000e+00 <br>
                      Linear solve converged due to CONVERGED_ATOL
                    iterations 1<br>
                     === start mymult ===<br>
                     === done mymult ===<br>
                     ( in rhs )<br>
                     ( leave rhs )<br>
                     ( in rhs )<br>
                     ( leave rhs )<br>
                      4 SNES Function norm 3.512070785520e-01 <br>
                     ++++++++++++ in jac shell +++++++++++<br>
                        0 KSP Residual norm 3.512070785520e-01 <br>
                     === start mymult ===<br>
                     === done mymult ===<br>
                        1 KSP Residual norm 0.000000000000e+00 <br>
                      Linear solve converged due to CONVERGED_ATOL
                    iterations 1<br>
                     === start mymult ===<br>
                     === done mymult ===<br>
                     ( in rhs )<br>
                     ( leave rhs )<br>
                     ( in rhs )<br>
                     ( leave rhs )<br>
                      5 SNES Function norm 2.769700293115e-01 <br>
                     ++++++++++++ in jac shell +++++++++++<br>
                        0 KSP Residual norm 2.769700293115e-01 <br>
                     === start mymult ===<br>
                     === done mymult ===<br>
                        1 KSP Residual norm 1.104778916974e-16 <br>
                      Linear solve converged due to CONVERGED_RTOL
                    iterations 1<br>
                     === start mymult ===<br>
                     === done mymult ===<br>
                     ( in rhs )<br>
                     ( leave rhs )<br>
                     ( in rhs )<br>
                     ( leave rhs )<br>
                      6 SNES Function norm 2.055345318150e-01 <br>
                     ++++++++++++ in jac shell +++++++++++<br>
                        0 KSP Residual norm 2.055345318150e-01 <br>
                     === start mymult ===<br>
                     === done mymult ===<br>
                        1 KSP Residual norm 0.000000000000e+00 <br>
                      Linear solve converged due to CONVERGED_ATOL
                    iterations 1<br>
                     === start mymult ===<br>
                     === done mymult ===<br>
                     ( in rhs )<br>
                     ( leave rhs )<br>
                     ( in rhs )<br>
                     ( leave rhs )<br>
                      7 SNES Function norm 1.267482220786e-01 <br>
                     ++++++++++++ in jac shell +++++++++++<br>
                        0 KSP Residual norm 1.267482220786e-01 <br>
                     === start mymult ===<br>
                     === done mymult ===<br>
                        1 KSP Residual norm 1.498679601680e-17 <br>
                      Linear solve converged due to CONVERGED_RTOL
                    iterations 1<br>
                     === start mymult ===<br>
                     === done mymult ===<br>
                     ( in rhs )<br>
                     ( leave rhs )<br>
                     ( in rhs )<br>
                     ( leave rhs )<br>
                      8 SNES Function norm 3.468150619264e-02 <br>
                     ++++++++++++ in jac shell +++++++++++<br>
                        0 KSP Residual norm 3.468150619264e-02 <br>
                     === start mymult ===<br>
                     === done mymult ===<br>
                        1 KSP Residual norm 5.944160522951e-18 <br>
                      Linear solve converged due to CONVERGED_RTOL
                    iterations 1<br>
                     === start mymult ===<br>
                     === done mymult ===<br>
                     ( in rhs )<br>
                     ( leave rhs )<br>
                     ( in rhs )<br>
                     ( leave rhs )<br>
                     ( in rhs )<br>
                     ( leave rhs )<br>
                     ( in rhs )<br>
                     ( leave rhs )<br>
                     ( in rhs )<br>
                     ( leave rhs )<br>
                     ( in rhs )<br>
                     ( leave rhs )<br>
                     ( in rhs )<br>
                     ( leave rhs )<br>
                     ( in rhs )<br>
                     ( leave rhs )<br>
                     ( in rhs )<br>
                     ( leave rhs )<br>
                     ( in rhs )<br>
                     ( leave rhs )<br>
                     ( in rhs )<br>
                     ( leave rhs )<br>
                     ( in rhs )<br>
                     ( leave rhs )<br>
                     ( in rhs )<br>
                     ( leave rhs )<br>
                     ( in rhs )<br>
                     ( leave rhs )<br>
                     === start mymult ===<br>
                     === done mymult ===<br>
                    Nonlinear solve did not converge due to
                    DIVERGED_LINE_SEARCH iterations 8<br>
                    Number of SNES iterations =     8<br>
                  </p>
                  <p>After each "Linear solve converged due to
                    CONVERGED_ATOL iterations", the code starts to do
                    mymult again. So I thought it did an extra
                    (unwanted) KSP iteration. I would like to ask if
                    this extra iteration could be disabled, or maybe I
                    am wrong about it.</p>
                  <p>Best regards,</p>
                  <p>Yi<br>
                  </p>
                  <div class="moz-cite-prefix">On 1/30/24 18:35, Barry
                    Smith wrote:<br>
                  </div>
                  <blockquote type="cite"
cite="mid:EAE20968-7A66-449E-8A14-6FBAD06CE0C3@petsc.dev">
                    <meta http-equiv="content-type"
                      content="text/html; charset=UTF-8">
                    <div><br>
                    </div>
                      How do I see a difference? What does "<span
                      style="font-family: Calibri, sans-serif;">hence
                      ruin my previous converged KSP result" mean? A
                      different answer at the end of the KSP solve?</span>
                    <div><br>
                    </div>
                    <div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">$ ./joe >
                          joe.basic</span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">~/Src/petsc/src/ksp/ksp/tutorials</span><span
style="font-variant-ligatures: no-common-ligatures; color: #c814c9"><b>
                            (barry/2023-09-15/fix-log-pcmpi=)</b></span><span
style="font-variant-ligatures: no-common-ligatures"> arch-fix-log-pcmpi</span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">$ ./joe -ksp_monitor
                          -ksp_converged_reason -snes_monitor >
                          joe.monitor</span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">~/Src/petsc/src/ksp/ksp/tutorials</span><span
style="font-variant-ligatures: no-common-ligatures; color: #c814c9"><b>
                            (barry/2023-09-15/fix-log-pcmpi=)</b></span><span
style="font-variant-ligatures: no-common-ligatures"> arch-fix-log-pcmpi</span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">$ diff joe.basic
                          joe.monitor </span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">0a1,36</span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>   0 SNES
                          Function norm 6.041522986797e+00 </span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>     0 KSP
                          Residual norm 6.041522986797e+00 </span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>     1 KSP
                          Residual norm 5.065392549852e-16 </span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>   Linear solve
                          converged due to CONVERGED_RTOL_NORMAL
                          iterations 1</span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>   1 SNES
                          Function norm 3.512662245652e+00 </span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>     0 KSP
                          Residual norm 3.512662245652e+00 </span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>     1 KSP
                          Residual norm 6.230314124713e-16 </span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>   Linear solve
                          converged due to CONVERGED_RTOL_NORMAL
                          iterations 1</span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>   2 SNES
                          Function norm 8.969285922373e-01 </span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>     0 KSP
                          Residual norm 8.969285922373e-01 </span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>     1 KSP
                          Residual norm 0.000000000000e+00 </span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>   Linear solve
                          converged due to CONVERGED_RTOL_NORMAL
                          iterations 1</span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>   3 SNES
                          Function norm 4.863816734540e-01 </span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>     0 KSP
                          Residual norm 4.863816734540e-01 </span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>     1 KSP
                          Residual norm 0.000000000000e+00 </span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>   Linear solve
                          converged due to CONVERGED_RTOL_NORMAL
                          iterations 1</span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>   4 SNES
                          Function norm 3.512070785520e-01 </span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>     0 KSP
                          Residual norm 3.512070785520e-01 </span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>     1 KSP
                          Residual norm 0.000000000000e+00 </span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>   Linear solve
                          converged due to CONVERGED_RTOL_NORMAL
                          iterations 1</span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>   5 SNES
                          Function norm 2.769700293115e-01 </span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>     0 KSP
                          Residual norm 2.769700293115e-01 </span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>     1 KSP
                          Residual norm 1.104778916974e-16 </span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>   Linear solve
                          converged due to CONVERGED_RTOL_NORMAL
                          iterations 1</span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>   6 SNES
                          Function norm 2.055345318150e-01 </span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>     0 KSP
                          Residual norm 2.055345318150e-01 </span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>     1 KSP
                          Residual norm 1.535110861002e-17 </span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>   Linear solve
                          converged due to CONVERGED_RTOL_NORMAL
                          iterations 1</span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>   7 SNES
                          Function norm 1.267482220786e-01 </span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>     0 KSP
                          Residual norm 1.267482220786e-01 </span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>     1 KSP
                          Residual norm 1.498679601680e-17 </span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>   Linear solve
                          converged due to CONVERGED_RTOL_NORMAL
                          iterations 1</span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>   8 SNES
                          Function norm 3.468150619264e-02 </span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>     0 KSP
                          Residual norm 3.468150619264e-02 </span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>     1 KSP
                          Residual norm 5.944160522951e-18 </span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures">>   Linear solve
                          converged due to CONVERGED_RTOL_NORMAL
                          iterations 1</span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures"><br>
                        </span></div>
                      <div
style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 16px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span
style="font-variant-ligatures: no-common-ligatures"><br>
                        </span></div>
                      <div><br>
                        <blockquote type="cite">
                          <div>On Jan 30, 2024, at 11:19 AM, Yi Hu <a
                              class="moz-txt-link-rfc2396E"
                              href="mailto:y.hu@mpie.de"
                              moz-do-not-send="true"><y.hu@mpie.de></a>
                            wrote:</div>
                          <br class="Apple-interchange-newline">
                          <div>
                            <meta charset="UTF-8">
                            <div class="WordSection1"
style="page: WordSection1; caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">
                              <div
style="margin: 0in 0in 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;"><span
                                  lang="EN-US">Dear PETSc team,<o:p></o:p></span></div>
                              <div
style="margin: 0in 0in 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;"><span
                                  lang="EN-US"><o:p> </o:p></span></div>
                              <div
style="margin: 0in 0in 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;"><span
                                  lang="EN-US">I am still trying to sort
                                  out my previous thread<span
                                    class="Apple-converted-space"> </span><a
href="https://lists.mcs.anl.gov/pipermail/petsc-users/2024-January/050079.html"
style="color: rgb(149, 79, 114); text-decoration: underline;"
                                    moz-do-not-send="true"
                                    class="moz-txt-link-freetext">https://lists.mcs.anl.gov/pipermail/petsc-users/2024-January/050079.html</a><span
                                    class="Apple-converted-space"> </span>using
                                  a minimal working example. However, I
                                  encountered another problem. Basically
                                  I combined the basic usage of SNES
                                  solver and shell matrix and tried to
                                  make it work. The jacobian of my snes
                                  is replaced by a customized
                                  MATOP_MULT. The minimal example code
                                  can be viewed here<span
                                    class="Apple-converted-space"> </span><a
href="https://github.com/hyharry/small_petsc_test/blob/master/test_shell_jac/ex1f.F90"
style="color: rgb(149, 79, 114); text-decoration: underline;"
                                    moz-do-not-send="true"
                                    class="moz-txt-link-freetext">https://github.com/hyharry/small_petsc_test/blob/master/test_shell_jac/ex1f.F90</a><o:p></o:p></span></div>
                              <div
style="margin: 0in 0in 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;"><span
                                  lang="EN-US"><o:p> </o:p></span></div>
                              <div
style="margin: 0in 0in 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;"><span
                                  lang="EN-US">When running with
                                  -ksp_monitor -ksp_converged_reason, it
                                  shows an extra mymult step, and hence
                                  ruin my previous converged KSP result.
                                  Implement a customized converged
                                  call-back also does not help. I am
                                  wondering how to skip this extra ksp
                                  iteration. Could anyone help me on
                                  this?<o:p></o:p></span></div>
                              <div
style="margin: 0in 0in 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;"><span
                                  lang="EN-US"><o:p> </o:p></span></div>
                              <div
style="margin: 0in 0in 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;"><span
                                  lang="EN-US">Thanks for your help.<o:p></o:p></span></div>
                              <div
style="margin: 0in 0in 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;"><span
                                  lang="EN-US"><o:p> </o:p></span></div>
                              <div
style="margin: 0in 0in 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;"><span
                                  lang="EN-US">Best wishes,<o:p></o:p></span></div>
                              <div
style="margin: 0in 0in 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;"><span
                                  lang="EN-US">Yi<o:p></o:p></span></div>
                            </div>
                            <br
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">
                            <br
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">
                            <hr
style="font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;"><span
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;">-------------------------------------------------</span><br
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">
                            <span
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;">Stay up to date and follow us on LinkedIn, Twitter and YouTube.</span><br
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">
                            <br
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">
                            <span
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;">Max-Planck-Institut für Eisenforschung GmbH</span><br
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">
                            <span
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;">Max-Planck-Straße 1</span><br
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">
                            <span
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;">D-40237 Düsseldorf</span><br
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">
                            <span
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;"> </span><br
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">
                            <span
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;">Handelsregister B 2533 </span><br
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">
                            <span
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;">Amtsgericht Düsseldorf</span><br
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">
                            <span
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;"> </span><br
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">
                            <span
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;">Geschäftsführung</span><br
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">
                            <span
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;">Prof. Dr. Gerhard Dehm</span><br
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">
                            <span
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;">Prof. Dr. Jörg Neugebauer</span><br
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">
                            <span
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;">Prof. Dr. Dierk Raabe</span><br
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">
                            <span
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;">Dr. Kai de Weldige</span><br
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">
                            <span
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;"> </span><br
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">
                            <span
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;">Ust.-Id.-Nr.: DE 11 93 58 514 </span><br
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">
                            <span
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;">Steuernummer: 105 5891 1000</span><br
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">
                            <br
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">
                            <br
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">
                            <span
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;">Please consider that invitations and e-mails of our institute are </span><br
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">
                            <span
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;">only valid if they end with …@mpie.de. </span><br
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">
                            <span
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;">If you are not sure of the validity please contact </span><a
                              href="mailto:rco@mpie.de"
style="color: rgb(149, 79, 114); text-decoration: underline; font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;"
                              moz-do-not-send="true"
                              class="moz-txt-link-freetext">rco@mpie.de</a><br
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">
                            <br
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">
                            <span
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;">Bitte beachten Sie, dass Einladungen zu Veranstaltungen und E-Mails</span><br
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">
                            <span
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;">aus unserem Haus nur mit der Endung …@mpie.de gültig sind. </span><br
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">
                            <span
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;">In Zweifelsfällen wenden Sie sich bitte an </span><a
                              href="mailto:rco@mpie.de"
style="color: rgb(149, 79, 114); text-decoration: underline; font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;"
                              moz-do-not-send="true"
                              class="moz-txt-link-freetext">rco@mpie.de</a><br
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">
                            <span
style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;">-------------------------------------------------</span></div>
                        </blockquote>
                      </div>
                      <br>
                    </div>
                  </blockquote>
                  <br>
                  <br>
                  <hr>
                  -------------------------------------------------<br>
Stay up to date and follow us on LinkedIn, Twitter and YouTube.<br>
                  <br>
                  Max-Planck-Institut für Eisenforschung GmbH<br>
                  Max-Planck-Straße 1<br>
                  D-40237 Düsseldorf<br>
                   <br>
                  Handelsregister B 2533 <br>
                  Amtsgericht Düsseldorf<br>
                   <br>
                  Geschäftsführung<br>
                  Prof. Dr. Gerhard Dehm<br>
                  Prof. Dr. Jörg Neugebauer<br>
                  Prof. Dr. Dierk Raabe<br>
                  Dr. Kai de Weldige<br>
                   <br>
                  Ust.-Id.-Nr.: DE 11 93 58 514 <br>
                  Steuernummer: 105 5891 1000<br>
                  <br>
                  <br>
Please consider that invitations and e-mails of our institute are <br>
                  only valid if they end with …@mpie.de. <br>
If you are not sure of the validity please contact <a class="moz-txt-link-abbreviated" href="mailto:rco@mpie.de">rco@mpie.de</a><br>
                  <br>
Bitte beachten Sie, dass Einladungen zu Veranstaltungen und E-Mails<br>
aus unserem Haus nur mit der Endung …@mpie.de gültig sind. <br>
                  In Zweifelsfällen wenden Sie sich bitte an <a class="moz-txt-link-abbreviated" href="mailto:rco@mpie.de">rco@mpie.de</a><br>
                  -------------------------------------------------<br>
                </div>
              </div>
            </blockquote>
          </div>
          <br>
        </div>
      </div>
    </blockquote>
  <BR />
<BR />
<HR />
-------------------------------------------------<BR />
Stay up to date and follow us on LinkedIn, Twitter and YouTube.<BR />
<BR />
Max-Planck-Institut für Eisenforschung GmbH<BR />
Max-Planck-Straße 1<BR />
D-40237 Düsseldorf<BR />
 <BR />
Handelsregister B 2533 <BR />
Amtsgericht Düsseldorf<BR />
 <BR />
Geschäftsführung<BR />
Prof. Dr. Gerhard Dehm<BR />
Prof. Dr. Jörg Neugebauer<BR />
Prof. Dr. Dierk Raabe<BR />
Dr. Kai de Weldige<BR />
 <BR />
Ust.-Id.-Nr.: DE 11 93 58 514 <BR />
Steuernummer: 105 5891 1000<BR />
<BR />
<BR />
Please consider that invitations and e-mails of our institute are <BR />
only valid if they end with …@mpie.de. <BR />
If you are not sure of the validity please contact rco@mpie.de<BR />
<BR />
Bitte beachten Sie, dass Einladungen zu Veranstaltungen und E-Mails<BR />
aus unserem Haus nur mit der Endung …@mpie.de gültig sind. <BR />
In Zweifelsfällen wenden Sie sich bitte an rco@mpie.de<BR />
-------------------------------------------------<BR />
</body>
</html>