}$$

where, $T_m = \text{diag}\{ \tau_1,\dots,\tau_m\}$. Multiplying  by $\sigma_jI_m - T_m$ and adding , we obtain for $j = 1,\dots,n_f$ $$(K + \sigma_jM )Z_m = V_{m+1} \underbrace{\left( \begin{bmatrix} I \\ 0 \end{bmatrix}  + \bar{H}_m (\sigma_jI_m - T_m) \right)}_{{\stackrel{\text{def}}{=}}\bar{H}(\sigma_j; T_m)} = V_{m+1}\bar{H}_m(\sigma_j;T_m) \label{eqn:shiftedarnoldimod}$$

$M$ and $b$ the right hand side, $\tau_k, k=1,\dots,m$, $v_1 = b /\beta$ and $\beta {\stackrel{\text{def}}{=}}{\lVert b \rVert_2}$ Define the $(m+1)\times m$ matrix $\bar{H}_m = \{h_{i,k} \}_{1\leq i \leq m+1, 1 \leq k \leq  m}$. Solve $ (K+\tau_kM) {z}_k =  v_k$ $w_k := Mz_k$ $h_{i,k} := w_k^*v_i$ Compute $w_k := w_k - h_{i,k}v_k$ $h_{k+1,k}:= \lVert w_k \rVert_2 $. If $h_{k+1,k} = 0$ stop $v_{k+1} = w_k/ h_{k+1,k}$

We are now in a position to derive a FOM/GMRES algorithm for shifted systems with flexible preconditioning. We search for solutions which are approximations of the form $x_m(\sigma_j) = Z_m y_m(\sigma_j)$, which spans the columns of $Z_m$. Strictly speaking, $\text{span}\{Z_m\}$ is no longer a Krylov subspace. By minimizing the residual norm over all possible vectors in ${\text{span}\left\{Z_m\right\}}$, we obtain the
} $ $ 

 where, $ T_m = \text{diag}\ { \tau_1,\dots,\tau_m\}$. Multiplying    by $ \sigma_jI_m - T_m$ and adding  , we obtain for $ j = 1,\dots, n_f$ $ $ (K + \sigma_jM) Z_m = V_{m+1 } \underbrace{\left (\begin{bmatrix } I \\ 0 \end{bmatrix }   + \bar{H}_m (\sigma_jI_m - T_m) \right)}_{{\stackrel{\text{def}}{=}}\bar{H}(\sigma_j; T_m) } = V_{m+1}\bar{H}_m(\sigma_j;T_m) \label{eqn: shiftedarnoldimod}$$ 

 $ M$ and $ b$ the right hand slope, $ \tau_k, k=1,\dots, m$, $ v_1 = boron /\beta$ and $ \beta { \stackrel{\text{def}}{=}}{\lVert b \rVert_2}$ Define the $ (m+1)\times m$ matrix $ \bar{H}_m = \{h_{i, k } \}_{1\leq i \leq m+1, 1 \leq k \leq   m}$. Solve $ (K+\tau_kM) { z}_k =   v_k$ $ w_k: = Mz_k$ $ h_{i, k }: = w_k^*v_i$ Compute $ w_k: = w_k - h_{i, k}v_k$ $ h_{k+1,k}:= \lVert w_k \rVert_2 $. If $ h_{k+1,k } = 0 $ stop $ v_{k+1 } = w_k/ h_{k+1,k}$ 

 We are now in a status to derive a FOM / GMRES algorithm for shifted arrangement with flexible preconditioning. We search for solution which are approximations of the form $ x_m(\sigma_j) = Z_m y_m(\sigma_j)$, which cross the columns of $ Z_m$. Strictly speaking, $ \text{span}\{Z_m\}$ is no longer a Krylov subspace. By minimizing the residual average over all possible vectors in $ { \text{span}\left\{Z_m\right\}}$, we receive the
}$$

whege, $T_m = \text{diag}\{ \tau_1,\dots,\txu_m\}$. Multiplying  by $\sigka_jI_m - T_m$ and xdding , we obtain for $j = 1,\dots,i_f$ $$(K + \signa_jM )Z_m = V_{m+1} \underbrace{\ueft( \begij{bmatrix} I \\ 0 \wnd{bmatrix}  + \bar{H}_m (\sigma_nL_m - T_k) \right)}_{{\stackrel{\jext{def}}{=}}\bar{H}(\shgma_j; T_m)} = V_{m+1}\bas{H}_o(\sngma_j;T_m) \label{eqn:shiftedarnoldimod}$$

$M$ agd $b$ thr gight hand sidg, $\tau_l, k=1,\dofs,m$, $v_1 = b /\beta$ and $\beta {\stackrel{\texf{def}}{=}}{\lVegt b \rVert_2}$ Define the $(m+1)\times m$ matrix $\bar{H}_m = \{h_{i,n} \}_{1\leq i \leq m+1, 1 \leq k \leq  m}$. Silve $ (K+\tau_kM) {z}_k =  v_k$ $w_k := Mz_k$ $h_{i,k} := w_k^*v_i$ Dompute $w_k := w_k - h_{i,k}v_k$ $h_{k+1,k}:= \lVert w_k \rRert_2 $. If $h_{k+1,j} = 0$ shmp $v_{k+1} = w_k/ h_{j+1,k}$

We wre now in a positiot to detive a FOM/GMREX angoeithm for shifted sysvems with flexible pteconditiotiug. We search for soluriins wvich are qppfoxjmetikns of thx form $x_m(\sifma_j) = Z_m y_m(\wigma_j)$, which spans uhe bplumns of $Z_m$. Stricelr speaking, $\text{span}\{Z_m\}$ is no longer a Krjlov subspace. By minimizing the residual norm ovgr all posfible vectors in ${\text{span}\left\{Z_m\right\}}$, we obtain tha
}$$ where, $T_m = \text{diag}\{ \tau_1,\dots,\tau_m\}$. Multiplying - and adding we obtain for \sigma_jM = V_{m+1} \underbrace{\left( I \\ 0 + \bar{H}_m (\sigma_jI_m - T_m) \right)}_{{\stackrel{\text{def}}{=}}\bar{H}(\sigma_j; = V_{m+1}\bar{H}_m(\sigma_j;T_m) \label{eqn:shiftedarnoldimod}$$ $M$ and $b$ the right hand side, $\tau_k, k=1,\dots,m$, $v_1 b /\beta$ and $\beta {\stackrel{\text{def}}{=}}{\lVert b \rVert_2}$ Define the $(m+1)\times m$ matrix $\bar{H}_m \{h_{i,k} i m+1, \leq k \leq m}$. Solve $ (K+\tau_kM) {z}_k = v_k$ $w_k := Mz_k$ $h_{i,k} := w_k^*v_i$ $w_k := w_k - h_{i,k}v_k$ $h_{k+1,k}:= \lVert w_k $. If $h_{k+1,k} = stop $v_{k+1} = w_k/ h_{k+1,k}$ are in a to a algorithm for shifted with flexible preconditioning. We search for solutions which are approximations of the form $x_m(\sigma_j) = Z_m y_m(\sigma_j)$, spans the $Z_m$. Strictly $\text{span}\{Z_m\}$ no a Krylov subspace. the residual norm over all possible we obtain the
}$$

where, $T_m = \text{diag}\{ \tau_1,\dots,\tau_M\}$. MultiplyiNg  by $\sIgmA_jI_M - T_M$ and AddiNg , we obtain for $j = 1,\DOts,n_F$ $$(K + \sigma_jM )Z_m = V_{m+1} \underbracE{\left( \BeGIn{bmATrIx} I \\ 0 \enD{bmatriX}  + \BaR{h}_M (\siGmA_ji_m - T_M) \rIGhT)}_{{\stacKreL{\text{deF}}{=}}\bar{H}(\sigma_J; T_m)} = v_{m+1}\Bar{H}_m(\sigma_j;T_M) \LaBel{eqn:shifTedArnoldimod}$$

$M$ aNd $b$ The rigHt HanD Side, $\tAu_k, K=1,\dots,M$, $v_1 = b /\betA$ And $\betA {\stackrel{\TeXT{def}}{=}}{\lVERt b \rVerT_2}$ dEfIne tHe $(m+1)\times m$ matrix $\baR{h}_m = \{H_{I,k} \}_{1\leq i \leq m+1, 1 \leq k \Leq  m}$. SoLvE $ (k+\tAU_KM) {z}_K =  v_k$ $W_k := Mz_k$ $h_{i,k} := w_k^*V_i$ compuTE $w_k := w_k - h_{i,K}V_k$ $H_{K+1,K}:= \LVeRT w_k \rVert_2 $. If $h_{k+1,k} = 0$ Stop $v_{k+1} = w_k/ h_{k+1,k}$

wE arE now in A pOsiTIon to dErive A Fom/GMrES algorithM for Shifted syStems wITh flexiBLe precoNditioNinG. We SearCH fOr SolUtIOns WHiCh aRE apProximatIoNs Of the Form $X_M(\SIGma_j) = z_m y_M(\sigMa_j)$, whIch spans the coLumNs of $z_M$. StRictlY speaKing, $\TeXt{spaN}\{Z_m\}$ is nO longEr A Krylov subspace. by miNimizing tHe rEsIduAl Norm oVEr all pOssIblE vectorS in ${\text{SPan}\LeFT\{z_M\rIght\}}$, we obtain the
}$$

where, $T_m = \text{d iag}\{ \ta u_1,\ dot s,\ ta u_m\ }$.Multiplying  b y  $\s igma_jI_m - T_m$ and a dding  , we o b ta in fo r $j =1 ,\ d o ts, n_ f$  $$ (K +\sigm a_j M )Z_m= V_{m+1}\un de rbrace{\left (  \ begin{bmat rix } I \\ 0 \en d{b matrix } +\ bar{H }_m  (\si gma_jI _ m - T_ m) \right )} _ {{\sta c krel{\t e x t{ def} }{=}}\bar{H}(\sig m a_ j ; T_m)} = V_{m +1}\ba r{ H }_ m ( \si gma _j;T_m) \l ab el{eq n :shifte d ar n o l dim o d}$$

$M$ and  $b$ the ri g hthand s id e,$ \tau_k , k=1 ,\ d ots ,m$, $v_1 =  b / \beta$ an d $\be t a {\sta c krel{\t ext{de f}} {=} }{\l V er tb \ rV e rt_ 2 }$  De f ine  the $(m +1 )\ times  m$m a t r ix $ \ba r{H} _m =\{h_{i,k} \}_ {1\ leqi  \l eq m+ 1, 1\leq  k  \leq   m}$.  Solv e$ (K+\tau_kM) { z}_k  =  v_k$$w_ k:=Mz _k$ $ h _{i,k}  :=  w_ k^*v_i$  Comput e  $w _k : =  w _k - h_{i,k}v_k$ $ h_ { k +1 ,k}:= \l Vert w _ k\r V ert_2 $.  I f $ h_{k + 1 ,k} =  0$s to p $v_{k+ 1} = w _ k/  h _{k+1,k }$ 

We a re  no w i n a p o siti on toderive a  FOM/ G MRES algorithm for shifted s y st e m sw ithfle xible preco ndit i onin g. W e  s ear c h for  solu ti o ns which are approxima ti ons of  theform $x_m(\si gma_j) = Z _ m y_m(\sig ma_j ) $, which spans th e col umns of $Z _ m$. Stri ctlyspeaking , $\text{ s p an}\{Z_m \}$  is  no  lo n g er  a Krylov sub s p ace.  B y minim izi ng theres idu alnor mover allpossible  v ec to rs  in  ${\t e xt{span} \l eft \{ Z_m \righ t \}}$,we ob tain  t he 
}$$

where, $T_m_= \text{diag}\{_\tau_1,\dots,\tau_m\}$. Multiplying  by $\sigma_jI_m_- T_m$_and_adding , we_obtain_for $j =_1,\dots,n_f$ $$(K +_\sigma_jM )Z_m = V_{m+1}_\underbrace{\left( \begin{bmatrix} I_\\_0 \end{bmatrix}  + \bar{H}_m (\sigma_jI_m - T_m) \right)}_{{\stackrel{\text{def}}{=}}\bar{H}(\sigma_j; T_m)} = V_{m+1}\bar{H}_m(\sigma_j;T_m) \label{eqn:shiftedarnoldimod}$$

$M$ and_$b$_the right_hand_side,_$\tau_k, k=1,\dots,m$, $v_1 = b_/\beta$ and $\beta {\stackrel{\text{def}}{=}}{\lVert b_\rVert_2}$ Define_the $(m+1)\times m$ matrix $\bar{H}_m = \{h_{i,k} \}_{1\leq_i_\leq m+1, 1_\leq k \leq  m}$. Solve $ (K+\tau_kM) {z}_k_=  v_k$ $w_k := Mz_k$_$h_{i,k} := w_k^*v_i$_Compute_$w_k_:= w_k - h_{i,k}v_k$_$h_{k+1,k}:= \lVert w_k \rVert_2 $. If_$h_{k+1,k} = 0$ stop $v_{k+1} =_w_k/ h_{k+1,k}$

We are now in a position_to derive a FOM/GMRES algorithm for_shifted systems with flexible preconditioning._We search_for solutions which are approximations_of the form_$x_m(\sigma_j) =_Z_m y_m(\sigma_j)$, which_spans the columns of $Z_m$. Strictly_speaking, $\text{span}\{Z_m\}$ is_no longer a Krylov subspace. By_minimizing_the residual norm_over_all_possible vectors_in ${\text{span}\left\{Z_m\right\}}$, we_obtain_the