Add proper copyright info and rewrite parts of code.

Author: genuinelucifer

inst/@sym/subsasgn.m

@@ -101,9 +101,9 @@
       end
 
     case '.'
-%      assert( isa(rhs, 'sym') || isa(rhs, 'symfun'))
-%      assert( ~isa(idx.subs, 'sym'))
-%      assert( ~isa(val, 'sym'))
+      assert( isa(rhs, 'sym'))
+      assert( ~isa(idx.subs, 'sym'))
+      assert( ~isa(val, 'sym'))
       val.(idx.subs) = rhs;
       out = val;
 

inst/@sym/sym.m

@@ -1,5 +1,6 @@
 %% Copyright (C) 2014-2016 Colin B. Macdonald
 %% Copyright (C) 2016 Lagu
+%% Copyright (C) 2017 Abhinav Tripathi
 %%
 %% This file is part of OctSymPy.
 %%

inst/@symfun/isequal.m

@@ -1,11 +1,56 @@
-function t = isequal(x, y)
+%% Copyright (C) 2017 Abhinav Tripathi
+%%
+%% This file is part of OctSymPy.
+%%
+%% OctSymPy 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.
+%%
+%% This software is distributed in the hope that it will be useful,
+%% but WITHOUT ANY WARRANTY; without even the implied warranty
+%% of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See
+%% the GNU General Public License for more details.
+%%
+%% You should have received a copy of the GNU General Public
+%% License along with this software; see the file COPYING.
+%% If not, see <http://www.gnu.org/licenses/>.
 
-  if (isa (x, 'symfun'))
-    x = x.sym;
-  end
-  if (isa (y, 'symfun'))
-    y = y.sym;
-  end
+%% -*- texinfo -*-
+%% @documentencoding UTF-8
+%% @defmethod  @@sym isequal (@var{f}, @var{g})
+%% @defmethodx @@sym isequal (@var{f}, @var{g}, @dots{})
+%% Test if contents of two or more symfuns are equal.
+%%
+%% Example:
+%% @example
+%% @group
+%% syms x
+%% f(x) = x;
+%% isequal (f(x), x)
+%%   @result{} 1
+%% @end group
+%% @end example
+%%
+%% @seealso{@@sym/isequal}
+%% @end defmethod
 
-  t = isequal(x, y);
-end
+function t = isequal (x, y, varargin)
+
+  t = isequal@sym (x, y, varargin{:});
+
+end
+
+
+%!test
+%! syms x;
+%! f(x) = x^2;
+%! g(x) = x^2;
+%! assert (isequal (f(x), g(x)))
+%! assert (isequal (f(x), x^2))
+
+%!test
+%! syms x y;
+%! f(x) = x^2;
+%! g(y) = y^2;
+%! assert (isequal (f(x), g(x)))

inst/@symfun/subsasgn.m

@@ -1,3 +1,32 @@
+%% Copyright (C) 2017 Abhinav Tripathi
+%%
+%% This file is part of OctSymPy.
+%%
+%% OctSymPy 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.
+%%
+%% This software is distributed in the hope that it will be useful,
+%% but WITHOUT ANY WARRANTY; without even the implied warranty
+%% of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See
+%% the GNU General Public License for more details.
+%%
+%% You should have received a copy of the GNU General Public
+%% License along with this software; see the file COPYING.
+%% If not, see <http://www.gnu.org/licenses/>.
+
+%% -*- texinfo -*-
+%% @documentencoding UTF-8
+%% @deftypeop  Method   @@sym {@var{f} =} subsasgn (@var{f}, @var{idx}, @var{rhs})
+%% @deftypeopx Operator @@sym {} {@var{f}(@var{i}) = @var{rhs}} {}
+%% @deftypeopx Operator @@sym {} {@var{f}(@var{i}, @var{j}) = @var{rhs}} {}
+%% @deftypeopx Operator @@sym {} {@var{f}(@var{i}:@var{j}) = @var{rhs}} {}
+%% @deftypeopx Operator @@sym {} {@var{f}(@var{x}) = @var{symexpr}} {}
+%% Assign value to a symbolic function [constructor].
+%%
+%% @end deftypeop
+
 function out = subsasgn (val, idx, rhs)
 
   switch idx.type
@@ -8,6 +37,7 @@
       %   x is idx.subs{1}
       % This also gets called for "syms f(x)"
 
+      disp('constructor called');
       all_syms = true;
       for i = 1:length(idx.subs)
         all_syms = all_syms && isa(idx.subs{i}, 'sym');
@@ -18,32 +48,18 @@
         all_Symbols = python_cmd (cmd, idx.subs);
       end
       if (all_syms && all_Symbols)
-	%% Make a symfun
+	      %% Make a symfun
         if (~isa(rhs, 'sym'))
           % rhs is, e.g., a double, then we call the constructor
           rhs = sym(rhs);
         end
         out = symfun(rhs, idx.subs);
-
-      else
-        %% Not symfun: e.g., f(double) = ..., f(sym(2)) = ...,
-        % convert any sym subs to double and do array assign
-        for i = 1:length(idx.subs)
-          if (isa(idx.subs{i}, 'sym'))
-            idx.subs{i} = double(idx.subs{i});
-          end
-        end
-	for i = 1:length(idx.subs)
-          if (~ sym.is_valid_index (idx.subs{i}))
-            error('OctSymPy:subsref:invalidIndices', ...
-                  'invalid indices: should be integers or boolean');
-          end
-	end
-        out = sym.mat_replace (val, idx.subs, sym(rhs));
       end
 
     case '.'
-      
+      assert (isa (rhs, 'symfun'))
+      assert (~ isa (idx.subs, 'symfun'))
+      assert (~ isa (val, 'symfun'))
       val.(idx.subs) = rhs;
       out = val;
 
@@ -54,4 +70,10 @@
       val
       error('broken');
   end
-end
+end
+
+
+%!test
+%! syms x;
+%! f(x) = x^2;
+%! assert (isa (f, 'symfun'))

inst/@symfun/symfun.m

@@ -1,4 +1,5 @@
 %% Copyright (C) 2014-2016 Colin B. Macdonald
+%% Copyright (C) 2017 Abhinav Tripathi
 %%
 %% This file is part of OctSymPy.
 %%