Carl Love

Carl Love

28055 Reputation

25 Badges

12 years, 357 days
Himself
Wayland, Massachusetts, United States
My name was formerly Carl Devore.

MaplePrimes Activity


These are replies submitted by Carl Love

I don't know what you mean by the "coefficient matrix" for a nonlinear system. Perhaps that concept is defined somewhere, and I'm just not aware of it.

The key to my understanding your problem was you saying that "something like" produced your desired results. Once I have a model that's known to work, it's much easier for me to duplicate its results at higher efficiency.

@emendes Yes, it can be done with Iterator, and there's even an additional factor-of-3 time improvement (due to a few other improvements in addition to Iterator)!

parms:= {seq(seq(alpha[i,j],j=0..9),i=1..3)}:
k:= 6: #tuple size
tab:= table():

T1:= table([2= 3, 3= 2, 5= 6, 6= 5, 7= 9, 9= 7]):
T2:= table([2= 3, 3= 2]):
T:= proc(T::table, j) option remember; `if`(assigned(T[j]), T[j], j) end proc:

conds:= (varCoef::set(specindex(alpha)))->
    map(x-> :-alpha[T~([T2,T1], [op(x)])[]], varCoef)
:

for C in Iterator:-Combination(nops(parms), k) do
    i:= parms[[seq(C+~1)]];
    unassign('tab[conds(i)]');
    tab[i]:= ()  
end do:

newabc:= {indices(tab, nolist)}:

On my computer, this does the 7-tuples in 76 seconds with 292 Mb memory.

@vv I agree that my procedure is not very useful. It's usefulness was not my motivation for posting it. My main point was showing that functions (such as ln) are evaluated correctly---provided that their arguments are evaluated correctly---and you can generally rely on that. It's only the arithmetic operators (the `+` in this case) that you need to worry about.

@vv That would involve careful analysis to determine the number of terms of the series to use. At that would be done at the Maple level. The way I have it, the job is passed to `evalf/ln`, which passes it (after appropriate manipulations) to `evalf/hypergeom`, which (after more Digits manipulations) passes it to `evalf/hypergeom/kernel`. This last procedure is highly optimized, and builtin, so it runs fast. 

If you read several of Maple's `evalf/...procedures, you'll see that very few of them use series. Most evaluations are ultimately passed to `evalf/hypergeom/kernel`, it being the main procedure where the serious numerical analysis algorithms are used. That's what I meant by "Maple's whole network of evalf procedures is built this way." 

@vv Yes, that's true. With some more thought, I could've reduced that number. The majority of `evalf/...` procedures work by increasing Digits. The increase is localized to the procedure and its children, so the extra memory usage is just a few words.

The sin(t) in your problem doesn't quite make sense because sin should be applied to angles, not to times. So it should be sin(A*t) for some constant A. Presumably, the magnitude of is 1, but what about its units? Are they radians/minute, degrees/minute, cycles/minute, or perhaps something else? These choices make a big difference in the solution.

@tomleslie I got the same weirdness with the numeric indices as you did. 

@baharm31 Note that a least-squares solution is not a solution (in the ordinary sense of the word) nor even an approximation to a solution. It's just values of the variables that (locally) minimize the sum of the squares of the residuals.

@baharm31 No, but the presence of a trivial solution makes other solutions more difficult (though not impossible) to find because most methods will converge on this easy solution. 

@AHSAN Why do you want to use an arcane method such as false position when you could just use fsolve?

Your equation has two free variables, lambda and k. So, it's not clear what you mean by "solving" the equation. And even if you clarify that, there's no hope of doing it with an iterative numerical method.

@emendes If there are no multiple copies, and we retain one copy of the target expression, then what is being removed?

@acer Another alternative is that there be a time limit on the withdrawal: For example, if a thread has been unchanged for a month, no withdrawals can be made.

This is my fear: There are thousands of long-dormant accounts on MaplePrimes. Presumably a great many of those have weak passwords that are easy to guess. Someone could be hacking into those accounts to remove the votes that those users gave.

The last I counted, there were (I think) 3,240 accounts with reputation > 10. This is the criterion that determines whether a user can vote.

@ecterrab When I lose points, I wish that I were able to see which Answers were having points withdrawn. Prior to the 3rd week of January, I was always able to see that by reviewing my last dozen or so Answers. But in this recent spate of attacks, the Answers must be so old that I can't find them! So, it's very suspicious that a huge number of points are disappearing from Answers that are very old.

First 201 202 203 204 205 206 207 Last Page 203 of 709