horzcat_domains.m File Reference

Horizontally concatenate vectors resulting from a multi-domain ODE integration. Eliminate non-unique points by excluding points in between each data set (for steps 1 .. N-1, remove last point, so it's right continuous*) useful for interpolation. Since it's varargin, call via argument expansion:

horzcat_domains(data.steps.ts) horzcat(dataCell{:})

Example:

horzcat_domains([1, 2, 3; 10, 20, 30], [3, 4, 5, 6; 30, 40, 50 60])

1 2 3 4 5 6 10 20 30 40 50 60 For now, this will reject any empty matrices, since that doesn't make sense since an initial condition is always present. Can test resulting thing with (there should be no indices)

[row, col] = find(diff(results, 1, 2) == 0); [row, col]

More...

Functions

function X = horzcat_domains (varargin)
 Horizontally concatenate vectors resulting from a multi-domain ODE integration. Eliminate non-unique points by excluding points in between each data set (for steps 1 .. N-1, remove last point, so it's right continuous*) useful for interpolation. Since it's varargin, call via argument expansion:

horzcat_domains(data.steps.ts) horzcat(dataCell{:})

Example:

horzcat_domains([1, 2, 3; 10, 20, 30], [3, 4, 5, 6; 30, 40, 50 60])

1 2 3 4 5 6 10 20 30 40 50 60 For now, this will reject any empty matrices, since that doesn't make sense since an initial condition is always present. Can test resulting thing with (there should be no indices)

[row, col] = find(diff(results, 1, 2) == 0); [row, col]

More...
 

Detailed Description

Horizontally concatenate vectors resulting from a multi-domain ODE integration. Eliminate non-unique points by excluding points in between each data set (for steps 1 .. N-1, remove last point, so it's right continuous*) useful for interpolation. Since it's varargin, call via argument expansion:

horzcat_domains(data.steps.ts) horzcat(dataCell{:})

Example:

horzcat_domains([1, 2, 3; 10, 20, 30], [3, 4, 5, 6; 30, 40, 50 60])

1 2 3 4 5 6 10 20 30 40 50 60 For now, this will reject any empty matrices, since that doesn't make sense since an initial condition is always present. Can test resulting thing with (there should be no indices)

[row, col] = find(diff(results, 1, 2) == 0); [row, col]

  • Right continuity, see [Grizzle p85]

Function Documentation

◆ horzcat_domains()

function X = horzcat_domains (   varargin)

Horizontally concatenate vectors resulting from a multi-domain ODE integration. Eliminate non-unique points by excluding points in between each data set (for steps 1 .. N-1, remove last point, so it's right continuous*) useful for interpolation. Since it's varargin, call via argument expansion:

horzcat_domains(data.steps.ts) horzcat(dataCell{:})

Example:

horzcat_domains([1, 2, 3; 10, 20, 30], [3, 4, 5, 6; 30, 40, 50 60])

1 2 3 4 5 6 10 20 30 40 50 60 For now, this will reject any empty matrices, since that doesn't make sense since an initial condition is always present. Can test resulting thing with (there should be no indices)

[row, col] = find(diff(results, 1, 2) == 0); [row, col]

  • Right continuity, see [Grizzle p85]