missing functions

[asobouti]

Hi Elmar,
I was going to use your marten_covariance function but I caught by two missed function check_boolean and scaled_distance. Can I find them here?

[litvinen]

Hi Alireza, I see these functions, they are there. Probably matlab just doesn't see them.

[litvinen]

function dist=scaled_distance(x1, x2, l, smooth, sqr)
% SCALED_DISTANCE Computes the scaled distance between points.
% DIST=SCALED_DISTANCE(X1, X2) computes the Euclidean distance between
% points given in X1 and X2.
%
% DIST=SCALED_DISTANCE(X1, X2, L) computes the Euclidean distance between
% X1 and X2 scaled by L where L can be scalar or a vector of the same
% dimension as the points in X1 and X2.
%
% DIST=SCALED_DISTANCE(X1, X2, L, SMOOTH) computes the smoothed
% Euclidean distance where the smoothing function is given by
% D_SMOOTH=SQRT(D+SMOOTH^2)-SMOOTH.
%
% DIST=SCALED_DISTANCE(X1, X2, L, SMOOTH, SQR) returns the square of the
% Eudclidean smoothed, scaled distance. This saves computing the square
% root if only the square is needed by the calling code anyway.
%
% Note: if an empty array is specified for a parameter the default value
% is used. DIST=SCALED_DISTANCE(X1, X2, [], [], TRUE) computes the square
% of the distance with parameters L and SMOOTH set to their defaults 1
% and 0.
%
% Example (run)
% d=linspace(-1,1);
% plot(d, scaled_distance(d, [], [], 0)); hold all;
% plot(d, scaled_distance(d, [], [], 1));
% plot(d, scaled_distance(d, [], [], 2));
% plot(d, scaled_distance(d, [], 0.3, 0));
% plot(d, scaled_distance(d, [], 0.3, 1));
% plot(d, scaled_distance(d, [], 0.3, 2)); hold off;
% axis equal;
% See also

% Elmar Zander
% Copyright 2010, Inst. of Scientific Computing, TU Braunschweig
%
% This program is free software: you can redistribute it and/or modify it
% under the terms of the GNU General Public License as published by the
% Free Software Foundation, either version 3 of the License, or (at your
% option) any later version.
% See the GNU General Public License for more details. You should have
% received a copy of the GNU General Public License along with this
% program. If not, see http://www.gnu.org/licenses/.

if nargin<5 || isempty(sqr)
sqr=false;
end
if nargin<4 || isempty(smooth)
smooth=0;
end
if nargin<3 || isempty(l)
l=1;
end

if isempty(x2)
dx=x1;
else
dx=x1-x2;
end

d=size(dx,1);
if isscalar(l)
dx=dx/l;
elseif length(l)==d
invl=diag(1./l);
dx=invl*dx;
else
error( 'statistics:scaled_distance', 'Size of cov_length vector doesn''t match dimension. Transposed?' );
end
dist2=sum( dx.*dx, 1);

if sqr
if smooth>0
dist=sqrt(dist2+smooth^2)-smooth;
dist=dist.*dist;
else
dist=dist2;
end
else
if smooth>0
dist=sqrt(dist2+smooth^2)-smooth;
else
dist=sqrt(dist2);
end
end

[litvinen]

function ok=check_boolean( ok, message, filename, varargin )
% CHECK_BOOLEAN Check whether condition is true on input.
% OK=CHECK_BOOLEAN( COND, MESSAGE, MFILENAME ) checks whether the given
% condition is true. If not an error message is printed and the program
% is aborted.
%
% Note: pass mfilename literally for the last argument (i.e. pass the
% return value of the buildin function 'mfilename' which tells you the
% name of the current script, and is thus exactly what you want.)
%
% Example (run)
% function my_function( str )
%
% check_boolean( strcmp(str,str(end:-1:1)), 'str must be a palindrome', mfilename );
%
% See also CHECK_RANGE, CHECK_CONDITION, CHECK_UNSUPPORTED_OPTIONS

% Elmar Zander
% Copyright 2007, Institute of Scientific Computing, TU Braunschweig.
%
% This program is free software: you can redistribute it and/or modify it
% under the terms of the GNU General Public License as published by the
% Free Software Foundation, either version 3 of the License, or (at your
% option) any later version.
% See the GNU General Public License for more details. You should have
% received a copy of the GNU General Public License along with this
% program. If not, see http://www.gnu.org/licenses/.

% save some time by returning immediately if ok is true
if ok
return
end

options=varargin2options( varargin );
[mode,options]=get_option( options, 'mode', 'debug' );
[depth,options]=get_option( options, 'depth', 1 );
check_unsupported_options( options, mfilename );

if isempty(filename)
filename='util:check';
end

switch mode
case 'error'
stack=dbstack('-completenames');
err_struct.message=sprintf( '%s: %s', filename, message );
err_struct.identifier=[filename ':check_failed'];
err_struct.stack=stack((1+depth):end);
% using 'rethrow' instead of 'error' suppresses incorrect
% reports of error position in 'check_condition' (the 'error'
% function disregards the stack for display of the error
% position)
rethrow( err_struct );
case 'warning'
if exist('dbstack') %#ok
fprintf('Warning: %s: %s\n', filename, message );
stack=dbstack('-completenames');
for i=(1+depth):length(dbstack)
%fprintf( ' In %s at line %d\n', stack(i).file, stack(i).name, stack(i).line );
fprintf( ' In %s at line %d\n', stack(i).file, stack(i).line, stack(i).name, stack(i).line );
end
else
warning([filename ':check_failed'], '%s: %s', filename, message );
end
case 'print'
if exist('dbstack') %#ok
stack=dbstack('-completenames');
i=(1+depth);
fprintf( 'Check failed in %s at line %d: %s\n', stack(i).file, stack(i).line, stack(i).name, stack(i).line, message );
else
fprintf( 'Check failed in %s: %s\n', filename, message );
end
case 'debug'
fprintf( '\nCheck failed in: %s\n', filename);
fprintf( 2, 'Reason: %s\n', message );
cmd=repmat( 'dbup;', 1, depth );
fprintf( 'Use the stack to get to the place the assertion failed to \n', cmd );
fprintf( 'investigate the error. Then press F5 to continue or stop debugging.\n' )
keyboard;
end

if nargout==0
clear ok;
end

[ezander]

Hi Alireza!
Those functions are in statistics/private and util respectively. If you run startup or sglib_startup from the main directory the paths should be set correctly and the functions should then be found automatically.
Hope that helps. Cheers
Elmar

[ezander]

@litvinen Hi Alex! Thanks for helping out! 👍