Differences between Octave and Matlab

From Octave
(Redirected from Compatibility)
Jump to navigation Jump to search

People often ask

I wrote some code for Matlab, and I want to get it running under Octave. Is there anything I should watch out for?

or alternatively

I wrote some code in Octave, and want to share it with Matlab users. Is there anything I should watch out for?

which is not quite the same thing. There are still a number of differences between Octave and Matlab, however in general differences between the two are considered as bugs. Octave might consider that the bug is in Matlab and do nothing about it, but generally functionality is almost identical. If you find an important functional difference between Octave behavior and Matlab, then you should send a description of this difference (with code illustrating the difference, if possible) to http://bugs.octave.org.

Furthermore, Octave adds a few syntactical extensions to Matlab that might cause some issues when exchanging files between Matlab and Octave users.

As both Octave and Matlab are under constant development, the information in this section is subject to change.

You should also look at the pages http://packages.octave.org and http://octave.sourceforge.net/docs.html that have a function reference that is up to date. You can use this function reference to see the number of octave functions that are available and their Matlab compatibility.

Graphical Curve Fitting Tool[edit]

Currently, Octave lacks a graphical curve fitting tool such as Matlab's curvefit tool. [1]

You can use polyfit or the packages signal or optim to achieve these task but these are not graphical tools.

Nested Functions[edit]

Octave has limited support for nested functions since version 3.8.0. That is

function y = foo (x)
  y = bar(x)
  function y = bar (x)
    y = ...;
  end
end

is equivalent to

function y = foo (x)
  y = bar(x)
endfunction

function y = bar (x)
  y = ...;
endfunction

The main difference with Matlab is a matter of scope. While nested functions have access to the parent function's scope in Matlab, no such thing is available in Octave, due to how Octave essentially “un-nests” nested functions.

The authors of Octave consider the nested function scoping rules of Matlab to be more problems than they are worth as they introduce difficult to find bugs as inadvertently modifying a variable in a nested function that is also used in the parent is particularly easy for those not attentive to detail.

Differences in core syntax[edit]

There are a few core Matlab syntaxes that are not accepted by Octave, these being

  • Some limitations on the use of function handles. The major difference is related to nested function scoping rules (as above) and their use with function handles.
  • Some limitations of variable argument lists on the LHS of an expression, though the most common types are accepted.
  • Matlab classdef object oriented programming is only partially supported, see classdef for details.

Differences in core functions[edit]

A large number of the Matlab core functions (i.e. those that are in the core and not a toolbox) are implemented, and certainly all of the commonly used ones. There are a few functions that aren't implemented, usually to do with specific missing Octave functionality (GUI, DLL, Java, ActiveX, DDE, web, and serial functions). Some of the core functions have limitations that aren't in the Matlab version. For example the sprandn function can not force a particular condition number for the matrix like Matlab can. Another example is that testing and the runtests function work differently in Matlab and Octave.

Just-In-Time compiler[edit]

Matlab includes a "Just-In-Time" compiler. This compiler allows the acceleration of for-loops in Matlab to almost native performance with certain restrictions. The JIT must know the return type of all functions called in the loops and so you can't include user functions in the loop of JIT optimized loops. Octave has a not fully functional JIT compiler. For this reason you must vectorize your code as much as possible. The MathWorks themselves have a good document discussing vectorization at http://www.mathworks.com/support/tech-notes/1100/1109.html.

Compiler[edit]

On a related point, there is no Octave compiler, and so you can't convert your Octave code into a binary for additional speed or distribution.

Graphic handles[edit]

The support for graphics handles is converging towards full compatibility. If you notice any incompatibilities, please report a bug.

GUI functions[edit]

The support for Matlab compatible GUI functions was added in Octave version 3.6.0 and is converging towards full compatibility. If you notice any incompatibilities, please report a bug.

Simulink[edit]

Octave itself includes no Simulink support. Typically the Simulink models lag research and are less flexible, so shouldn't really be used in a research environment. However, some Matlab users that try to use Octave complain about this lack.

MEX-Files[edit]

Octave includes an API to the Matlab MEX interface. However, as MEX is an API to the internals of Matlab and the internals of Octave differ from Matlab, there is necessarily a manipulation of the data to convert from a MEX interface to the Octave equivalent. This is notable for all complex matrices, where Matlab stores complex arrays as real and imaginary parts, whereas Octave respects the C99/C++ standards of co-locating the real/imag parts in memory. Also due to the way Matlab allows access to the arrays passed through a pointer, the MEX interface might require copies of arrays (even non complex ones).

Block comments[edit]

Block comments denoted by #{ and #} markers (or %{ and %}) are supported by Octave with some limitations. The major limitation is that block comments are not supported within [] or {}.

Mat-File format[edit]

There are some differences in the mat v5 file format accepted by Octave. Matlab recently introduced the "-V7.3" save option which is an HDF5 format which is particularly useful for 64-bit platforms where the standard Matlab format can not correctly save variables. Octave accepts HDF5 files, but is not yet compatible with the "-v7.3" versions produced by Matlab.

Although Octave can load inline function handles saved by Matlab, it can not yet save them.

Finally, some multi-byte Unicode characters aren't yet treated in mat-files.

Profiler[edit]

Thanks to Daniel Kraft's 2011 Google Summer of Code project, Octave has a profiler since version 3.6.0.

Toolboxes[edit]

Octave is a community project and so the toolboxes that exist are donated by those interested in them through Octave Forge. These might be lacking in certain functionality relative to the Matlab toolboxes, and might not exactly duplicate the Matlab functionality or interface.

Short-circuit & and | operators[edit]

The & and | operators in Matlab short-circuit when included in a condition (e.g. an if or while statement) and not otherwise. In Octave only the && and || short circuit. Note that this means that

if (a | b)
  ...
end

and

t = a | b;
if (t)
  ...
end

have different semantics in Matlab. This is really a Matlab bug, but there is too much code out there that relies on this behavior to change it. Prefer the && and || operators in if statements if possible.

Note that the difference with Matlab is also significant when either argument is a function with side effects or if the first argument is a scalar and the second argument is an empty matrix. For example, note the difference between

t = 1 | [];          ## results in [], so...
if (t) 1, end        ## in if ([]), this is false.

and

if (1 | []) 1, end   ## short circuits so condition is true.

In the latter case, Octave displays since version 4.0.0 a warning:

 warning: Matlab-style short-circuit operation performed for operator |

Another case that is documented in the Matlab manuals is that

t = [1, 1] | [1, 2, 3];          ## error
if ([1, 1] | [1, 2, 3]) 1, end   ## OK

Also Matlab requires the operands of && and || to be scalar values but Octave does not (it just applies the rule that for an operand to be considered true, every element of the object must be nonzero or logically true).

Finally, note the inconsistence of thinking of the condition of an if statement as being equivalent to all(X(:)) when X is a matrix. This is true for all cases EXCEPT empty matrices:

if ([0, 1]) == if (all ([0, 1]))   ==>  i.e., condition is false.
if ([1, 1]) == if (all ([1, 1]))   ==>  i.e., condition is true.

However,

if ([])

is not the same as

if (all ([]))

because, despite the name, the all is really returning true if none of the elements of the matrix are zero, and since there are no elements, well, none of them are zero. This is an example of vacuous truth. But, somewhere along the line, someone decided that if ([]) should be false. The Mathworks probably thought it just looks wrong to have [] be true in this context even if you can use logical gymnastics to convince yourself that "all" the elements of an empty matrix are nonzero. Octave however duplicates this behavior for if statements containing empty matrices.

Solvers for singular, under- and over-determined matrices[edit]

Matlab's solvers as used by the operators mldivide \ and mrdivide /, use a different approach than Octave's in the case of singular, under-, or over-determined matrices. In the case of a singular matrix, Matlab returns the result given by the LU decomposition, even though the underlying solver has flagged the result as erroneous. Octave has made the choice of falling back to a minimum norm solution of matrices that have been flagged as singular which arguably is a better result for these cases.

In the case of under- or over-determined matrices, Octave continues to use a minimum norm solution, whereas Matlab uses an approach that is equivalent to

function x = mldivide (A, b)
  m = rows (A);
  [Q, R, E] = qr (A);
  x = [A \ b, E(:, 1:m) * (R(:, 1:m) \ (Q' * b))]
end

While this approach is certainly faster and uses less memory than Octave's minimum norm approach, this approach seems to be inferior in other ways.

A numerical question arises: how big can the null space component become, relative to the minimum-norm solution? Can it be nicely bounded, or can it be arbitrarily big? Consider this example:

m = 10;
n = 10000;
A = ones (m, n) + 1e-6 * randn (m, n);
b = ones (m, 1) + 1e-6 * randn (m, 1);
norm (A \ b)

while Octave's minimum-norm values are about 3e-2, Matlab's results are 50-times larger. For another issue, try this code:

m = 5;
n = 100;
j = floor (m * rand (1, n)) + 1;
b = ones (m, 1);
A = zeros (m, n);
A(sub2ind(size(A),j,1:n)) = 1;
x = A \ b;
[~,p] = sort (rand (1, n));
y = A(:,p) \ b;
norm (x(p) - y)

It shows that unlike in Octave, mldivide in Matlab is not invariant with respect to column permutations. If there are multiple columns of the same norm, permuting columns of the matrix gets you different result than permuting the solution vector. This will surprise many users.

Since the mldivide \ and mrdivide / operators are often part of a more complex expression, where there is no room to react to warnings or flags, it should prefer intelligence (robustness) to speed, and so the Octave developers are firmly of the opinion that Octave's approach for singular, under- and over-determined matrices is a better choice than Matlab's.

Octave extensions[edit]

The extensions in Octave over MATLAB syntax are very useful, but might cause issues when sharing with Matlab users. A list of the major extensions that should be avoided to be compatible with Matlab are:

Comments in Octave can be marked with #. This allows POSIX systems to have the first line as #! octave -q and mark the script itself executable. MATLAB doesn't have this feature due to the absence of comments starting with #".

Code blocks like if, for, while, etc can be terminated with block specific terminations like endif. MATLAB doesn't have this and all blocks must be terminated with end.

Octave has a lisp-like unwind_protect block that allows blocks of code that terminate in an error to ensure that the variables that are touched are restored. You can do something similar with try/catch combined with rethrow (lasterror ()) in Matlab, however rethrow and lasterror are only available in Octave 2.9.10 and later. MATLAB 2008a also introduced OnCleanUp that is similar to unwind_protect, except that the object created by this function has to be explicitly cleared in order for the cleanup code to run.

Note that using try/catch combined with rethrow (lasterror ()) cannot guarantee that global variables will be correctly reset, as it won't catch user interrupts with Ctrl-C. For example

global a
a = 1;
try
  _a = a;
  a = 2
  while true
  end
catch
  fprintf ('caught interrupt\n');
  a = _a;
  rethrow (lasterror());
end

compared to

global a
a = 1;
unwind_protect
  _a = a;
  a = 2
  while true
  end
unwind_protect_cleanup
  fprintf ('caught interrupt\n');
  a = _a;
end

Typing Ctrl-C in the first case returns the user directly to the prompt, and the variable a is not reset to the saved value. In the second case the variable a is reset correctly. Therefore Matlab gives no safe way of temporarily changing global variables.

Indexing can be applied to all objects in Octave and not just variables. Therefore sin(x)(1:10) for example is perfectly valid in Octave but not Matlab. To do the same in Matlab you must do y = sin(x); y = y([1:10]);

Octave has the operators ++, –-, -=, +=, *=, etc. As MATLAB doesn't, if you are sharing code these should be avoided.

Character strings in Octave can be denoted with double or single quotes. There is a subtle difference between the two in that escaped characters like \n (newline), \t (tab), etc are interpreted in double quoted strings but not single quoted strings. This difference is important on Windows platforms where the \ character is used in path names, and so single quoted strings should be used in paths. MATLAB doesn't have double quoted strings and so they should be avoided if the code will be transferred to a MATLAB user.

What features are unique to Octave?[edit]

Although most of the Octave language will be familiar to Matlab users, it has some unique features of its own.

Functions defined on the command-line[edit]

Functions can be defined by entering code on the command line, a feature not supported by Matlab. For example, you may type:

>> function s = hello_string (to_who)
> ## Say hello
> if nargin<1, to_who = "World"; end
> s = ["Hello ",\
>      to_who];
> endfunction
>> hello_string ("Moon")
ans = Hello Moon

As a natural extension of this, functions can also be defined in script files (m-files whose first non-comment line isn't function out = foo (...))

Note: MATLAB R2016b added the ability to define functions in script files.

Comments with #[edit]

The pound character, #, may be used to start comments, in addition to %. See the previous example. The major advantage of this is that as # is also a comment character for unix script files, any file that starts with a string like #! /usr/bin/octave -q will be treated as an octave script and be executed by octave.

Strings delimited by double quotes "[edit]

In 2016, Matlab introduced String Arrays, that are initialized by using double quoted strings, and are not implemented in Octave yet. In Octave double-quoted strings include backslash interpretation (like C++, C, and Perl) while single quoted are uninterpreted (like Matlab and Perl).

Line continuation by backslash[edit]

Lines can be continued with a backslash, \, in addition to three points ... as in Matlab.

Informative block closing[edit]

You may close function, for, while, if, ... blocks with endfunction, endfor, endwhile, ... keywords in addition to using end. As with Matlab, the end (or endfunction) keyword that marks the end of a function defined in a .m file is optional.

Coherent syntax[edit]

Indexing other things than variables is possible, as in:

>> [3 1 4 1 5 9](3)
ans = 4
>> cos([0 pi pi/4 7])(3)
ans = 0.70711

In Matlab, it is for example necessary to assign the intermediate result cos([0 pi pi/4 7]) to a variable before it can be indexed again.

Exclamation mark as not operator[edit]

The exclamation mark ! (aka “Bang!”) is a negation operator, just like the tilde ~:

>> if ! strcmp (program_name, "octave"),
>   "It's an error"
> else
>   "It works!"
> end
ans = It works!

Note however that Matlab uses the ! operator for shell escapes, for which Octave requires using the system command.

Increment and decrement operators[edit]

If you like the ++, += etc operators, rejoice! Octave includes the C-like increment and decrement operators ++ and -- in both their prefix and postfix forms, in addition to +=, -=, *=, /=, ^=,.+=,.-=,.*=, ./= and .^=.

For example, to pre-increment the variable x, you would write ++x. This would add one to x and then return the new value of x as the result of the expression. It is exactly the same as the expression x = x + 1.

To post-increment a variable x, you would write x++. This adds one to the variable x, but returns the value that x had prior to incrementing it. For example, if x is equal to 2, the result of the expression x++ is 2, and the new value of x is 3.

For matrix and vector arguments, the increment and decrement operators work on each element of the operand.

Unwind-protect[edit]

In addition to try-catch blocks, Octave supports an alternative form of exception handling modeled after the unwind-protect form of Lisp. The general form of an unwind_protect block looks like this:

unwind_protect
  body
unwind_protect_cleanup
  cleanup
end_unwind_protect

Where body and cleanup are both optional and may contain any Octave expressions or commands. The statements in cleanup are guaranteed to be executed regardless of how control exits body.

The unwind_protect statement is often used to reliably restore the values of global variables that need to be temporarily changed.

Matlab can be made to do something similar with their onCleanup function that was introduced in 2008a. Octave also has onCleanup since version 3.4.0.

Built-in ODE and DAE solvers[edit]

Octave includes LSODE, DASSL and DASPK for solving systems of stiff ordinary differential and differential-algebraic equations. These functions are built in to the interpreter.

Do-Until loop structure[edit]

Similar to the do-while loop in C and C++, Octave allows a do-until loop which does not exist in Matlab:

x = 0
do
  x += 1;
until (x == 10)

Broadcasting[edit]

Borrowed from other languages, octave broadcasting allows easy and readable vectorization.

f = (1:0.1:2);
# put angular frequencies on the first dimension to prepare broadcasting
omega = 2 * pi * f(:);
# time is already on the second dimension (row vector)
t = 0:0.02:2;
# the resulting s will be a 2-dimensional array
s = sin(omega .* t);
# which can be displayed as
pcolor(t, f, s)
xlabel("t (s)")
ylabel("f (Hz)")

Note: Automatic expansion of dimensions was added to MATLAB R2016b.

Documentation strings[edit]

Octave allows extensive formatting of the help string of functions using Texinfo. The effect on the online documentation is relatively small, but makes the help string of functions conform to the help of Octave’s own functions. However, the effect on the appearance of printed or online documentation will be greatly improved.

Test functions[edit]

Octave allows to add self-tests to user defined functions. Tests are put after function definition in specially commented block.

function mult = a(val)
  mult = val.*2;
endfunction
%!test
%! assert (a(3), 6);

Such a function can be tested for valid outputs by following code:

 >> test a
 PASSES 1 out of 1 test

Demonstration Functions[edit]

Example code block can be part of function file in a similar manner as test functions. For example to run demo for function multinom of package specfun, use:

  demo multinom

Powerful assert[edit]

Function assert have extended input possibilities.

How to Programmatically Identify If You're Running Octave[edit]

Because you may want to write cross-platform compatible code but make use of Octave-specific features when they're available, it may be necessary to programmatically detect the program environment. The following function can be used to detect whether it is being run in Octave. Due to the persistent variable it can be called repeatedly without a heavy performance hit.

   %% Return: true if the environment is Octave.
   %%
   function retval = isOctave
     persistent cacheval;  % speeds up repeated calls
   
     if isempty (cacheval)
       cacheval = (exist ("OCTAVE_VERSION", "builtin") > 0);
     end
   
     retval = cacheval;
   end

External links[edit]