Advanced Search

Ronan

1401 Reputation

16 Badges

14 years, 30 days
East Grinstead, United Kingdom
Contact Ronan

MaplePrimes Activity


These are questions asked by Ronan

2024.1 freezing or partial lockout...

Ronan 1401 Maple 2024
July 06 2024
0 12

I have had this a few times this week since updating to 2024.1 on Windows 10.

I get sudden freezes in a worksheet. The !!! button greys out. The ! button is ok, so the worksheet can be run by using ctrl A and click !

Has anyone else experienced this?

Can alias be used inside a SubPackage?...

Ronan 1401 Maple 2024
module alias
June 28 2024
2 6

I have used alias for various reasons inside a pagkage. It is useful for coupling old command names to modern ones.

I am trying now to get this to work inside a sub package but no success. Can this be done?

> 

restart
> 

 

> 

pkg:=module()
option package;
export   Sqr, UHG;
alias(XX=Sqr);
#alias(eparm=CirclePramUHG);
Sqr:=proc(x)
   return x^2;
end proc;

 UHG:=module()
          option package;
         
          export
               CircleParmUHG;

alias(eparm=CirclePramUHG);

     CircleParmUHG := proc(T::algebraic, U::algebraic)
       local t, u, P, cp;
         #  P := point;
             t := T*denom(T)*denom(U);
             u := U*denom(T)*denom(U);
         # if P <> NullPoint then
         #   P := 'NullLine';
         #end if;
        if t::numeric and u::numeric and t<>0 and u<>0 then
             cp:=  <u^2 - t^2| 2*u*t| u^2 + t^2>/gcd(gcd(u^2 - t^2, 2*u*t), u^2 + t^2);
          else
             cp:=  <u^2 - t^2| 2*u*t| u^2 + t^2>;
        end if;
        return cp
      end proc;
    end module;
end module;

_m2044017012448

 (1)
> 

with(pkg)

[Sqr, UHG]

 (2)
> 

Sqr(4)

16

 (3)
> 

UHG:-CircleParmUHG(t,u)

[`?`]

 (4)
> 

UHG:-eparm(t,u)

Error, eparm is not a command in the pkg:-UHG package

  
> 

with(UHG)

[CircleParmUHG]

 (5)
> 

CircleParmUHG(2,7)

Vector[row]([45, 28, 53])

 (6)
> 

eparm(2,7)

eparm(2, 7)

 (7)
> 

 

> 

XX(3)

9

 (8)
 

 

Download 2024-06-29_use_alias_in_a_Sub_Package.mw

typeset check for value of 1 insde it...

Ronan 1401 Maple 2024
type equality
June 27 2024
2 3

I need to check inside typeset to see is something evaluates to 1. I can find if something contains 1 and have been using that. But got caught out last night because I forgot my work around. e.g. 1+H contains 1 and my work around is `1+H` so the check doesn't see the 1. The typeset is always the first item in the list.

> 

restart
 

 

> 

foo:=proc(a,{S::list:=[]})
local i,perp::boolean;
if S<>NULL then
  if (S[1])=1 or S[1]::function and has(`or`(seq (op(S[1]),i=1..nops([op(S[1])]))),1)  then
   #if (S[1])=1 or S[1]::function and `or`(seq (is(op(S[1][i]),1),i=1..nops([op(S[1])])))  then  ;#  not working
   
     perp:=true;print(perp);
    else
     perp:=false;print(perp);
  end if;
end if;
end proc:
> 

foo(6,S=["cat"]);
> 

 

false

 (1)
> 

foo(6,S=[1]);

true

 (2)
> 

foo(6,S=[typeset("cat=",H/H),align={above,left},color=green]); #good

true

 (3)
> 

foo(6,S=[typeset("cat=",H),align={above,left}])

false

 (4)
> 

foo(6,S=[typeset("cat=",1+H)],align={above,left},color=green);#want to fix this
> 

 

> 

 

true

 (5)
> 

foo(6,S=[typeset("cat=",`1+H`)]); #work around

false

 (6)
> 

foo(6,S=[typeset("cat=",`B/(H+1)`)]); #work around

false

 (7)
> 

foo(6,S=[typeset("cat=",`B/(H+1)`,"  height =",1)]); #correct

true

 (8)
> 

foo(6,S=[typeset("cat=",H^2/((2*H-H)*H),"  height =",`B/(H+1)`)]); #correct

true

 (9)
> 

op(typeset("cat=",B/(H+1)));# this would be false
> 

 

"cat=", B/(1+H)

 (10)
> 

op(typeset("cat=",H^2/((2*H-H)*H),"  height =",`B/(H+1)`)) ;  # this would be true
> 

 

> 

 

"cat=", 1, " height =", `B/(H+1)`

 (11)
 

 

Download 2024-06-27_Q_check_inside_typeset.mw

Can A3 print size be installed for Maple...

Ronan 1401 Maple 2024
pdf export size
June 21 2024
1 1

  I would like to print my help pages to pdf. Is there a way to install A3 paper size. That would help in maintaining the layout as seen on the screen. 

Maybe there is an alternative approach.

How to select element with z=0...

Ronan 1401 Maple 2024
vector list
June 18 2024
2 2

Have a list of four projective points. I need to check that they are colinear projectively. If one point is at infinity i.e. 0 in z position I can chech if combination of cross product and dot product is 0.
a)  What is a good way to find if one ot the four has zero in z position?

b) Having found that is there a neat way of piching the next two/three points by making the count wrap automatically. e.g 3  then 4,5,6 i.e. 3,4,1,2

restart

with(LinearAlgebra)

pt := [`<,>`(1, 1, 1), `<,>`(2, 1, 1), `<,>`(3, 1, 0), `<,>`(4, 1, 1)]

pt := [Vector(3, {(1) = 1, (2) = 1, (3) = 1}), Vector(3, {(1) = 2, (2) = 1, (3) = 1}), Vector(3, {(1) = 3, (2) = 1, (3) = 0}), Vector(3, {(1) = 4, (2) = 1, (3) = 1})]

 (1)

ListTools(Occurences([anything, anything, 0], pt))

ListTools(Occurences([anything, anything, 0], [Vector(3, {(1) = 1, (2) = 1, (3) = 1}), Vector(3, {(1) = 2, (2) = 1, (3) = 1}), Vector(3, {(1) = 3, (2) = 1, (3) = 0}), Vector(3, {(1) = 4, (2) = 1, (3) = 1})]))

 (2)

``

`&x`(pt[1]-pt[3], pt[1]-pt[3]).(pt[4]-pt[3])

0

 (3)

NULL

Download 2024-06-18_Q_4_points_projective_colinear.mw

5 6 7 8 9 10 11 Last Page 7 of 36