LYAPACK: lp_e(i1,i2,i3,i4,i5,i6,i7,i8,i9,i10,i11,i12,i13,i14,i15,i16,... - File Exchange

No BSD License

Highlights from
LYAPACK

... % Low rank Cholesky factor Newton method (explicit version LRCF-NM and
... % LRCF-ADI for solving the stable Lyapunov equation:
... % Using updated QR factorizations, this routine computes efficiently
as_l(tr,X) % Solves linear systems with the real, symmetric, negative definite matrix A,
as_l_d % Deletes the global data that has been generated by 'as_l_i'.
as_l_i
as_m(tr,X) % Evaluates matrix-matrix products with the real, symmetric matrix A,
as_m_d % Deletes the global data that has been generated by 'as_m_i'.
as_m_i(A)
as_pre(A,B,C) % Preprocessing of the system
as_pst(X,iprm) % Postprocessing. This routine is the counterpart of 'as_pre'. It
as_s(tr,X,i) % Solves shifted linear systems with the real, symmetric, negative definite
as_s_d(p) % Deletes the global data that has been generated by 'as_s_i'.
as_s_i(p)
au_l(tr,X) % Solves linear systems with the real matrix A or its transposed A':
au_l_d % Deletes the global data that has been generated by 'au_l_i'.
au_l_i
au_m(tr,X) % Evaluates matrix-matrix products with the real matrix A or its
au_m_d % Deletes the global data that has been generated by 'au_m_i'.
au_m_i(A)
au_pre(A,B,C) % Preprocessing of the system
au_pst(X,iprm) % Postprocessing. This routine is the counterpart of 'au_pre'. It
au_qmr_ilu_l(tr,X) % Solves linear systems with the real matrix A or its transpose A' using
au_qmr_ilu_l_d % Deletes the global data that has been generated by 'au_qmr_ilu_l_i'.
au_qmr_ilu_l_i % Generates the data used in 'au_qmr_ilu_l'. Data are stored in global
au_qmr_ilu_m(tr,X) % Evaluates matrix-matrix products with the real matrix A or its
au_qmr_ilu_m_d % Deletes the global data that has been generated by 'au_qmr_ilu_m_i'.
au_qmr_ilu_m_i(A,mc,max_it_qm... % Generates the data used in 'au_qmr_ilu_m'. Data are stored in global
au_qmr_ilu_s(tr,X,i) % Solves shifted linear systems with the matrix A or its transpose A' by
au_qmr_ilu_s_d(p) % Deletes the global data that has been generated by 'au_qmr_ilu_s_i'
au_qmr_ilu_s_i(p) % Generates the data used in 'au_qmr_ilu_s'. Data are stored in global
au_s(tr,X,i) % Solves shifted linear systems with the real matrix A or its
au_s_d(p) % Deletes the global data that has been generated by 'au_s_i'.
au_s_i(p)
fdm_2d_matrix(n0,fx_str,fy_st... % Generates the stiffness matrix A for the finite difference
fdm_2d_vector(n0,f_str)
lp_arn_m(name,Bf,Kf,k,r)
lp_arn_p(name,Bf,Kf,k,r)
lp_dspmr(name,B,C,ZB,ZC,max_o... % Model reduction method DSPMR (Dominant Subspace Projection Model
lp_e(i1,i2,i3,i4,i5,i6,i7,i8,... % Auxiliary routine for the execution of strings with MATLAB
lp_gnorm(Gs,m,q) % Given the transfer function sample matrix Gs, which contains
lp_lgfrq(min_freq,max_freq,no... % Delivers a vector of "frequency sampling points" which are
lp_lrsrm(name,B,C,ZB,ZC,max_o... % Model reduction method LRSRM (Low Rank Square Root Method)
lp_mnmx(rw,l0) % Suboptimal solution of the ADI minimax problem. The delivered parameter
lp_nrm( tp, name, Bf, Kf, G, ... % Computes efficiently either of the following norms provided that
lp_para(name,Bf,Kf,l0,kp,km,b... % Estimation of suboptimal ADI shift parameters for the matrix
lp_prm(A,E) % Reordering of SPARSE standard/generalized systems
lp_rcnrm( name, B, C0, R, Z ) % Computes efficiently the Riccati norm when the approximate solution
lp_s(p,set) % Computation of the maximal magnitude of the rational ADI function over
lp_trfia(freq,A,B,C,D,E) % Computes the transfer function for systems
msns_l(tr,X) % Solves linear systems with the symmetric, negative definite matrix A,
msns_l_d % Deletes the global data that has been generated by 'msns_l_i'.
msns_l_i
msns_m(tr,X) % Evaluates matrix-matrix products with the symmetric matrix A,
msns_m_d % Deletes the global data that has been generated by 'msns_m_i'.
msns_m_i(M,MU,N)
msns_pre(M,N,B,C) % Preprocessing of the system
msns_pst(Z,MU,iprm) % Postprocessing. This routine is the counterpart of 'msns_pre'. It
msns_s(tr,X,i) % Solves shifted linear systems with the symmetric, negative definite
msns_s_d(p) % Deletes the global data that has been generated by 'msns_s_i'.
msns_s_i(p)
munu_l(tr,X) % Solves linear systems with the real matrix A or its transposed A':
munu_l_d % Deletes the global data that has been generated by 'munu_l_i'.
munu_l_i
munu_m(tr,X) % Evaluates matrix-matrix products with the real matrix A or its
munu_m_d % Deletes the global data that has been generated by 'munu_m_i'.
munu_m_i(M,ML,MU,N)
munu_pre(M,N,B,C) % Preprocessing of the system
munu_pst(ZB,ZC,ML,MU,iprm) % Postprocessing. This routine is the counterpart of preprocessing
munu_s(tr,X,i) % Solves shifted linear systems with the real matrix A or its
munu_s_d(p) % Deletes the global data that has been generated by 'munu_s_i'.
munu_s_i(p)
demo_l1.m
demo_m1.m
demo_m2.m
demo_r1.m
demo_u1.m
demo_u2.m
demo_u3.m
lstartup.m
View all files

from LYAPACK by Volker Mehrmann
LYAPACK toolbox provides solutions for certain large scale problems related to Lyapunov equations.

lp_e(i1,i2,i3,i4,i5,i6,i7,i8,i9,i10,i11,i12,i13,i14,i15,i16,...

function s = lp_e(i1,i2,i3,i4,i5,i6,i7,i8,i9,i10,i11,i12,i13,i14,i15,i16,...
             i17,i18,i19,i20)
%
%  Auxiliary routine for the execution of strings with MATLAB 
%  expressions, which depend on numerical values.
%
%  Input:
%
%    i1,...,i20 strings or numerical values.    
%
%  Output:
%
%    s          string.
%
%
%  LYAPACK 1.0 (Thilo Penzl, November 1999)

ni = nargin;

if ni>=1, if ischar(i1), s = i1; else, s = num2str(i1); end, end 
if ni>=2, if ischar(i2), s = [s,i2]; else, s = [s,num2str(i2)]; end, end 
if ni>=3, if ischar(i3), s = [s,i3]; else, s = [s,num2str(i3)]; end, end 
if ni>=4, if ischar(i4), s = [s,i4]; else, s = [s,num2str(i4)]; end, end 
if ni>=5, if ischar(i5), s = [s,i5]; else, s = [s,num2str(i5)]; end, end 
if ni>=6, if ischar(i6), s = [s,i6]; else, s = [s,num2str(i6)]; end, end 
if ni>=7, if ischar(i7), s = [s,i7]; else, s = [s,num2str(i7)]; end, end 
if ni>=8, if ischar(i8), s = [s,i8]; else, s = [s,num2str(i8)]; end, end 
if ni>=9, if ischar(i9), s = [s,i9]; else, s = [s,num2str(i9)]; end, end 
if ni>=10, if ischar(i10), s = [s,i10]; else, s = [s,num2str(i10)]; end, end 
if ni>=11, if ischar(i11), s = [s,i11]; else, s = [s,num2str(i11)]; end, end 
if ni>=12, if ischar(i12), s = [s,i12]; else, s = [s,num2str(i12)]; end, end 
if ni>=13, if ischar(i13), s = [s,i13]; else, s = [s,num2str(i13)]; end, end 
if ni>=14, if ischar(i14), s = [s,i14]; else, s = [s,num2str(i14)]; end, end 
if ni>=15, if ischar(i15), s = [s,i15]; else, s = [s,num2str(i15)]; end, end 
if ni>=16, if ischar(i16), s = [s,i16]; else, s = [s,num2str(i16)]; end, end 
if ni>=17, if ischar(i17), s = [s,i17]; else, s = [s,num2str(i17)]; end, end 
if ni>=18, if ischar(i18), s = [s,i18]; else, s = [s,num2str(i18)]; end, end 
if ni>=19, if ischar(i19), s = [s,i19]; else, s = [s,num2str(i19)]; end, end 
if ni>=20, if ischar(i20), s = [s,i20]; else, s = [s,num2str(i20)]; end, end 
if ni>=21, error('Too many input arguments.'); end

Highlights from LYAPACK

Highlights from
LYAPACK