Smooth Takagi decomp. of a complex symmetric matrix function

Numerical computation of a smooth Takagi decomposition of an n-by-n complex symmetric matrix valued function of one real parameter

This function numerically computes the Takagi decomposition of an n-by-n complex symmetric matrix-valued function FUN of one real parameter. The singular values are arranged in decreasing order and are assumed to be distinct and non-zero for all values of the parameter. It is possible to continue a subset of the singular values and corresponding singular vectors.

A typical call to takagiCont is:

[Tout,Uout,Sout,flag]=takagiCont(f,tspan,params)

See the script example_takagiCont.m for an example of how to use the function.

The command "help takagiCont" displays information and functionality of the software.

Please cite the references below if you use this software.

Authors: Alessandra Papini and Alessandro Pugliese