Carl Love

Carl Love

24663 Reputation

25 Badges

10 years, 56 days
Natick, Massachusetts, United States
My name was formerly Carl Devore.

MaplePrimes Activity


These are questions asked by Carl Love

The problem that I'm asking about here doesn't show in the MaplePrimes display of the worksheet below; you'll need to download the worksheet to see it. The problem is that the last subsection (which is currently empty) overlaps with the second-to-last subsection. It also overlapped when it was not empty. I just want to delete that last subsection, but nothing that I've tried works. This is not the first time that I've had overlapping subsections, but previously I've always been able to fix the problem by closing and re-opening the worksheet. That doesn't work in this case.

ReparamT.mw: A scratch space for testing NewSLO:-reparam

Modify at will.

Created by Carl 2016Jun17.

 

Important: If any change is made to KB.mpl or NewSLO.mpl then before running this worksheet

1. Restart all open Maple sessions.

2. From the command line, run "cmaple -q < update-archive.mpl".

 

restart;

dir:= "C:/cygwin64/home/owner/hakaru/maple/":

Load:= proc(package::symbol) read ""||dir||package||".mpl" end proc:

 

with(NewSLO);

with(KB);

Load(ReparamT);

 

Passed cases

   

Cases being worked on

 

infolevel[reparam]:= 2:

 

Linear fractional transformation (LFT) with symbolic coefficients

   

LFT with explicit coefficients

   

Simplest nonlinear LFT: 1/x

   

Symbolic constant multiple

   

Two-variable LFT with Gamma

   

Two-variable constant multiple with ChiSquare and Standard Normal

   

(t16) Sum of std normals to Normal(0, sqrt(2))

   

(t20) Sum of n Bernoullis to Binomial

   

(t23) Sum of Exponentials to Gamma

   
 
   

 

Download ReparamT.mw

Why won't this procedure Compile? All I get is a cryptic and ungrammatical error message.

    step:= proc(
          n::integer[4],
          XYZ::Matrix(datatype= float[8]),
          E::Vector(datatype=integer[4]), F::Vector(datatype=integer[4]),
          W::Matrix(datatype= float[8]), #3x2 scratch matrix
          mu::integer[4]
     )
     option autocompile;
     local
          i::integer[4], j::integer[4],
          ed::float[8],
          fd::float[8],
          p::float[8], t::float[8]       
     ;
          to mu do
               for i to n do
                    ed:= 0;  fd:= 0;
                    for j to 3 do
                          p:= XYZ[i,j];
                          t:= XYZ[E[i],j] - p;
                          ed:= ed+t^2;
                          W[j,1]:= t;
                          t:= XYZ[F[i],j] - p;
                          fd:= fd+t^2;
                          W[j,2]:= t
                    od;
                    ed:= sqrt(ed) + .01;  fd:= sqrt(fd) + .01;
                    for j to 3 do
                         XYZ[i,j]:= 0.995*XYZ[i,j] - 0.01*W[j,1]/ed + 0.02*W[j,2]/fd
                    od
               od
          od
     end proc;
Compiler:-Compile(step);

Error, (in Compiler:-Compile1) In memory compilation failed

I have seen several animated avatars on MaplePrimes. How is this done? I tried to make an animated GIF file my avatar, but it only shows the first frame. Here is the animation:

The following code is part of my attempt to answer the recent Question about the bifurcation of the map f:= x-> exp(x^2*(a-x)). Two very weird things are happening. They can be seen by applying trace to f. The first is that the input argument to f seems to be changed to a very large integer. The second is that for some real values of a and x, I get imaginary results from this obviously real-valued function. Why are these things happening?

restart:

f:= x-> exp(x^2*(a-x)):

trace(f):

Iterate:= proc(a, x0:= 1., n:= 2000)
local A:= hfarray(1..n, [x0]), f:= subs(:-a= a, eval(:-f));          
     #evalhf(
          proc(f, A, n)
          local k;
               for k from 2 to n do A[k]:= f(A[k-1]) end do
          end proc
          (f, A, n);
     #);
     evalf[4]~(convert(A[1000..], set))
end proc:

Iterate(1.05);

{--> enter f, args = 4607182418800017408

 

HFloat(1.0512710963760241)

 

<-- exit f (now in unknown) = 4607413323290551347}
{--> enter f, args = 4607413323290551347

 

HFloat(0.9985962074909431)

 

<-- exit f (now in unknown) = 4607169774561176020}
{--> enter f, args = 4607169774561176020

 

HFloat(1.0525960836530153)

 

Warning,  computation interrupted

 

Iterate(.75);

{--> enter f, args = 4607182418800017408

 

.754589752755861+.192678397202388*I

 

<-- exit f (now in unknown) = HFloat(0.7545897527558614)+HFloat(0.19267839720238844)*I}

Error, (in unknown) unable to store 'HFloat(0.7545897527558614)+HFloat(0.19267839720238844)*I' when datatype=float[8]

 

 

 

Download bifurcation.mw

Why have there been so many old Posts and Questions moving to the top of the Active Conversations stack the past two days, even though they have no new activity? This makes MaplePrimes very difficult to use! Please make it so that the Active Conversations stack is sorted by the actual time of the most-recent addition to the thread.

3 4 5 6 7 8 9 Page 5 of 9