mirror of
https://github.com/lcn2/calc.git
synced 2025-08-19 01:13:27 +03:00
a31078bbec
Some folks might think: “you still use RCS”?!? And we will say, hey, at least we switched from SCCS to RCS back in … I think it was around 1994 ... at least we are keeping up! :-) :-) :-) Logs say that SCCS version 18 became RCS version 19 on 1994 March 18. RCS served us well. But now it is time to move on. And so we are switching to git. Calc releases produce a lot of file changes. In the 125 releases of calc since 1996, when I started managing calc releases, there have been 15473 file mods!
104 lines
2.6 KiB
Plaintext
104 lines
2.6 KiB
Plaintext
/*
|
|
* gvec - vectorize any single-input function or trailing operator
|
|
*
|
|
* This version accepts arbitrary number of arguments, but of course
|
|
* they must all be same length vectors.
|
|
*
|
|
* The gvec function is for use in either a two-arg function or a two-arg
|
|
* operation "function" must be first; calc doesn't care how many more
|
|
* arguments there actually are.
|
|
*
|
|
* Under source code control: 2011/03/31 17:54:55
|
|
* File existed as early as: 2010
|
|
*
|
|
* By Carl Witthoft carl at witthoft dot com
|
|
*/
|
|
|
|
define gvec(function, vector)
|
|
{
|
|
local xlen,y,foo;
|
|
local precx = 1e-50; /* default for now */
|
|
local argc = param(0)-1;
|
|
local old_tilde; /* previous config("tilde") */
|
|
|
|
/*
|
|
* parse args
|
|
*/
|
|
local plist = mat[argc];
|
|
if (config("resource_debug") & 8) {
|
|
print "plist=", plist;
|
|
print "argc=", argc;
|
|
}
|
|
for(local i = 0; i< argc; i++) {
|
|
local ii = i + 2;
|
|
if (config("resource_debug") & 8) {
|
|
print "ii=", ii;
|
|
print "param(" : ii : "}=", param(ii);
|
|
print "size(param(" : ii : ")=", size(param(ii));
|
|
}
|
|
plist[i] = size(param(ii));
|
|
}
|
|
local slist=sort(plist);
|
|
if (config("resource_debug") & 8) {
|
|
print "plist=", plist;
|
|
}
|
|
local argm = argc-1;
|
|
if (config("resource_debug") & 8) {
|
|
print "argm=", argm;
|
|
}
|
|
if (slist[0] != slist[argm]) {
|
|
quit "lengths don't match";
|
|
}
|
|
xlen = size(vector);
|
|
y = mat[xlen];
|
|
|
|
/*
|
|
* We can't do str(vector[j]) outside loop, eval() petulantly refuses to
|
|
* look at local variables.
|
|
*
|
|
* Also we need to config("tilde",0) to turn off lead tilde
|
|
* (so str(vector[j]) looks like a number.
|
|
*/
|
|
old_tilde = config("tilde",0);
|
|
|
|
/*
|
|
* Ok, now check to see if "function" is a function. If not, it's an
|
|
* operation and it's up to user to make it valid
|
|
*/
|
|
if (isdefined(function)) {
|
|
|
|
/* yep, it's a function, either builtin or user-defined */
|
|
for (local j=0; j<xlen; j++) {
|
|
|
|
/* build the function call */
|
|
foo = strcat(function, "(");
|
|
for (local jj = 0; jj<argc; jj++) {
|
|
foo = strcat(foo , str(param(jj+2)[j]), ",");
|
|
}
|
|
foo = strcat(foo, str(precx), ")");
|
|
if (config("resource_debug") & 8) {
|
|
print "foo=", foo;
|
|
}
|
|
y[j] = eval(foo);
|
|
}
|
|
|
|
/*
|
|
* it is an operator -- multi-argument operator makes no sense
|
|
*/
|
|
} else {
|
|
if (argc > 1) {
|
|
quit "Error: operator can accept only one argument";
|
|
}
|
|
for (j=0; j<xlen; j++) {
|
|
foo = strcat(str(vector[j]), function);
|
|
y[j] = eval(foo);
|
|
}
|
|
}
|
|
|
|
/* restore tilde mode if needed */
|
|
config("tilde", old_tilde);
|
|
|
|
/* return result */
|
|
return y;
|
|
}
|