# HG changeset patch # User Michele Zaffalon # Date 1433911330 25200 # Node ID 507ccf8f10ed4500359fc48795bece5be535c8c4 # Parent 7fa170cc14fec67c0ef9d96eef2ae8faca77793c ellipke.m: Modify to accept row vectors (bug #45283). * ellipke.m: Don't reshape input m to a column vector. diff -r 7fa170cc14fe -r 507ccf8f10ed scripts/specfun/ellipke.m --- a/scripts/specfun/ellipke.m Tue Jun 09 16:25:45 2015 -0700 +++ b/scripts/specfun/ellipke.m Tue Jun 09 21:42:10 2015 -0700 @@ -94,7 +94,6 @@ endif sz = size (m); - m = m(:); if (! isreal (m)) error ("ellipke: M must be real"); elseif (any (m > 1)) @@ -127,8 +126,8 @@ mult_k = 1./sqrt (1 - m(idx_neg)); mult_e = sqrt (1 - m(idx_neg)); m(idx_neg) = -m(idx_neg) ./ (1 - m(idx_neg)); - a = ones (sum (idx), 1); b = sqrt (1 - m(idx)); + a = ones (size (b)); c = sqrt (m(idx)); f = 0.5; sum = f*c.^2;