Structure-Preserving Model Order Reduction for Index One Port-Hamiltonian Descriptor Systems
We develop optimization-based structure-preserving model order reduction (MOR) methods for port-Hamiltonian (pH) descriptor systems of differentiation index one. Descriptor systems in pH form permit energy-based modeling and intuitive coupling of physical systems across different physical domains, scales, and accuracies. This makes pH models well-suited building-blocks for component-wise modeling of large system networks. In this context, it is often necessary to preserve the pH structure during MOR. We discuss current projection-based and structure-preserving MOR algorithms for pH systems and present a new optimization-based framework for that task. The benefits of our method include a simplified treatment of algebraic constraints and often a higher accuracy of the resulting reduced-order model, which is demonstrated by several numerical examples.



\subsection{PH-preserving MOR techniques}%
Traditional methods which directly retain the pH form in the reduction process are based on Galerkin projections, see, e.g., \cite{Antoulas2020, MehU22} for surveys. The original state $x(t)$ is approximated by ${x(t)\approx V x_r}$,
where the columns of ${V\in\R^{n\times r}}$ form a basis for a suitably chosen subspace of dimension $r$. For instance, in the pH-ODE case with $\pE= I_n$, the ROM coefficient matrices are computed by
\begin{equation*}
	\begin{array}{lll}
		\rJ=U^\T\pJ U, & \rG = U^\T\pG, & \rQ = V^\T\pQ V, \\ \rR=U^\T\pR U, & \rP = U^\T\pP, & \rN=\pN, \quad \rS=\pS,
	\end{array}
\end{equation*}
where ${U = QV(V^\T QV)^{-1}}$, which clearly enforces structure preservation. In \cite{Gugercin2012} an adaptation (IRKA-PH) of the well-known iterative rational Krylov algorithm (IRKA) was proposed that iteratively updates $V$ to fulfill a subset of $\cH_2$ optimality conditions via tangential interpolation of the original transfer function. While IRKA does not ensure stability (or even passivity) of the ROM a priori, it leads to (locally) $\cH_2$ optimal models upon convergence. IRKA-PH, on the other hand, preserves the pH structure and thus produces passive ROMs which, however, generally only fulfill a subset of the $\cH_2$ optimality conditions. Consequently, no $\cH_2$ optimality is achieved in general. The matrix $V$ can also be chosen in order to approximate the Dirac structure of the original model, resulting in the effort- and flow-constraint method \cite{Polyuga2012}. The extension of IRKA-PH and Dirac structure approximation to pH-DAEs is addressed in \cite{HauMM19,BeaGM21}. 

\subsection{Passivity-preserving MOR techniques}%
Positive real balanced truncation (PRBT) for ODEs is a well-studied MOR method, see  \cite{GugA04} and the references therein for a survey. The method is based on computing a minimal solution $X_{\min}$ of \eqref{eq:KYP_ineq} (in the sense of the Loewner ordering in symmetric matrices) and a minimal solution $Y_{\min}$ of the \emph{dual KYP inequality} 
\begin{equation*}
    \mat{-AY - YA^\T & B-YC^\T \\ B^\T-C Y & D+D^\T} \geq 0.
\end{equation*}
The solutions $X_{\min}$ and $Y_{\min}$ are then used to transform the system to a positive-real balanced realization where the transformed minimal solution of the KYP inequalities $\hat{X}_{\min}$ and $\hat{Y}_{\min}$ are equal and diagonal, i.\,e.,  $\hat{X}_{\min} = \hat{Y}_{\min} = \text{diag}(\eta_1,\ldots,\eta_n)$. Then the states corresponding to the small positive real characteristic values $\eta_i$ can be truncated. Finally, PRBT admits an a priori error bound in the gap metric that is derived in \cite{GuiO14}.

Another passivity-preserving MOR technique is  achieved via spectral factorization \cite{BreU21}. It initially requires a factorization 
\begin{equation}\label{eq:specfac}
	\mathcal{W}(X)=\mat{L&M}^\T\mat{L&M}
\end{equation} 
with $L\in\R^{k\times n}$, $M\in\R^{k\times m}$ and where $\text{rank}\mat{L&M}$ is as small as possible. If the pair $(A,B)$ is \emph{stabilizable} and $X = X_{\min}$ is the minimal solution of $\mathcal{W}(X) \ge 0$, then the resulting spectral factor 
system $(A,B,L,M)$
may be reduced via traditional (unstructured) MOR techniques, such as IRKA or balanced truncation (BT), and the passive ROM is obtained from the reduced spectral factor. 

For all passivity-preserving reduction methods, the passive ROM may eventually be transformed back to a pH representation by applying Theorem \ref{thm:KYP}.



\section{Introduction}

\file{elsarticle.cls} is a thoroughly re-written document class
for formatting \LaTeX{} submissions to Elsevier journals.
The class uses the environments and commands defined in \LaTeX{} kernel
without any change in the signature so that clashes with other
contributed \LaTeX{} packages such as \file{hyperref.sty},
\file{preview-latex.sty}, etc., will be minimal.
\file{elsarticle.cls} is primarily built upon the default
\file{article.cls}.  This class depends on the following packages
for its proper functioning:

\begin{enumerate}
\item \file{natbib.sty} for citation processing;
\item \file{geometry.sty} for margin settings;
\item \file{fleqn.clo} for left aligned equations;
\item \file{graphicx.sty} for graphics inclusion;
\item \file{txfonts.sty} optional font package, if the document is to
  be formatted with Times and compatible math fonts;
\item \file{hyperref.sty} optional packages if hyperlinking is
  required in the document;
\item \file{endfloat.sty} optional packages if floats to be placed at
 end of the PDF.
\end{enumerate}

All the above packages (except some optional packages) are part of any
standard \LaTeX{} installation. Therefore, the users need not be
bothered about downloading any extra packages.  Furthermore, users are
free to make use of \textsc{ams} math packages such as
\file{amsmath.sty}, \file{amsthm.sty}, \file{amssymb.sty},
\file{amsfonts.sty}, etc., if they want to.  All these packages work in
tandem with \file{elsarticle.cls} without any problems.

\section{Major Differences}

Following are the major differences between \file{elsarticle.cls}
and its predecessor package, \file{elsart.cls}:

\begin{enumerate}[\textbullet]
\item \file{elsarticle.cls} is built upon \file{article.cls}
while \file{elsart.cls} is not. \file{elsart.cls} redefines
many of the commands in the \LaTeX{} classes/kernel, which can
possibly cause surprising clashes with other contributed
\LaTeX{} packages;

\item provides preprint document formatting by default, and
optionally formats the document as per the final
style of models $1+$, $3+$ and $5+$ of Elsevier journals;

\item some easier ways for formatting \verb+list+ and
\verb+theorem+ environments are provided while people can still
use \file{amsthm.sty} package;

\item \file{natbib.sty} is the main citation processing package
  which can comprehensively handle all kinds of citations and
works perfectly with \file{hyperref.sty} in combination with
\file{hypernat.sty};

\item long title pages are processed correctly in preprint and
  final formats.

\end{enumerate}

\section{Installation}

The package is available at author resources page at Elsevier
(\url{http://www.elsevier.com/locate/latex}).
It can also be found in any of the nodes of the Comprehensive
\TeX{} Archive Network (\textsc{ctan}), one of the primary nodes
being
\url{http://tug.ctan.org/tex-archive/macros/latex/contrib/elsarticle/}.
Please download the \file{elsarticle.dtx} which is a composite
class with documentation and \file{elsarticle.ins} which is the
\LaTeX{} installer file. When we compile the
\file{elsarticle.ins} with \LaTeX{} it provides the class file,
\file{elsarticle.cls} by
stripping off all the documentation from the \verb+*.dtx+ file.
The class may be moved or copied to a place, usually,
\verb+$TEXMF/tex/latex/elsevier/+,
or a folder which will be read                   
by \LaTeX{} during document compilation.  The \TeX{} file
database needs updation after moving/copying class file.  Usually,
we use commands like \verb+mktexlsr+ or \verb+texhash+ depending
upon the distribution and operating system.


\section{Usage}\label{sec:usage}
The class should be loaded with the command:

\begin{vquote}
 \documentclass[<options>]{elsarticle}
\end{vquote}

\noindent where the \verb+options+ can be the following:


\begin{description}

\item [{\tt\color{verbcolor} preprint}]  default option which format the
  document for submission to Elsevier journals.

\item [{\tt\color{verbcolor} review}]  similar to the \verb+preprint+
option, but increases the baselineskip to facilitate easier review
process.

\item [{\tt\color{verbcolor} 1p}]  formats the article to the look and
feel of the final format of model 1+ journals. This is always single
column style.

\item [{\tt\color{verbcolor} 3p}] formats the article to the look and
feel of the final format of model 3+ journals. If the journal is a two
column model, use \verb+twocolumn+ option in combination.

\item [{\tt\color{verbcolor} 5p}] formats for model 5+ journals. This
is always of two column style.

\item [{\tt\color{verbcolor} authoryear}] author-year citation style of
\file{natbib.sty}. If you want to add extra options of
\file{natbib.sty}, you may use the options as comma delimited strings
as arguments to \verb+\biboptions+ command. An example would be:
\end{description}

\begin{vquote}
 \biboptions{longnamesfirst,angle,semicolon}
\end{vquote}

\begin{description}
\item [{\tt\color{verbcolor} number}] numbered citation style. Extra options
  can be loaded with\linebreak \verb+\biboptions+ command.

\item [{\tt\color{verbcolor} sort\&compress}] sorts and compresses the
numbered citations. For example, citation [1,2,3] will become [1--3].

\item [{\tt\color{verbcolor} longtitle}] if front matter is unusually long, use
  this option to split the title page across pages with the correct
placement of title and author footnotes in the first page.

\item [{\tt\color{verbcolor} times}] loads \file{txfonts.sty}, if
available in the system to use Times and compatible math fonts.

\item [{\tt\color{verbcolor} reversenotenum}] Use alphabets as
author--affiliation linking labels and use numbers for author
footnotes. By default, numbers will be used as author--affiliation
linking labels and alphabets for author footnotes. 

\item [{\tt\color{verbcolor} lefttitle}] To move title and
author/affiliation block to flushleft. \verb+centertitle+ is the
default option which produces center alignment.

\item [{\tt\color{verbcolor} endfloat}] To place all floats at the end
of the document.

\item [{\tt\color{verbcolor} nonatbib}] To unload natbib.sty.

\item [{\tt\color{verbcolor} doubleblind}] To hide author name, 
affiliation, email address etc. for double blind refereeing purpose.

\item[] All options of \file{article.cls} can be used with this
  document class.

\item[] The default options loaded are \verb+a4paper+, \verb+10pt+,
  \verb+oneside+, \verb+onecolumn+ and \verb+preprint+.

\end{description}

\section{Frontmatter}

There are two types of frontmatter coding:
\begin{enumerate}[(1)]
\item each author is
connected to an affiliation with a footnote marker; hence all
authors are grouped together and affiliations follow;
\pagebreak
\item authors of same affiliations are grouped together and the
relevant affiliation follows this group. 
\end{enumerate}

An example of coding the first type is provided below.

\begin{vquote}
 \title{This is a specimen title\tnoteref{t1,t2}}
 \tnotetext[t1]{This document is the results of the research
    project funded by the National Science Foundation.}
 \tnotetext[t2]{The second title footnote which is a longer 
    text matter to fill through the whole text width and 
    overflow into another line in the footnotes area of the 
    first page.}
\end{vquote}

\begin{vquote}
\author[1]{Jos Migchielsen\corref{cor1}%
  \fnref{fn1}}
\ead{J.Migchielsen@elsevier.com}

\author[2]{CV Radhakrishnan\fnref{fn2}}
\ead{cvr@sayahna.org}

\author[3]{CV Rajagopal\fnref{fn1,fn3}}
\ead[url]{www.stmdocs.in}

 \cortext[cor1]{Corresponding author}
 \fntext[fn1]{This is the first author footnote.}
 \fntext[fn2]{Another author footnote, this is a very long 
   footnote and it should be a really long footnote. But this 
   footnote is not yet sufficiently long enough to make two 
   lines of footnote text.}
 \fntext[fn3]{Yet another author footnote.}

 \affiliation[1]{organization={Elsevier B.V.}, 
                 addressline={Radarweg 29},
                 postcode={1043 NX}, 
                 city={Amsterdam}, 
                 country={The Netherlands}}

 \affiliation[2]{organization={Sayahna Foundation},
                 addressline={JWRA 34, Jagathy}, 
                 city={Trivandrum}
                 postcode={695014}, 
                 country={India}}
\end{vquote}

\begin{vquote}

 \affiliation[3]{organization={STM Document Engineering 
                               Pvt Ltd.},
                 addressline={Mepukada, Malayinkil}, 
                 city={Trivandrum}
                 postcode={695571}, 
                 country={India}}
\end{vquote}

The output of the above \TeX{} source is given in Clips~\ref{clip1} and
\ref{clip2}. The header portion or title area is given in
Clip~\ref{clip1} and the footer area is given in Clip~\ref{clip2}.

\deforange{blue!70}
\src{Header of the title page.}
\includeclip{1}{130 612 477 707}{1psingleauthorgroup.pdf
\deforange{orange}

\deforange{blue!70}
\src{Footer of the title page.}
\includeclip{1}{93 135 499 255}{1pseperateaug.pdf
\deforange{orange}

Most of the commands such as \verb+\title+, \verb+\author+,
\verb+\affiliation+ are self explanatory.  Various components are
linked to each other by a label--reference mechanism; for
instance, title footnote is linked to the title with a footnote
mark generated by referring to the \verb+\label+ string of
the \verb=\tnotetext=.  We have used similar commands
such as \verb=\tnoteref= (to link title note to title);
\verb=\corref= (to link corresponding author text to
corresponding author); \verb=\fnref= (to link footnote text to
the relevant author names).  \TeX{} needs two compilations to
resolve the footnote marks in the preamble part.  
Given below are the syntax of various note marks and note texts.


\begin{vquote}
  \tnoteref{<label(s)>}
  \corref{<label(s)>}
  \fnref{<label(s)>}
  \tnotetext[<label>]{<title note text>}
  \cortext[<label>]{<corresponding author note text>}
  \fntext[<label>]{<author footnote text>}
\end{vquote}

\noindent where \verb=<label(s)>= can be either one or more comma
delimited label strings. The optional arguments to the
\verb=\author= command holds the ref label(s) of the address(es)
to which the author is affiliated while each \verb=\affiliation=
command can have an optional argument of a label. In the same
manner, \verb=\tnotetext=, \verb=\fntext=, \verb=\cortext= will
have optional arguments as their respective labels and note text
as their mandatory argument.

The following example code provides the markup of the second type
of author-affiliation.

\begin{vquote}
\author{Jos Migchielsen\corref{cor1}%
  \fnref{fn1}}
\ead{J.Migchielsen@elsevier.com}
 \affiliation[1]{organization={Elsevier B.V.}, 
                 addressline={Radarweg 29},
                 postcode={1043 NX}, 
                 city={Amsterdam}, 
                 country={The Netherlands}}

\author{CV Radhakrishnan\fnref{fn2}}
\ead{cvr@sayahna.org}
 \affiliation[2]{organization={Sayahna Foundation},
                 addressline={JWRA 34, Jagathy}, 
                 city={Trivandrum}
                 postcode={695014}, 
                 country={India}}

\author{CV Rajagopal\fnref{fn1,fn3}}
\ead[url]{www.stmdocs.in}
 \affiliation[3]{organization={STM Document Engineering 
                               Pvt Ltd.},
                 addressline={Mepukada, Malayinkil}, 
                 city={Trivandrum}
                 postcode={695571}, 
                 country={India}}
\end{vquote}

\vspace*{-.5pc}

\begin{vquote}
\cortext[cor1]{Corresponding author}
\fntext[fn1]{This is the first author footnote.}
\fntext[fn2]{Another author footnote, this is a very long 
  footnote and it should be a really long footnote. But this 
  footnote is not yet sufficiently long enough to make two lines 
  of footnote text.}
\end{vquote}

The output of the above \TeX{} source is given in Clip~\ref{clip3}.

\deforange{blue!70}
\src{Header of the title page..}
\includeclip{1}{119 563 468 709}{1pseperateaug.pdf
\deforange{orange}

Clip~\ref{clip4} shows the output after giving \verb+doubleblind+ class option. 

\deforange{blue!70}
\src{Double blind article}
\includeclip{1}{124 567 477 670}{elstest-1pdoubleblind.pdf
\deforange{orange}

\vspace*{-.5pc}
The frontmatter part has further environments such as abstracts and
keywords.  These can be marked up in the following manner:

\begin{vquote}
 \begin{abstract}
  In this work we demonstrate the formation of a new type of 
  polariton on the interface between a ....
 \end{abstract}
\end{vquote} 

\vspace*{-.5pc}
\begin{vquote}
 \begin{keyword}
  quadruple exiton \sep polariton \sep WGM
 \end{keyword}
\end{vquote}

\noindent Each keyword shall be separated by a \verb+\sep+ command.
\textsc{msc} classifications shall be provided in 
the keyword environment with the commands
\verb+\MSC+. \verb+\MSC+ accepts an optional
argument to accommodate future revisions.
eg., \verb=\MSC[2008]=. The default is 2000.\looseness=-1

\subsection{New page}
Sometimes you may need to give a page-break and start a new page after
title, author or abstract. Following commands can be used for this
purpose.

\begin{vquote}
  \newpageafter{title}
  \newpageafter{author}
  \newpageafter{abstract}
\end{vquote}


\begin{itemize}
\leftskip-2pc
\item [] {\tt\color{verbcolor} \verb+\newpageafter{title}+} typeset the title alone on one page.

\item [] {\tt\color{verbcolor} \verb+\newpageafter{author}+}  typeset the title
and author details on one page.

\item [] {\tt\color{verbcolor} \verb+\newpageafter{abstract}+}
typeset the title,
author details and abstract \& keywords one one page.

\end{itemize}

\section{Floats}
{Figures} may be included using the command, \verb+\includegraphics+ in
combination with or without its several options to further control
graphic. \verb+\includegraphics+ is provided by \file{graphic[s,x].sty}
which is part of any standard \LaTeX{} distribution.
\file{graphicx.sty} is loaded by default. \LaTeX{} accepts figures in
the postscript format while pdf\LaTeX{} accepts \file{*.pdf},
\file{*.mps} (metapost), \file{*.jpg} and \file{*.png} formats. 
pdf\LaTeX{} does not accept graphic files in the postscript format. 

The \verb+table+ environment is handy for marking up tabular
material. If users want to use \file{multirow.sty},
\file{array.sty}, etc., to fine control/enhance the tables, they
are welcome to load any package of their choice and
\file{elsarticle.cls} will work in combination with all loaded
packages.

\section[Theorem and ...]{Theorem and theorem like environments}

\file{elsarticle.cls} provides a few shortcuts to format theorems and
theorem-like environments with ease. In all commands the options that
are used with the \verb+\newtheorem+ command will work exactly in the same
manner. \file{elsarticle.cls} provides three commands to format theorem or
theorem-like environments: 

\begin{vquote}
 \newtheorem{thm}{Theorem}
 \newtheorem{lem}[thm]{Lemma}
 \newdefinition{rmk}{Remark}
 \newproof{pf}{Proof}
 \newproof{pot}{Proof of Theorem \ref{thm2}}
\end{vquote}

The \verb+\newtheorem+ command formats a
theorem in \LaTeX's default style with italicized font, bold font
for theorem heading and theorem number at the right hand side of the
theorem heading.  It also optionally accepts an argument which
will be printed as an extra heading in parentheses. 

\begin{vquote}
  \begin{thm} 
   For system (8), consensus can be achieved with 
   $\|T_{\omega z}$
   ...
     \begin{eqnarray}\label{10}
     ....
     \end{eqnarray}
  \end{thm}
\end{vquote}  

Clip~\ref{clip5} will show you how some text enclosed between the
above code\goodbreak \noindent looks like:

\vspace*{6pt}
\deforange{blue!70}
\src{{\ttfamily\color{verbcolor}\expandafter\@gobble\string\\ newtheorem}}
\includeclip{2}{1 1 453 120}{jfigs.pdf}
\deforange{orange}

The \verb+\newdefinition+ command is the same in
all respects as its\linebreak \verb+\newtheorem+ counterpart except that
the font shape is roman instead of italic.  Both
\verb+\newdefinition+ and \verb+\newtheorem+ commands
automatically define counters for the environments defined.

\vspace*{6pt}
\deforange{blue!70}
\src{{\ttfamily\color{verbcolor}\expandafter\@gobble\string\\ newdefinition}}
\includeclip{1}{1 1 453 105}{jfigs.pdf}
\deforange{orange}

The \verb+\newproof+ command defines proof environments with
upright font shape.  No counters are defined. 

\vspace*{6pt}
\deforange{blue!70}
\src{{\ttfamily\color{verbcolor}\expandafter\@gobble\string\\ newproof}}
\includeclip{3}{1 1 453 65}{jfigs.pdf}
\deforange{orange}

Users can also make use of \verb+amsthm.sty+ which will override
all the default definitions described above.

\section[Enumerated ...]{Enumerated and Itemized Lists}
\file{elsarticle.cls} provides an extended list processing macros
which makes the usage a bit more user friendly than the default
\LaTeX{} list macros.   With an optional argument to the
\verb+\begin{enumerate}+ command, you can change the list counter
type and its attributes.

\begin{vquote}
 \begin{enumerate}[1.]
 \item The enumerate environment starts with an optional
   argument `1.', so that the item counter will be suffixed
   by a period.
 \item You can use `a)' for alphabetical counter and '(i)' for
   roman counter.
  \begin{enumerate}[a)]
    \item Another level of list with alphabetical counter.
    \item One more item before we start another.
\end{vquote}

\deforange{blue!70}
\src{List -- Enumerate}
\includeclip{4}{1 1 453 185}{jfigs.pdf}
\deforange{orange}

Further, the enhanced list environment allows one to prefix a
string like `step' to all the item numbers.  

\begin{vquote}
 \begin{enumerate}[Step 1.]
  \item This is the first step of the example list.
  \item Obviously this is the second step.
  \item The final step to wind up this example.
 \end{enumerate}
\end{vquote}

\deforange{blue!70}
\src{List -- enhanced}
\includeclip{5}{1 1 313 83}{jfigs.pdf}
\deforange{orange}


\section{Cross-references}
In electronic publications, articles may be internally
hyperlinked. Hyperlinks are generated from proper
cross-references in the article.  For example, the words
\textcolor{black!80}{Fig.~1} will never be more than simple text,
whereas the proper cross-reference \verb+\ref{tiger}+ may be
turned into a hyperlink to the figure itself:
\textcolor{blue}{Fig.~1}.  In the same way,
the words \textcolor{blue}{Ref.~[1]} will fail to turn into a
hyperlink; the proper cross-reference is \verb+\cite{Knuth96}+.
Cross-referencing is possible in \LaTeX{} for sections,
subsections, formulae, figures, tables, and literature
references.

\section[Mathematical ...]{Mathematical symbols and formulae}

Many physical/mathematical sciences authors require more
mathematical symbols than the few that are provided in standard
\LaTeX. A useful package for additional symbols is the
\file{amssymb} package, developed by the American Mathematical
Society. This package includes such oft-used symbols as
$\lesssim$ (\verb+\lesssim+), $\gtrsim$ (\verb+\gtrsim+)  or 
$\hbar$ (\verb+\hbar+). Note that your \TeX{}
system should have the \file{msam} and \file{msbm} fonts installed. If
you need only a few symbols, such as $\Box$ (\verb+\Box+), you might try the
package \file{latexsym}.

Another point which would require authors' attention is the
breaking up of long equations.  When you use
\file{elsarticle.cls} for formatting your submissions in the 
\verb+preprint+ mode, the document is formatted in single column
style with a text width of 384pt or 5.3in.  When this document is
formatted for final print and if the journal happens to be a double column
journal, the text width will be reduced to 224pt at for 3+
double column and 5+ journals respectively. All the nifty 
fine-tuning in equation breaking done by the author goes to waste in
such cases.  Therefore, authors are requested to check this
problem by typesetting their submissions in final format as well
just to see if their equations are broken at appropriate places,
by changing appropriate options in the document class loading
command, which is explained in section~\ref{sec:usage},
\nameref{sec:usage}. This allows authors to fix any equation breaking
problem before submission for publication.
\file{elsarticle.cls} supports formatting the author submission
in different types of final format.  This is further discussed in
section \ref{sec:final}, \nameref{sec:final}.

\enlargethispage*{\baselineskip}

\subsection*{Displayed equations and double column journals}

Many Elsevier journals print their text in two columns. Since
the preprint layout uses a larger line width than such columns,
the formulae are too wide for the line width in print. Here is an
example of an equation  (see equation 6) which is perfect in a
single column preprint format:

In normal course, articles are prepared and submitted in single column
format even if the final printed article will come in a double column
format journal. Here the problem is that when the article is typeset by
the typesetters for paginating and fit within the single column width,
they have to break the lengthy equations and align them properly. Even
if most of the tasks in preparing your proof is automated, the equation
breaking and aligning requires manual judgement, hence this task is manual.
When there comes a manual operation that area is error prone. Author
needs to check that equation pretty well. 

However if authors themselves break the equation to the single column
width typesetters need not want to touch these area and the proof authors
get will be without any errors.

\setlength\Sep{6pt}
\src{See equation (6)}
\deforange{blue!70}
\includeclip{4}{105 500 500 700}{1psingleauthorgroup.pdf}
\deforange{orange}
                 	
\noindent When this document is typeset for publication in a
model 3+ journal with double columns, the equation will overlap
the second column text matter if the equation is not broken at
the appropriate location.

\vspace*{6pt}
\deforange{blue!70}
\src{See equation (6) overprints into second column}
\includeclip{3}{59 421 532 635}{elstest-3pd.pdf}
\deforange{orange}
\vspace*{6pt}

\noindent The typesetter will try to break the equation which
need not necessarily be to the liking of the author or as it
happens, typesetter's break point may be semantically incorrect.
Therefore, authors may check their submissions for the incidence
of such long equations and break the equations at the correct
places so that the final typeset copy will be as they wish.

\section{Bibliography}

Three bibliographic style files (\verb+*.bst+) are provided ---
\file{elsarticle-num.bst}, \file{elsarticle-num-names.bst} and
\file{elsarticle-harv.bst} --- the first one can be used for the
numbered scheme, second one for numbered with new options of 
\file{natbib.sty}. The third one is for the author year
scheme.

In \LaTeX{} literature, references are listed in the
\verb+thebibliography+ environment.  Each reference is a
\verb+\bibitem+ and each \verb+\bibitem+ is identified by a label,
by which it can be cited in the text:

\verb+\bibitem[Elson et al.(1996)]{ESG96}+ is cited as
\verb+\citet{ESG96}+. 

\noindent In connection with cross-referencing and
possible future hyperlinking it is not a good idea to collect
more that one literature item in one \verb+\bibitem+.  The
so-called Harvard or author-year style of referencing is enabled
by the \LaTeX{} package \file{natbib}. With this package the
literature can be cited as follows:

\begin{enumerate}[\textbullet]
\item Parenthetical: \verb+\citep{WB96}+ produces (Wettig \& Brown, 1996).
\item Textual: \verb+\citet{ESG96}+ produces Elson et al. (1996).
\item An affix and part of a reference:
\verb+\citep[e.g.][Ch. 2]{Gea97}+ produces (e.g. Governato et
al., 1997, Ch. 2).
\end{enumerate}

In the numbered scheme of citation, \verb+\cite{<label>}+ is used,
since \verb+\citep+ or \verb+\citet+ has no relevance in the numbered
scheme.  \file{natbib} package is loaded by \file{elsarticle} with
\verb+numbers+ as default option.  You can change this to author-year
or harvard scheme by adding option \verb+authoryear+ in the class
loading command.  If you want to use more options of the \file{natbib}
package, you can do so with the \verb+\biboptions+ command, which is
described in the section \ref{sec:usage}, \nameref{sec:usage}.  For
details of various options of the \file{natbib} package, please take a
look at the \file{natbib} documentation, which is part of any standard
\LaTeX{} installation.

In addition to the above standard \verb+.bst+ files, there are 10
journal-specific \verb+.bst+ files also available.
Instruction for using these \verb+.bst+ files can be found at 
\href{http://support.stmdocs.in/wiki/index.php?title=Model-wise_bibliographic_style_files}
{http://support.stmdocs.in}

\section[Graphical ...]{Graphical abstract and highlights}
A template for adding graphical abstract and highlights are available
now. This will appear as the first two pages of the PDF before the
article content begins.

\pagebreak
Please refer below to see how to code them.

\begin{vquote}
....
....

\end{abstract}

\begin{graphicalabstract}
\end{graphicalabstract}

\begin{highlights}
\item Research highlight 1
\item Research highlight 2
\end{highlights}

\begin{keyword}
....
....
\end{vquote}

\section{Final print}\label{sec:final}

The authors can format their submission to the page size and margins
of their preferred journal.  \file{elsarticle} provides four
class options for the same. But it does not mean that using these
options you can emulate the exact page layout of the final print copy. 


\lmrgn=3em
\begin{description}
\item [\texttt{1p}:] $1+$ journals with a text area of
384pt $\times$ 562pt or 13.5cm $\times$ 19.75cm or 5.3in $\times$
7.78in, single column style only.

\item [\texttt{3p}:] $3+$ journals with a text area of 468pt
$\times$ 622pt or 16.45cm $\times$ 21.9cm or 6.5in $\times$
8.6in, single column style.

\item [\texttt{twocolumn}:] should be used along with 3p option if the
journal is $3+$ with the same text area as above, but double column
style. 

\item [\texttt{5p}:] $5+$ with text area of 522pt $\times$
682pt or 18.35cm $\times$ 24cm or 7.22in $\times$ 9.45in,
double column style only.
\end{description}

Following pages have the clippings of different parts of
the title page of different journal models typeset in final
format.

Model $1+$ and $3+$  will have the same look and
feel in the typeset copy when presented in this document. That is
also the case with the double column $3+$ and $5+$ journal article
pages. The only difference will be wider text width of
higher models.  Here are the specimen single and double column journal 
pages.



\begin{comment}
\begin{center}
\hypertarget{bsc}{}
\hyperlink{sc}{
{\bf [Specimen single column article -- Click here]}
}


\hypertarget{bsc}{}
\hyperlink{dc}{
{\bf [Specimen double column article -- Click here]}
}
\end{center}
\end{comment}

\vspace*{-.5pc}

\enlargethispage*{\baselineskip}

\src{}\hypertarget{sc}{}
\deforange{blue!70}
\hyperlink{bsc}{\includeclip{1}{88 120 514 724}{elstest-1p.pdf}}
\deforange{orange}

\src{}\hypertarget{dc}{}
\deforange{blue!70}
\hyperlink{bsc}{\includeclip{1}{27 61 562 758}{elstest-5p.pdf}}
\deforange{orange}

~\hfill $\Box$

\end{document}


\section{}
\label{}






\section{}
\label{}







\section{}
\label{}






\section{Introduction}%
\label{sec:introduction}
\input{Introduction}

\section{Preliminaries}%
\label{sec:preliminaries}
\input{Prelims}

\section{A new optimization-based approach}%
\label{sec:our_method}

\input{pauls_approach}
\input{tims_approach}

\section{Numerical examples}%
\label{sec:numerical_examples}

\input{Experimental_setup.tex}

\input{Results.tex}

\section{Conclusion}%
\label{sec:conclusion}

We have presented two optimization-based methods for structure-preserving MOR of pH-DAEs. These make it possible to compute accurate ROMs with respect to either the $\hinf$ or the $\htwo$ norm. The main benefits compared to state-of-the-art methods are the simplified treatment of the algebraic equations, which can be incorporated into the parameterized ROM in a structure-preserving way without increasing its state dimension. Furthermore, our methods are data-driven, such that no transformations to the FOM system matrices are required. Nonetheless, we have shown, how transformations can be applied in order to obtain an accurate estimate of the feedthrough, which is essential in the $\htwo$ case. Finally, our numerical experiments show that the optimization-based methods often lead to a higher accuracy (especially in the $\hinf$ norm).

We are currently investigating the application of our method to higher index pH-DAEs. These may have improper transfer functions, which are not currently supported in our parameterization.

\bibliographystyle{elsarticle-num} 


\subsection{$\hinf$ approximation}%
\label{sub:ourmethod_hinf}

The algorithm presented in \cite{Schwerdtner2020} to obtain a good $\hinf$ approximation of a pH-FOM with transfer function $\pHtf$ using a parametrized low-order system with transfer function $\pHtfr( \cdot, \theta)$ is based on minimizing the objective function
\begin{align}
  \begin{split}
  \pushleft{\loss(\theta;\pHtf,\gamma,\mathcal{S}) :=} \\ \pushright{\frac{1}{\gamma}\sum\limits_{s_i\in \mathcal{S}}
  \sum\limits_{j=1}^{m}{\left({\left[\sigma_j \left(\pHtf(s_i)-\pHtfr(s_i,\theta)\right)-\gamma\right]}_+\right)}^2}
\end{split}
  \label{eq:loss}
\end{align}
with respect to $\theta$, where 
\begin{align*}
  {[ \cdot ]}_+:  \R \rightarrow [0,\infty), \quad x \mapsto 
  \begin{cases}
    x & \text{if } x\ge 0,\\
    0 & \text{if } x<0,
  \end{cases}
\end{align*}
for decreasing values of $\gamma > 0$. Here $\mathcal{S} \subset \ri \R$ is a set of sample points, at which the original and reduced transfer functions are evaluated and $\sigma_j(\cdot)$ denotes the $j$-th singular value of its matrix argument.

The justification for using 
$\loss$ 
as a surrogate for the $\hinf$ error is that $\loss(\cdot;H,\gamma,\mathcal{S})$ attains its global minimum (at zero), when all singular values of the error transfer function at all sample points $s_i$ are below the threshold~$\gamma$. Therefore, if the sample points are chosen appropriately, then $\loss(\theta;\pHtf,\gamma,\mathcal{S}) = 0$ is an indication for ${{\|\pHtf - \pHtf_r(\cdot,\theta)\|}_{\cH_\infty} \le \gamma}$. In \cite{Schwerdtner2021}, an adaptive sampling procedure is introduced which ensures that the sample points are appropriately distributed along the imaginary axis based on the given FOM and status of the optimization.

In this work, we use Algorithm~\ref{alg:bisection} to determine an $\hinf$ approximation via a bisection procedure, which determines the minimal value for $\gamma$ (up to a relative tolerance $\epsilon_1$) at which a minimization of $\loss$ with respect to $\theta$ terminates at zero. The tolerance $\epsilon_2$ that is used in line $6$ of the algorithm is the maximum value of $\loss$, which is still numerically interpreted as zero, such that $\gamma$ is reduced in the subsequent bisection step. The sample points are updated after each bisection step because the adaptive sampling update rule introduced in~\cite{Schwerdtner2021} depends on the current value of $\gamma$.

\begin{algorithm}[tbh]
  \LinesNumbered
  \SetAlgoLined
  \DontPrintSemicolon
  \SetKwInOut{Input}{Input}\SetKwInOut{Output}{Output}
  \Input{FOM transfer function $\pHtf\in\mathcal{RH}_{\infty}^{\enu \times \enu}$,
  initial ROM transfer function $\pHtfr( \cdot, \theta_0) \in \rhinf^{\enu \times \enu}$ with parameter $\theta_0 \in \R^{n_\theta}$,
  initial sample point set $\mathcal{S} \subset \mathrm{i}\R$,
  upper bound $\gamma_{\rm u} > 0$,
  bisection tolerance $\varepsilon_1 > 0$,
  termination tolerance $\varepsilon_2 > 0$
}
  \Output{Reduced pH-ODE of order $r$}
  Set $j:=0$ and $\gamma_{\rm l}:=0$.\;
  \While{$(\gamma_{\rm u}-\gamma_{\rm l})/(\gamma_{\rm u}+\gamma_{\rm l}) > \varepsilon_1$}{
    Set $\gamma:=(\gamma_{\rm u}+\gamma_{\rm l})/2$.\;
    Update the sample set $\mathcal{S}$ using~\cite[Alg.~3.1]{Schwerdtner2021}. \;
    Solve the minimization problem $\alpha := \min_{\theta\in \R^{n_\theta}}\loss(\theta;\pHtf,\gamma,\mathcal{S})$ with minimizer $\theta_{j+1} \in \R^{n_\theta}$, initialized at $\theta_j$. \;
  \eIf{$\alpha > \varepsilon_2$}{
    Set $\gamma_{\rm l}:=\gamma$.\;
    }{
    Set $\gamma_{\rm u}:=\gamma$.\;
  }
  Set $j:=j+1$.
  }
  Construct the ROM with $\theta_j$ as in Lemma~\ref{lem:PHParam}.
  \caption{SOBMOR-$\hinf$}
  \label{alg:bisection}
\end{algorithm}

The benefits of using this approach instead of directly minimizing the $\hinf$ norm were discussed in detail in~\cite[Remark~3.3]{Schwerdtner2020}. The main reasons for using~\eqref{eq:loss} instead of the $\hinf$ norm are the differentiability of~$\loss$ with respect to~$\theta$, the local convergence of the method, and the prohibitive computational costs as well as reliability issues of the $\hinf$ norm computation (for the large-scale error system) inside an optimization loop.

\subsection{$\htwo$ approximation}%
\label{sub:ourmethod_htwo}
In order to obtain a finite $\htwo$ error, the polynomial parts of the FOM and the ROM transfer function must be equal. When using projection-based methods on ODE models this \emph{feedthrough matching} is automatic. In the DAE case, this is typically obtained by preserving the algebraic part, i.\,e., by including the null-space of the $E$-matrix in the projection matrices. For systems with multiple algebraic constraints this is undesirable, and a reduction of the subsystem corresponding to the  algebraic constraints may be necessary, see, for instance, \cite{MehS05}.

Another remedy (used for pH-DAEs with index one in \cite{BeaGM21}) is to compute the polynomial part $\pHtf_{\rm p}(s) \equiv \mathcal{D}_0$ of the FOM before the reduction and include it in the feedthrough terms $\rS$ and $\rN$ of the ROM, which we will use here as well. The direct computation of $\mathcal{D}_0$ (see Section \ref{sec:preliminaries}) may, however, require transformations of the FOM. Alternatively, $\mathcal{D}_0$ can also be estimated by sampling the transfer function of the FOM at sufficiently large $s\in\C$ in an iterative manner as proposed in \cite{Schwerdtner2020a}. We then decompose it in its symmetric and skew-symmetric part, respectively, i.\,e., 
\begin{align}
  S_0 &:= \frac{1}{2}\left(\mathcal{D}_0^\T+\mathcal{D}_0\right), \label{eq:S0} \\ 
  N_0 &:= \frac{1}{2}\left(\mathcal{D}_0^\T-\mathcal{D}_0\right). \label{eq:N0} 
\end{align}
If the ROM is parameterized as in Lemma~\ref{lem:PHParam} then the $\cH_2$ error ${\Vert \pHtf - \pHtfr(\cdot,\theta)\Vert}_{\cH_2}$ is only well-defined if we have that 
\begin{align}\label{eq:SN_match}
	\rS(\theta)  &= S_0,\quad \rN(\theta) = N_0,
\end{align}
since otherwise $\pHtf-\pHtfr(\cdot,\theta) \notin \mathcal{RH}_{2}^{m \times m}$. Consequently,  we first have to fix all parameters in  $\theta_N$ such that ${\rN(\theta) = N_0}$. Since we indirectly parameterize $\rS(\theta)$ via $\theta_W$, we first analyze which parameters in $\theta_W$ have an impact on $\rS(\theta)$. For this, consider a separation of ${\theta_W \in \R^{n_W}}$ into ${\theta_W =: \mat{\theta_{W_1}^{\T}, \theta_{W_2}^{\T}}^\T}$, where ${\theta_{W_1} \in \R^{n_W - m(m+1)/2}}$ and ${\theta_{W_2} \in \R^{m(m+1)/2}}$. Then we can decompose 
\begin{align}
	\vtu(\theta_W) = \mat{\Xi_1 & \Xi_2 \\ 0 & \Xi_3},
\end{align}
where the matrices $\Xi_1, \,\Xi_2$ depend only on $\theta_{W_1}$, and $\Xi_3$ depends only on~$\theta_{W_2}$. Consequently, ${\rS(\theta) = \Xi_3\Xi_3^\T}$ only depends on~$\theta_{W_2}$ and we can set $\theta_{W_2}$ such that ${\rS(\theta) = S_0}$. The remaining parameters $\theta_{W_1}$ may still be subject to optimization and it holds that ${W(\theta) \geq 0}$ for all ${\theta_{W_1} \in \R^{n_W - m(m+1)/2}}$. Consequently, for minimizing the $\cH_2$ error, the parameter vector which is subject to optimization reduces to $\theta := \begin{bmatrix}
	\theta_J^\T, \theta_{W_1}^\T, \theta_Q^\T, \theta_{\pG}^\T
\end{bmatrix}^\T.$
\begin{rem}
  Note that $S_0 \geq 0$ always holds since the implicit pH-ODE~\eqref{eq:FOM_decpl} has the same transfer function $\pHtf$ as the original pH-DAE and is therefore positive real, i.\,e., we have that
  \begin{equation*}
    \Phi(\mathsf{i}\omega) := \pHtf(\mathsf{i}\omega) + \pHtf(\mathsf{i}\omega)^{\mathsf{H}} \geq 0, 
  \end{equation*}
  for all $\omega \in \R$ and consequently, $\lim_{\omega\rightarrow\infty}\Phi(\mathsf{i}\omega) = 2S_0\geq 0$.
\end{rem} 

Now we can formulate the $\cH_2$ optimization problem in the pole-residue framework originally proposed in \cite{Beattie2009a} for unstructured LTI systems and extended to pH-ODE systems in \cite{Moser2020}. Assume that
$(\rJ(\theta)-\rR(\theta))\rQ(\theta)$ is diagonalizable and consider the spectral decomposition
\begin{equation} \label{eq:red_evp}
	(\rJ(\theta)-\rR(\theta))\rQ(\theta) Z(\theta) = Z(\theta)\Lambda(\theta),
\end{equation}
where $\Lambda(\theta) = \text{diag}(\reig_1(\theta),\,...\,,\reig_r(\theta))$ and $Z(\theta)$ contains the right eigenvectors as columns. If the eigenvalues $\reig_i(\theta)~\in~\C$, $i=1,\,\ldots,\,r$ are simple, then the transfer function $\pHtfr(\cdot,\theta)$ may be represented by the partial fraction expansion ${\pHtfr(s,\theta) = \sum_{i=1}^{r} \frac{\rl_i(\theta)\rr_i(\theta)^\T}{s-\reig_i(\theta)}+\rS(\theta)-\rN(\theta)}$, where ${\rl_i(\theta),\,\rr_i(\theta)\in \C^m}$ with
\begin{align*}
	\rl_i(\theta) &= (\rG(\theta)+\rP(\theta))^\T\rQ(\theta) Z(\theta) e_i, \\
	\rr_i(\theta) &= (\rG(\theta)-\rP(\theta))^\T Z(\theta)^{-\T}e_i,
\end{align*}
and where $e_i$ denotes the $i$-th standard basis vector of $\R^r$. Assuming that $\eqref{eq:SN_match}$ holds, we have that	
\begin{equation}
	\begin{aligned}
		{\Vert \pHtf-\pHtfr(\cdot,\theta) \Vert}_{\cH_2}^2 = &{\Vert \pHtf_{\rm sp} \Vert}_{\cH_2}^2 \\	
				& - 2\sum_{i=1}^{r}\rl_i(\theta)^\T \pHtf_{\rm sp}(-\reig_i(\theta))\rr_i(\theta) \\ 
					& + \sum_{j,k=1}^{r} \frac{\rl_j(\theta)^\T\rl_k(\theta)\rr_k(\theta)^\T\rr_j(\theta)}{-\reig_j(\theta)-\reig_k(\theta)},
	\end{aligned}
\end{equation}
as shown in \cite[Theorem 2.1]{Beattie2009a}. Since ${\Vert \pHtf_{\rm sp} \Vert}_{\cH_2}^2$ does not depend on $\theta$, it can be neglected in the optimization. Consequently, we define the objective functional
\begin{align*}
\losst(\theta;\pHtf) &:={\Vert \pHtf-\pHtfr(\cdot,\theta) \Vert}_{\cH_2}^2 - {\Vert \pHtf_{\rm sp} \Vert}_{\cH_2}^2 \\
&:= \hat{\mathcal{F}}\Big(\big[\rl_1(\theta)^\T,\ldots ,\rl_r(\theta)^\T,\rr_1(\theta)^\T,\ldots\\ & \quad\quad\quad\quad \quad\rr_r(\theta)^\T,\reig_1(\theta)^\T,\ldots,\reig_r(\theta)^\T\big]^\T\Big) \\
&= (\hat{\mathcal{F}} \circ q)(\theta),
\end{align*}
where
\begin{align*}
	q(\theta) &:= [\rl_1(\theta)^\T,\ldots,\rl_r(\theta)^\T,\rr_1(\theta)^\T,\ldots \\ & \quad\quad\quad\quad\rr_r(\theta)^\T,\reig_1(\theta),\ldots,\reig_r(\theta)]^\T \in \C^{n_q}.
\end{align*}
This functional can be evaluated efficiently because it only requires the solution of the reduced-order eigenvalue problem in \eqref{eq:red_evp} as well as $r$ evaluations of $\pHtf_{\rm sp}$ at $-\reig_i(\theta)$. 
The eigenvalues $\reig_i(\theta)$ and rank-one residues $\rl_i(\theta)\rr_i(\theta)^\T$ are functions of the parameter vector $\theta$.
If $\bar{\theta} \in \R^{n_\theta}$ is chosen such that all eigenvalues are simple, then $\mathcal{F}$ is differentiable in a neighborhood of $\bar{\theta}$. 
Its derivative is obtained by applying the chain rule, i.\,e., with the differentiation operator $\text{D}$ we obtain
\begin{equation*}
 \text{D}\losst(\bar{\theta}) = \left(\nabla \losst(\bar{\theta})\right)^\T = \text{D}\hat{\mathcal{F}}(q(\bar{\theta}))\cdot \text{D}q(\bar{\theta}),
\end{equation*}
with
\begin{multline*}
	\text{D}\hat{\losst}(q(\bar{\theta})) = \left[\text{D}_{b_1}\hat{\losst}(q(\bar{\theta})),\ldots,\text{D}_{b_r}\hat{\losst}(q(\bar{\theta})), \ldots \right. \\ \left. \text{D}_{c_1}\hat{\losst}(q(\bar{\theta})),\ldots, \text{D}_{c_r}\hat{\losst}(q(\bar{\theta})), \ldots\right. \\ \left.\text{D}_{\lambda_1}\hat{\losst}(q(\bar{\theta})),\ldots,\text{D}_{\lambda_r}\hat{\losst}(q(\bar{\theta}))\right]\in\C^{1\times n_q},
\end{multline*}
and 
\begin{align*}
	\text{D}q(\bar{\theta}) &= \left[\text{D}_{\theta_1}{q(\bar{\theta})},\ldots,\text{D}_{\theta_{n_\theta}}{q(\bar{\theta})}\right]\in\C^{n_q\times n_{\theta}}.
\end{align*}
For all $i=1,\,\ldots,\,r$ it holds that
\begin{align*}
    \text{D}_{\rr_i} \hat{\losst}(q(\bar{\theta})) &= 2\rl_i(\bar{\theta})^\T\big(\pHtfr(-\reig_i(\bar{\theta}))-\pHtf(-\reig_i(\bar{\theta}))\big), \\
    \text{D}_{\rl_i} \hat{\losst}(q(\bar{\theta})) &= 2\rr_i(\bar{\theta})^\T\big(\pHtfr(-\reig_i(\bar{\theta}))-\pHtf(-\reig_i(\bar{\theta}))\big)^\T, \\
	\text{D}_{\lambda_i} \hat{\losst}(q(\bar{\theta})) &= -2\rl_i(\bar{\theta})^\T\big(\pHtfr^{\prime}(-\reig_i(\bar{\theta}))-\pHtf^{\prime}(-\reig_i(\bar{\theta}))\big)\rr_i(\bar{\theta}),
\end{align*}
and we refer  to~\cite{Moser2020} for the differentiation of $\text{D} q$. 
\begin{rem}
	The partial derivatives in $\mathrm{D}q(\bar{\theta})$ may be computed efficiently with block-wise expressions. For instance, the derivative $\mathrm{D}_{\theta_G}\rl_i(\bar{\theta}) \in \C^{m \times r\cdot m}$ can be computed as
	\begin{align*}
		\text{D}_{\theta_G}\rl_i(\bar{\theta}) &=  \mat{z_i(\bar{\theta})^\T\rQ(\bar{\theta}) & 0 & 0 \\ 0 & \ddots & 0 \\ 0 & 0 & z_i(\bar{\theta})^\T\rQ(\bar{\theta})} \\ &= I_m \otimes  z_i(\bar{\theta})^\T\rQ(\bar{\theta}),
	\end{align*}
	where $z_i(\bar{\theta})\in\C^r$ denotes the $i$-th column in $Z(\bar{\theta})$. This is also the case for more complex derivatives that involve the differentiation of the eigenvalue problem in \eqref{eq:red_evp}.
\end{rem}
Here, we highlight some important advantages of the pole-residue framework compared to recently proposed methods that are formulated in the Lyapunov framework (see~\cite{Sato2018,Jiang2019}), in particular for pH-DAEs. These methods require the solution of large-scale Lyapunov equations for the evaluation of $\losst$ and its gradient. Currently no structure-preserving Lyapunov-based methods exist for pH-DAEs; see \cite{SchM22_ppt} for new Lyapunov-based formulations of pH-DAEs. If the strictly proper part of the transfer function can be easily decoupled from the constant polynomial part, then for pH-DAEs with index one, the existing methods for pH-ODEs may be applied to this part.
However, if the splitting into the strictly proper and polynomial part has first to be computed via a factorization method, then the sparsity patterns of the original pH-DAE may be lost
which complicates the repetitive solution of Lyapunov equations for these systems in the large-scale setting. We highlight that the pole-residue framework only requires evaluations of $\pHtf_{\rm sp}$. Since we have that
\begin{equation}
	\pHtf_{\rm sp}(s) = \pHtf(s) - (S_0-N_0),
\end{equation}
for all $s\in\C$, we may work directly with the sparse matrices of the original pH-DAE and do not require the solution of large-scale Lyapunov equations.
\begin{algorithm}[tbh]
	\LinesNumbered
	\SetAlgoLined
	\DontPrintSemicolon
	\SetKwInOut{Input}{Input}\SetKwInOut{Output}{Output}
	\Input{FOM transfer function $\pHtf\in\mathcal{RH}_{\infty}^{\enu \times \enu}$, reduced order $r\in\mathbb{N}$.
	}
	\Output{Reduced pH-ODE of order $r$}
	Compute $S_0,\,N_0$ as in \eqref{eq:S0}--\eqref{eq:N0}.\;
	Initialize $\theta_0$ s.t. $\rS(\theta_0)=S_0$, $\rN(\theta_0)=N_0$.\;
	Solve \vspace{-12pt}\begin{align*}
	    \theta_{\rm fin} &= \argmin\limits_{\substack{\theta\in \R^{n_\theta}}} \losst(\theta; \pHtf) \\
	    &\text{s.t. } \rS(\theta)=S_0,\, \rN(\theta)=N_0.
	\end{align*}\vspace{-18pt} \; 
	Construct the ROM with $\theta_{\rm fin}$ as in Lemma \ref{lem:PHParam}.\;
	\caption{PROPT-$\mathcal{H}_2$}
	\label{alg:H2OPT}
\end{algorithm}

Since the $\cH_2$ optimization problem is non-convex, the choice of the initial parameter vector $\theta_0$ will generally impact the fidelity of the final ROM obtained by Algorithm~\ref{alg:H2OPT}. Simple initialization strategies are, for instance, choosing $\theta_0$ randomly or using IRKA-PH (see~\cite{Moser2020, Sato2018}), which generally converges very quickly. Here, we propose another approach that may use \emph{unstructured} ROMs for initialization which is based on the following parameterization.
\begin{lem}\label{lem:pHInit}
	 Let $(\widetilde{A},\widetilde{B},\widetilde{C},\widetilde{D})$ be a ROM of state-space dimension $r$ such that $\widetilde{D} = S_0-N_0$ and such that $\widetilde{A}$ has all its eigenvalues in the open left half of the complex plane. Let  $\theta_G \in\R^{r\cdot m}$ and $\theta_K \in\R^{r\cdot p}$ be two parameter vectors and define the matrix-valued functions
	 \begin{align*}
	    \iG(\theta_G) &:= \vtf_m(\theta_G), \\
	 	\rK(\theta_K) &:= \vtf_q(\theta_K).
	 \end{align*}
 	 Let $\iQ(\theta_K)>0$ solve the Lyapunov equation
 	 \begin{equation}
 	 	\widetilde{A}^\T\iQ(\theta_K)+\iQ(\theta_K)\widetilde{A}+\rK(\theta_K)\rK(\theta_K)^\T = 0,
 	 \end{equation}
  	and define 
  	\begin{align*}
  		\iJ(\theta_K)&=\frac{1}{2}\left(\widetilde{A} \iQ(\theta_K)^{-1}-\iQ(\theta_K)^{-1}\widetilde{A}^\T\right), \\
  		\iR(\theta_K)&=-\frac{1}{2}\left(\widetilde{A} \iQ(\theta_K)^{-1}+\iQ(\theta_K)^{-1}\widetilde{A}^\T\right).
  	\end{align*}
  	Then the parametric system
  	\begin{align}
  		\label{eq:pHParamInit}
  		\pHsysr(\theta_G,\theta_K):
  		\begin{cases}
  			\!\begin{aligned}
  				\dot x_r(t)= &\big(\iJ(\theta_K)-\iR(\theta_K))\iQ(\theta_K\big)\rx(t) \\ &\quad \quad + \iG(\theta_G) u(t),
  			\end{aligned}\\
  			\!\begin{aligned}
  		 	\ry(t)= &\iG(\theta_G)^\T \iQ(\theta_K) \rx(t)\\ &\quad \quad+ (S_0-N_0) u(t)
  			\end{aligned}\\
  		\end{cases}
  	\end{align}
  	is a pH-ODE system with
  	\begin{equation*}
  		\big(\iJ(\theta_K)-\iR(\theta_K)\big)\iQ(\theta_K) = \widetilde{A}.
  	\end{equation*}  	 
\end{lem}
Let $\widetilde{H}$ denote the transfer function of the (possibly unstructured) ROM $\big(\widetilde{A},\widetilde{B},\widetilde{C},\widetilde{D}\big)$ with ${\widetilde{\pHtf}(s)=\sum\limits_{i=1}^{r} \frac{\widetilde{\rl}_i\widetilde{\rr}_i^\T}{s-\widetilde{\reig}_i}+S_0-N_0}$ and ${\widetilde{\rl}_i,\,\widetilde{\rr}_i\in \C^m}$. 

Based on the parameterization in Lemma \ref{lem:pHInit}, we can then compute an initial pH model by minimizing the weighted sum of squared errors between the residuals in the Frobenius norm,  i.\,e.,  
\begin{equation*}
	\losst_0(\theta_G,\theta_K) := \sum_{i=1}^{r} \frac{1}{|\widetilde{\reig}_i|} \left\Vert \widetilde{\rl}_i\widetilde{\rr}_i^\T - \rl_i(\theta_G,\theta_K)\rr_i(\theta_G,\theta_K)^\T\right\Vert_{\rm F}^2,
\end{equation*}
where
\begin{align*}
	\rl_i(\theta_G,\theta_K) &= \iG(\theta_G)^\T\iQ(\theta_K) \widetilde{Z} e_i, \\
	\rr_i(\theta_G,\theta_K) &= \iG(\theta_G)^\T \widetilde{Z}^{-\T}e_i,
\end{align*}
for $i=1,\,\ldots,\,r$ and $\widetilde{Z}$ is, 
again under a diagonalizability assumption, obtained from the spectral decomposition
 \begin{equation*} 
	\widetilde{A} \widetilde{Z} = \widetilde{Z}\widetilde{\Lambda},
\end{equation*}
with $\widetilde{\Lambda} = \text{diag}(\widetilde{\reig}_1,\,...\,,\widetilde{\reig}_r)$.

Note that the computation of the gradient of $\losst_0$ is very simple, since it does not involve a differentiation of the eigenvalues or eigenvectors. While the partial gradients of $\rl_i(\cdot)$ and $\rr_i(\cdot)$ with respect to $\theta_G$ are straightforward, the partial gradients of $\iQ(\cdot)$ with respect to the $l$-th entry in $\theta_K$ is the solution of the (reduced-order) Lyapunov equation
\begin{multline*}
	\widetilde{A}^\T \frac{\partial\iQ(\theta_K)}{\partial \theta_{K,l}} + \frac{\partial\iQ(\theta_K)}{\partial \theta_{K,l}}\widetilde{A} \\ + \vtf_q(e_l)\rK(\theta_K)^\T+\rK(\theta_K)\vtf_q(e_l)^\T= 0,
\end{multline*}
where $e_l$ denotes the $l$-th standard basis vector of $\R^{r\cdot p}$. As the number of optimization parameters is reduced to $r(p+m)$, this initialization generally converges very quickly. In combination with Algorithm~\ref{alg:H2OPT}, this enables a \emph{two-step} approach with a more restrictive (yet simpler) pre-optimization of only the residuals and a subsequent (more complex) optimization of all system matrices.
\begin{rem}
	Note that the sample-based SOBMOR method can be tuned to compute a ROM with small $\htwo$ error as well. Instead of using $\loss$ in conjunction with the bisection method outlined in Algorithm~\ref{alg:bisection}, the integral in~\eqref{eq:htwonorm} can be approximated by means of an adaptive quadrature rule (see \cite[Algorithm~1]{Gonnet2010} for a template method). In this way, it is possible to compute the $\htwo$ error and its gradient with respect to the free ROM parameters in terms of the error transfer function at specific sample points. This makes it possible to use the same optimization techniques as in~\cite{Schwerdtner2020}. In particular, \cite[Theorem~3.1]{Schwerdtner2020} for the gradient computation can be reused. 
	
	We denote this method by SOBMOR-$\htwo$. Further details regarding the implementation of the adaptive integration are provided in the Appendix.
\end{rem}

2022-06-06T02:15:37.000
2206
2206.01608
en
https://arxiv.org/abs/2206.01608
We develop optimization-based structure-preserving model order reduction (MOR) methods for port-Hamiltonian (pH) descriptor systems of differentiation index one. Descriptor systems in pH form permit energy-based modeling and intuitive coupling of physical systems across different physical domains, scales, and accuracies. This makes pH models well-suited building-blocks for component-wise modeling of large system networks. In this context, it is often necessary to preserve the pH structure during MOR. We discuss current projection-based and structure-preserving MOR algorithms for pH systems and present a new optimization-based framework for that task. The benefits of our method include a simplified treatment of algebraic constraints and often a higher accuracy of the resulting reduced-order model, which is demonstrated by several numerical examples.
Optimization and Control (math.OC); Systems and Control (eess.SY); Dynamical Systems (math.DS)
Structure-Preserving Model Order Reduction for Index One Port-Hamiltonian Descriptor Systems
Qwen/Qwen-72B
1. YES
2. YES