Difference between revisions of "Hess reorder.m"
(Sync syntax/arguments/outputs with current Spinach source) |
|||
| Line 1: | Line 1: | ||
{{DISPLAYTITLE:hess_reorder.m}} __NOTOC__ | {{DISPLAYTITLE:hess_reorder.m}} __NOTOC__ | ||
| + | |||
A utility to reorder a Hessian for use in [[grape.m|GRAPE]] optimal control problems. Swaps the Hessian ordering between control channel and time step variables being first or second. The Hessian elements correspond to the elements of a control waveform array ordered as: | A utility to reorder a Hessian for use in [[grape.m|GRAPE]] optimal control problems. Swaps the Hessian ordering between control channel and time step variables being first or second. The Hessian elements correspond to the elements of a control waveform array ordered as: | ||
<center> | <center> | ||
| Line 11: | Line 12: | ||
==Syntax== | ==Syntax== | ||
| − | + | hess=hess_reorder(hess,K,N) | |
| − | == | + | ==Outputs== |
| − | + | hess - reordered Hessian with N first then K. | |
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | + | david.goodwin@inano.au.dk | |
| + | ilya.kuprov@weizmann.ac.il | ||
| − | + | ==See also== | |
| − | |||
| − | |||
| − | |||
[[Optimal control module]] | [[Optimal control module]] | ||
''Version 1.9, authors: [[David Goodwin]]'' | ''Version 1.9, authors: [[David Goodwin]]'' | ||
| + | |||
| + | ==Parameters== | ||
| + | |||
| + | hess - Hessian matrix, | ||
| + | |||
| + | K - number of control channels. | ||
| + | |||
| + | N - number of time steps. | ||
Revision as of 15:03, 5 April 2026
A utility to reorder a Hessian for use in GRAPE optimal control problems. Swaps the Hessian ordering between control channel and time step variables being first or second. The Hessian elements correspond to the elements of a control waveform array ordered as:
\(\begin{bmatrix}X_1 & X_2 & \cdots & X_N\\ Y_1 & Y_2 & \cdots & Y_N\\ Z_1 & Z_2 & \cdots & Z_N\end{bmatrix}\)
where X,Y,Z are different control channels and the index enumerates the time discretization points. A re-ordered Hessian will correspond to a waveform ordered as:
\(\begin{bmatrix}X_1 & Y_1 & Z_1\\ X_2 & Y_2 & Z_2\\ \vdots & \vdots & \vdots\\ X_N & Y_N & Z_N\end{bmatrix}\)
Syntax
hess=hess_reorder(hess,K,N)
Outputs
hess - reordered Hessian with N first then K.
david.goodwin@inano.au.dk ilya.kuprov@weizmann.ac.il
See also
Version 1.9, authors: David Goodwin
Parameters
hess - Hessian matrix,
K - number of control channels.
N - number of time steps.
