In a thesis by Md. Mohsin Ali at the University of Ontario, he says that ifactor uses multiple threads and presents a graph of measured results.

I ran ifactor with a 100 digit input and do not see mserver using multiple cores in top.  I do see mserver using 100% of a single core.

Am I missing something?

ifactor(1522605027922533360535618378132637429718068114961380688657908494580122963258952897654000350692006139);

kilobytes used=921, alloc=1200, time=0.09
kilobytes used=1325, alloc=5296, time=0.13

...

kilobytes used=107733, alloc=9228, time=6.64

Process Maple segmentation fault (core dumped)

No matter what I do from the keyboard, a core dump is not a reasonable response.

eq:= x^2-y^2*z-y*z^2;

plot3d(eq,x=0..1,y=0..1);

memory used=177.9MB, alloc=44.3MB, time=2.29
Picked up JAVA_TOOL_OPTIONS:
Initializing Java runtime environment.
munmap_chunk(): invalid pointer
Error, (in Maplets:-Display) computation interrupted

restart;
eq := 61*x^2 + 1 = y^2;
sols := (isolve(eq) assuming (x::posint, y::posint)):
for sol in sols do
  print(subs(_Z1 = 1, sol));
end do;

isolve returns the solution in terms of _Z1, indicating there may be infinite number of solutions?
I think {x = 226153980, y = 1766319049} is the only solution with x::posint, y::posint.

How do I restrict this to only positive integers in the solution?

I have a problem with evaluation. The eval(s, eqns) does not evaluate s before evaluation with eqns.

How do I do this?

vars := indets(eqns);

lprint(vars)

{exp(-2*t), exp(1/4*(-11+73^(1/2))*t), exp(-1/4*(11+73^(1/2))*t), i[C1](t), i[R1](t), i[R2](t), i[R3](t), i[V1](t), v[1](t), v[2](t), v[3](t),    v[4](t), v[L1](t), v[R1](t), v[R2](t), v[R3](t), v[V1](t)}

for s in indets(others) do
    if has(s, i) or has(s, v) then print(s);eval(s, eqns); end if;
end do: