
        $ where $$\begin{aligned}
        T=\frac{1}{h^{2}} \cdot \operatorname{tridiag}(-1,2,-1) \in \mathbb{R}^{p \times p}, \quad F=\frac{1}{h} \cdot \operatorname{tridiag}(0,1,-1) \in \mathbb{R}^{p \times p},
        \end{aligned}$$ and $
        E=\operatorname{diag}\left(1, p+1, \cdots, p^{2}-p+1\right)
      $ in which $\otimes$ denotes Kronecker product and $h=\frac{1}{p+1}$ the discretization meshsize.\

           $p$                              
  -- ----- ---------- ---------- ---------- ----------
           16         32         48         64
     IT    870        3100       6759       11181
     CPU   0.74       85.95      884.22     3409.13
     Err   3.02e-08   8.86e-08   2.30e-05   9.62e-08
     IT    35         38         39         40
     CPU   0.07       0.28       0.93       3.35
     Err   1.06e-6    4.00e-06   5.07e-06   1.68e-05
     IT    26         29         31         32
     CPU   0.02       0.23       0.88       3.30
     Err   1.97e-07   9.11e-07   3.50e-07   6.07e-07
     IT    2          2          2          2
     CPU   0.02       0.21       0.83       3.19
     Err   2.33e-13   1.40e-12   4.72e-12   1.18e-11

  : Numerical results for GMRES to solve Example \[ex1\]

\[tab1\]

In this example, we take $S=\diag(B A^{-1}B^{T}).$ Table \[tab1\] shows the iteration counts and the CPU time of the GMRES with the preconditioners $\mathcal{P}_{D},$ ${\mathcal{P}_1}$ and $\mathcal
$ where $ $ \begin{aligned } 
         T=\frac{1}{h^{2 } } \cdot \operatorname{tridiag}(-1,2,-1) \in \mathbb{R}^{p \times p }, \quad F=\frac{1}{h } \cdot \operatorname{tridiag}(0,1,-1) \in \mathbb{R}^{p \times p }, 
         \end{aligned}$$ and $ 
         E=\operatorname{diag}\left(1, p+1, \cdots, p^{2}-p+1\right) 
       $ in which $ \otimes$ denotes Kronecker product and $ h=\frac{1}{p+1}$ the discretization meshsize.\ 

            $ p$                              
   -- ----- ---------- ---------- ---------- ---------- 
            16          32          48          64 
      IT     870         3100        6759        11181 
      CPU    0.74        85.95       884.22      3409.13 
      Err    3.02e-08    8.86e-08    2.30e-05    9.62e-08 
      IT     35          38          39          40 
      CPU    0.07        0.28        0.93        3.35 
      Err    1.06e-6     4.00e-06    5.07e-06    1.68e-05 
      IT     26          29          31          32 
      CPU    0.02        0.23        0.88        3.30 
      Err    1.97e-07    9.11e-07    3.50e-07    6.07e-07 
      information technology     2           2           2           2 
      central processing unit    0.02        0.21        0.83        3.19 
      Err    2.33e-13    1.40e-12    4.72e-12    1.18e-11 

  : numeral results for GMRES to resolve Example \[ex1\ ] 

 \[tab1\ ] 

 In this exemplar, we take $ S=\diag(B A^{-1}B^{T}).$ Table \[tab1\ ] show the iteration count and the CPU time of the GMRES with the preconditioners $ \mathcal{P}_{D},$ $ { \mathcal{P}_1}$ and $ \mathcal

        $ where $$\begin{aligned}
        T=\frac{1}{h^{2}} \eeot \opxratornzme{tridixg}(-1,2,-1) \in \mathbb{R}^{p \times p}, \quad H=\frax{1}{h} \cdit \operatorname{tridiag}(0,1,-1) \in \mathbh{R}^{p \timew p},
        \end{aligned}$$ akb $
        Z=\o'eratorname{diag}\keft(1, p+1, \cdods, p^{2}-p+1\right)
      $ iu which $\otimes$ denotes Kronecker proquct anc $j=\frac{1}{p+1}$ the dissretpzwtioh meshsize.\

           $p$                              
  -- ----- ---------- ---------- ---------- ----------
           16         32         48         64
     LT    870        3100       6759       11181
     CPU   0.74       85.95      884.22     3409.13
     Err   3.02e-08   8.86e-08   2.30e-05   9.62e-08
     NT    35         38         39         40
     CPU   0.07       0.28       0.93       3.35
     Err   1.06e-6    4.00e-06   5.07e-06   1.68e-05
     IT    26         29         31         32
     CPU   0.02       0.23       0.88       3.30
     Efr   1.97x-07   9.11e-07   3.50e-07   6.07x-07
     IT    2          2          2          2
     CPU   0.02       0.21       0.83       3.19
     Err   2.33e-13   1.40e-12   4.72e-12   1.18e-11

  : Numtricam results for GMRES to wolve Example \[ex1\]

\[tab1\]

In this exaiple, we take $S=\diag(B A^{-1}B^{T}).$ Table \[tab1\] shows the iteradion rojntw wve hhe CPU time of the GMRES with the preconditijhets $\mathcal{P}_{D},$ ${\matmcal{P}_1}$ and $\mathcal
$ where $$\begin{aligned} T=\frac{1}{h^{2}} \cdot \operatorname{tridiag}(-1,2,-1) \in p}, F=\frac{1}{h} \cdot \in \mathbb{R}^{p \times p+1, p^{2}-p+1\right) $ in $\otimes$ denotes Kronecker and $h=\frac{1}{p+1}$ the discretization meshsize.\ $p$ ----- ---------- ---------- ---------- ---------- 16 32 48 64 IT 870 3100 6759 CPU 0.74 85.95 884.22 3409.13 Err 3.02e-08 8.86e-08 2.30e-05 9.62e-08 IT 35 38 40 0.07 0.93 Err 1.06e-6 4.00e-06 5.07e-06 1.68e-05 IT 26 29 31 32 CPU 0.02 0.23 0.88 3.30 Err 9.11e-07 3.50e-07 6.07e-07 IT 2 2 2 2 0.02 0.21 0.83 3.19 2.33e-13 1.40e-12 4.72e-12 1.18e-11 : results GMRES to Example \[tab1\] this example, we $S=\diag(B A^{-1}B^{T}).$ Table \[tab1\] shows the iteration counts and the CPU time of the GMRES with the $\mathcal{P}_{D},$ ${\mathcal{P}_1}$

        $ where $$\begin{aligned}
        T=\frac{1}{h^{2}} \cdOt \operatorName{tRidIag}(-1,2,-1) \In \MathBb{R}^{p \Times p}, \quad F=\fraC{1}{H} \cdoT \operatorname{tridiag}(0,1,-1) \in \MathbB{R}^{P \TimeS P},
        \eNd{aliGned}$$ and $
        e=\OpERAtoRnAmE{diAg}\LEfT(1, p+1, \cdoTs, p^{2}-P+1\right)
      $ iN which $\otimEs$ dEnOtes KroneckeR PrOduct and $h=\fRac{1}{P+1}$ the discretiZatIon mesHsIze.\

           $P$                              
  -- ----- ---------- ---------- ---------- ----------
           16         32         48         64
     iT    870        3100       6759       11181
     CPU   0.74       85.95      884.22     3409.13
     err   3.02E-08   8.86e-08   2.30e-05   9.62e-08
     It    35         38         39         40
     CPU   0.07       0.28       0.93       3.35
     ErR   1.06E-6    4.00e-06   5.07e-06   1.68e-05
     IT    26         29         31         32
     cPU   0.02       0.23       0.88       3.30
     Err   1.97e-07   9.11e-07   3.50e-07   6.07E-07
     It    2          2          2          2
     cPU   0.02       0.21       0.83       3.19
     Err   2.33E-13   1.40E-12   4.72e-12   1.18e-11

  : NumeRICaL resUlts for GMRES to solVE EXAmple \[ex1\]

\[tab1\]

In thIs examPlE, We TAKe $S=\DiaG(B A^{-1}B^{T}).$ Table \[TaB1\] showS The iterATiON COunTS and the CPU timE of the GMRES WIth The preCoNdiTIoners $\MathcAl{p}_{d},$ ${\maThcal{P}_1}$ and $\maThcaL

        $ where $$\begin{ aligned}
      T =\f ra c{1} {h^{ 2}} \cdot \ope r ator name{tridiag}(-1,2,-1)  \in\m a thbb { R} ^{p \ times p } ,\ q uad  F =\ fra c{ 1 }{ h} \c dot  \opera torname{tr idi ag }(0,1,-1) \i n  \ mathbb{R}^ {p\times p},
      \e nd {al i gned} $$and $ 
      E=\o peratorna me { diag}\ l eft(1,p + 1,  \cd ots, p^{2}-p+1\ri g ht ) 
      $ in wh ich $\ ot i me s $  de not es Kroneck er  prod u ct and$ h= \ f r ac{ 1 }{p+1}$ the d iscretizati o n m eshsiz e. \
         $ p$                          
  - -  ------ ------- -- --- --- --- - -- - -- -- ---  - - --- - -- --         1 6      32   48       6 4
     IT   870      3100     675 9     11181 
    C PU   0.74      85. 95      8 84. 22    3409. 1 3
    Er r  3.02e- 08   8. 8 6e- 08 2. 30e-05   9.62e-08  I T    35      3 8       3 9     4 0
    CP U   0.07      0. 2 8     0.9 3     3. 35    Err   1. 06e-6   4.00e -06 5.07e-06   1.6 8 e-05
     IT 2 6       2 9         3 1      32   C P U   0 .02   0 .23       0.88     3.30
    E rr   1.97e-07    9.11e-0 7  3.50e-0 7  6 .0 7 e-07
     IT   2         2       2         2     CPU0.02     0.2 1   0. 83       3.19    Er r   2.3 3e- 13   1. 40e -12   4.7 2e -12   1.1 8e-11

 :  N um er ica l res u lts forGM RES  t o s olveE xample  \[ex 1\]
\ [t a b1\ ]

In t h is e xamp le ,we t ake  $ S=\di ag(B A^{ -1}B^{T }).$ Tabl e \ [ tab1 \]  s hows th e iteration c ou nts and th eCPU  timeo f  the GMR ES with the preconditio n ers $\m ath cal{P }_{D },$ ${\ma thc al{P}_ 1}$ and $\ mathca l

 _ _   _$ where_$$\begin{aligned}
_ __  _ T=\frac{1}{h^{2}} \cdot_\operatorname{tridiag}(-1,2,-1) \in \mathbb{R}^{p \times_p}, \quad F=\frac{1}{h}_\cdot_\operatorname{tridiag}(0,1,-1) \in \mathbb{R}^{p \times p},
        \end{aligned}$$ and_$
_ ___   E=\operatorname{diag}\left(1, p+1,_\cdots, p^{2}-p+1\right)
   _ _$ in which $\otimes$ denotes Kronecker product and_$h=\frac{1}{p+1}$_the discretization meshsize.\

_        _ $p$    _  ___   _     _     _    
  --_----- ---------- ---------- ---------- ----------
 _    _ _  16  _  _ _ 32 _     _ 48 _     __64
  ___IT _  870__ __   3100  __   6759   _   11181
__   CPU_  0.74  _    85.95 _ _ _884.22     3409.13
     Err_  3.02e-08   8.86e-08_  2.30e-05 _ 9.62e-08
__  _IT_ _ 35     _ _ 38    _   _39_     _  40
   _ CPU   0.07___    0.28_     _0.93  _    3.35
  __ Err   1.06e-6 __ 4.00e-06   5.07e-06 __1.68e-05
_    IT_   26  _      29  _  _   31 ___     32
    _CPU _ 0.02 _     0.23    _  0.88    _  3.30
     Err   1.97e-07_  9.11e-07   3.50e-07  _6.07e-07
     IT _  2 __      _2___  _  _ 2 _ _  _ _ 2
 __  CPU   0.02       0.21     _ 0.83    _  3.19
  __ Err _ 2.33e-13  _1.40e-12   4.72e-12   1.18e-11

 _: Numerical results for GMRES to solve Example \[ex1\]

\[tab1\]

In_this_example, we_take $S=\diag(B A^{-1}B^{T}).$ Table_\[tab1\] shows_the iteration_counts and the_CPU time_of the_GMRES with_the preconditioners $\mathcal{P}_{D},$ ${\mathcal{P}_1}$ and $\mathcal