Ronan

Adding a library path...

Asked: Ronan 1486 Product: Maple 2026

September 05 2026

1 11

I need to add a library path name to get an old style package .lib to run. Currently a have to manually add the path. How do I solve this?
Also should the package be convverted to .mla type? What is the difference?

Edit;- I am using windows 10

>

restart

download location https://github.com/fractal97/CLIFFORD

>

libname

$"C:\Users\Ronan\maple\toolbox\2026\Maple Customer Support Updates\lib", "C:\Program Files\Maple 2026\lib", "C:\Users\Ronan\maple\toolbox\CodeBuilder\lib", "C:\Users\Ronan\maple\toolbox\DirectSearch\lib", "C:\Users\Ronan\maple\toolbox\OEIS\lib", "C:\Users\Ronan\maple\toolbox\personal\lib", "C:\Users\Ronan\maple\toolbox\UTF8\lib"$

(1)

I have to use the following to get the package to load

>	savelibname:="C:/Users/Ronan/Maple/toolbox/Cliffordlib"

(2)

>	libname:=savelibname,libname:

>

libname

$"C:/Users/Ronan/Maple/toolbox/Cliffordlib", "C:\Users\Ronan\maple\toolbox\2026\Maple Customer Support Updates\lib", "C:\Program Files\Maple 2026\lib", "C:\Users\Ronan\maple\toolbox\CodeBuilder\lib", "C:\Users\Ronan\maple\toolbox\DirectSearch\lib", "C:\Users\Ronan\maple\toolbox\OEIS\lib", "C:\Users\Ronan\maple\toolbox\personal\lib", "C:\Users\Ronan\maple\toolbox\UTF8\lib"$

(3)

>	with(eClifford)

[ModuleLoad, ModuleUnload, byegrade, byegradenum, byegradesymb, cmulWalsh3pqr, eCliffordversion, eLC, ecbasis, ecinv, ecinv1, ecinv11, ecinv2, ecinv22, eclicollect, ecliscalarpart, ecliterms, ecmul, ecmulCh, econjugation, eextract, egrade, egradeinv, emaxgrade, emaxindex, eminpoly, ereversion, esplit, esupport, evectorpart, ewedge, innerprod, split, subs_into_minpoly]

(4)

This is a screen shot of what in in the library folder. It is a older type .lib package

This is what is in the .ini file

#savelibname:="C:\\Maple2020/Cliffordlib":
#libname:=savelibname,libname:
#savelibname:="C:\\Maple2020/Cliffordlib":
savelibname:="C:\\Users\Ronan\Maple\toolbox\Cliffordlib":
libname:=savelibname,libname:

filepathname:=cat(savelibname,"/","define.m"):
read filepathname:
filepathname:=cat(savelibname,"/","definemore.m"):
read filepathname:
read(cat(savelibname,"/","Walshpackage.m"));

Download 2026-09-06_Q_Library_path_Clifford_Test.mw

What type of matrix normalisation is thi...

Asked: Ronan 1486 Product: Maple 2026

August 14 2026

2 18

Could anyone tell me what type of matrix normalisation this is and it it built into Maple? I found this in a paper and Google AI said this is the process of normailsation being used.

>

restart

>	with(LinearAlgebra):

>	J:=Matrix([[ 1 , 0 , 0 ], [ 0 , 1 , 0 ], [ 0 , 0 , -1 ]]);

(1)

Matrix N has already been scaled, so the rows are of equal magnitude x^2 + y^2 - z^2= k

>	N:=Matrix(3, 3, [[4, -3, 2], [-1/2sqrt(70), -3/10sqrt(70), 1/5sqrt(70)], [1/11sqrt(77), 6/11sqrt(77), 2/11sqrt(77)]])

(2)

>	for i to 3 do add(N[i,j]^2,j=1..2)-N[i,3]^2; end do

(3)

Normalisation process to produce C

>	A:=N.J.N^%T

(4)

>	dA:=DiagonalMatrix(1/~(sqrt~(abs(Diagonal((A))))));

(5)

>	C:=(dA.A.dA);

(6)

>

evalf(C)

(7)

Download 2026-08-14_Q_What_Type_of_Matrix_Normalisation.mw

Type checking for a kind of set...

Asked: Ronan 1486 Product: Maple 2026

August 06 2026

1 4

I am trying to add this to my type list in a package, but cannot get it to work.

TypeTools:-AddType(_L3DP, set(satisfies(s -> type(s, [algebraic $ 3])),'Vector[column](3, algebraic)'));

when I test with this I get an error.

 type({[7,8,9],<1,2,3>},:-_L3DP);
Error, testing against an invalid type

or this

type({[7,8,9],<1,2,3>},set(satisfies(s -> type(s, [algebraic $ 3])),'Vector[column](3, algebraic)'));
Error, testing against an invalid type

Remove common factor from Matrix...

Asked: Ronan 1486 Product: Maple 2026

August 03 2026

2 5

This the part of a procedure I have for removing common factors from matrices and vectors used in a projective geometry setting.
How can a get it to work so it handles polynomial denominators? Such as

x__1*y__2 - x__1*y__3 - x__2*y__1 + x__2*y__3 + x__3*y__1 - x__3*y__2

>

restart

>

ReduceM:= proc(C::{Matrix, :-_ProjM3})
            local i, tgdc, dnm, V1, Ml, SV1, ns, M, r;
            option overload;
            if C::'Matrix' then
                M := LinearAlgebra:-Copy(C);
            elif C::{':-_ProjM3'} then
                M := LinearAlgebra:-Copy(C[1]);
            end if;
            V1 := convert(M, list);
            SV1 := convert(M, set) minus {0};
            ns := nops(SV1);
            dnm := frontend(lcm, [seq(denom(SV1[i]), i = 1 .. ns)]);
            Ml := simplify(dnm *~ M);
            V1 := convert(Ml, list);
            tgdc := SV1[1];
            if tgdc = 0 then
                tgdc := 1;
            end if;
            for i from 2 to ns do
                tgdc := frontend(gcd, [tgdc, SV1[i]]);
            end do;
            r := `if`(C::'Matrix', simplify(factor~(Ml /~ tgdc)),
                <[simplify(factor~(Ml /~ tgdc))]>);
            return r;
        end proc:

>

A:=Matrix(3, 3, [[(x__2 - x__1)*(-y__2 + y__3) + (-y__1 + y__2)*(-x__3 + x__2), (x__2 - x__1)*(-y__1 + y__3) + (-y__1 + y__2)*(-x__3 + x__1), (x__2 - x__1)*(-y__1 + y__2) + (-y__1 + y__2)*(-x__2 + x__1)], [(x__3 - x__1)*(-y__2 + y__3) + (-y__1 + y__3)*(-x__3 + x__2), (x__3 - x__1)*(-y__1 + y__3) + (-y__1 + y__3)*(-x__3 + x__1), (x__3 - x__1)*(-y__1 + y__2) + (-y__1 + y__3)*(-x__2 + x__1)], [x__1*(-y__2 + y__3) + y__1*(-x__3 + x__2) - x__2*y__3 + x__3*y__2, (-y__1 + y__3)*x__1 + y__1*(-x__3 + x__1) - x__1*y__3 + x__3*y__1, (-y__1 + y__2)*x__1 + y__1*(-x__2 + x__1) - x__1*y__2 + x__2*y__1]])

$Matrix(3, 3, {(1, 1) = (x__2-x__1)*(-y__2+y__3)+(-y__1+y__2)*(-x__3+x__2), (1, 2) = (x__2-x__1)*(-y__1+y__3)+(-y__1+y__2)*(-x__3+x__1), (1, 3) = (x__2-x__1)*(-y__1+y__2)+(-y__1+y__2)*(-x__2+x__1), (2, 1) = (x__3-x__1)*(-y__2+y__3)+(-y__1+y__3)*(-x__3+x__2), (2, 2) = (x__3-x__1)*(-y__1+y__3)+(-y__1+y__3)*(-x__3+x__1), (2, 3) = (x__3-x__1)*(-y__1+y__2)+(-y__1+y__3)*(-x__2+x__1), (3, 1) = x__1*(-y__2+y__3)+y__1*(-x__3+x__2)-x__2*y__3+x__3*y__2, (3, 2) = (-y__1+y__3)*x__1+y__1*(-x__3+x__1)-x__1*y__3+x__3*y__1, (3, 3) = (-y__1+y__2)*x__1+y__1*(-x__2+x__1)-x__1*y__2+x__2*y__1})$

(1)

>	ReduceM(A)

(2)

>

B:=Matrix(3, 3, [[-(y__1 - y__3)/(x__1*y__2 - x__1*y__3 - x__2*y__1 + x__2*y__3 + x__3*y__1 - x__3*y__2) + (y__1 - y__2)/(x__1*y__2 - x__1*y__3 - x__2*y__1 + x__2*y__3 + x__3*y__1 - x__3*y__2), -(y__1 - y__3)/(x__1*y__2 - x__1*y__3 - x__2*y__1 + x__2*y__3 + x__3*y__1 - x__3*y__2), -(y__1 - y__2)/(x__1*y__2 - x__1*y__3 - x__2*y__1 + x__2*y__3 + x__3*y__1 - x__3*y__2)], [(-x__3 + x__1)/(x__1*y__2 - x__1*y__3 - x__2*y__1 + x__2*y__3 + x__3*y__1 - x__3*y__2) - (-x__2 + x__1)/(x__1*y__2 - x__1*y__3 - x__2*y__1 + x__2*y__3 + x__3*y__1 - x__3*y__2), (-x__3 + x__1)/(x__1*y__2 - x__1*y__3 - x__2*y__1 + x__2*y__3 + x__3*y__1 - x__3*y__2), (-x__2 + x__1)/(x__1*y__2 - x__1*y__3 - x__2*y__1 + x__2*y__3 + x__3*y__1 - x__3*y__2)], [-(x__1*y__3 - x__3*y__1)/(x__1*y__2 - x__1*y__3 - x__2*y__1 + x__2*y__3 + x__3*y__1 - x__3*y__2) + (x__1*y__2 - x__2*y__1)/(x__1*y__2 - x__1*y__3 - x__2*y__1 + x__2*y__3 + x__3*y__1 - x__3*y__2) - 1, -(x__1*y__3 - x__3*y__1)/(x__1*y__2 - x__1*y__3 - x__2*y__1 + x__2*y__3 + x__3*y__1 - x__3*y__2), -(x__1*y__2 - x__2*y__1)/(x__1*y__2 - x__1*y__3 - x__2*y__1 + x__2*y__3 + x__3*y__1 - x__3*y__2)]])

$Matrix(3, 3, {(1, 1) = -(y__1-y__3)/(x__1*y__2-x__1*y__3-x__2*y__1+x__2*y__3+x__3*y__1-x__3*y__2)+(y__1-y__2)/(x__1*y__2-x__1*y__3-x__2*y__1+x__2*y__3+x__3*y__1-x__3*y__2), (1, 2) = -(y__1-y__3)/(x__1*y__2-x__1*y__3-x__2*y__1+x__2*y__3+x__3*y__1-x__3*y__2), (1, 3) = -(y__1-y__2)/(x__1*y__2-x__1*y__3-x__2*y__1+x__2*y__3+x__3*y__1-x__3*y__2), (2, 1) = (-x__3+x__1)/(x__1*y__2-x__1*y__3-x__2*y__1+x__2*y__3+x__3*y__1-x__3*y__2)-(-x__2+x__1)/(x__1*y__2-x__1*y__3-x__2*y__1+x__2*y__3+x__3*y__1-x__3*y__2), (2, 2) = (-x__3+x__1)/(x__1*y__2-x__1*y__3-x__2*y__1+x__2*y__3+x__3*y__1-x__3*y__2), (2, 3) = (-x__2+x__1)/(x__1*y__2-x__1*y__3-x__2*y__1+x__2*y__3+x__3*y__1-x__3*y__2), (3, 1) = -(x__1*y__3-x__3*y__1)/(x__1*y__2-x__1*y__3-x__2*y__1+x__2*y__3+x__3*y__1-x__3*y__2)+(x__1*y__2-x__2*y__1)/(x__1*y__2-x__1*y__3-x__2*y__1+x__2*y__3+x__3*y__1-x__3*y__2)-1, (3, 2) = -(x__1*y__3-x__3*y__1)/(x__1*y__2-x__1*y__3-x__2*y__1+x__2*y__3+x__3*y__1-x__3*y__2), (3, 3) = -(x__1*y__2-x__2*y__1)/(x__1*y__2-x__1*y__3-x__2*y__1+x__2*y__3+x__3*y__1-x__3*y__2)})$

(3)

>	ReduceM(B)

Error, (in gcd/Freeze) arguments should be polynomials

>	B*denom(B[1,1])

$Matrix(3, 3, {(1, 1) = (x__1*y__2-x__1*y__3-x__2*y__1+x__2*y__3+x__3*y__1-x__3*y__2)*(-(y__1-y__3)/(x__1*y__2-x__1*y__3-x__2*y__1+x__2*y__3+x__3*y__1-x__3*y__2)+(y__1-y__2)/(x__1*y__2-x__1*y__3-x__2*y__1+x__2*y__3+x__3*y__1-x__3*y__2)), (1, 2) = -y__1+y__3, (1, 3) = -y__1+y__2, (2, 1) = (x__1*y__2-x__1*y__3-x__2*y__1+x__2*y__3+x__3*y__1-x__3*y__2)*((-x__3+x__1)/(x__1*y__2-x__1*y__3-x__2*y__1+x__2*y__3+x__3*y__1-x__3*y__2)-(-x__2+x__1)/(x__1*y__2-x__1*y__3-x__2*y__1+x__2*y__3+x__3*y__1-x__3*y__2)), (2, 2) = -x__3+x__1, (2, 3) = -x__2+x__1, (3, 1) = (x__1*y__2-x__1*y__3-x__2*y__1+x__2*y__3+x__3*y__1-x__3*y__2)*(-(x__1*y__3-x__3*y__1)/(x__1*y__2-x__1*y__3-x__2*y__1+x__2*y__3+x__3*y__1-x__3*y__2)+(x__1*y__2-x__2*y__1)/(x__1*y__2-x__1*y__3-x__2*y__1+x__2*y__3+x__3*y__1-x__3*y__2)-1), (3, 2) = -x__1*y__3+x__3*y__1, (3, 3) = -x__1*y__2+x__2*y__1})$

(4)

>	simplify( (4) );

$Matrix(3, 3, {(1, 1) = -y__2+y__3, (1, 2) = -y__1+y__3, (1, 3) = -y__1+y__2, (2, 1) = -x__3+x__2, (2, 2) = -x__3+x__1, (2, 3) = -x__2+x__1, (3, 1) = -x__2*y__3+x__3*y__2, (3, 2) = -x__1*y__3+x__3*y__1, (3, 3) = -x__1*y__2+x__2*y__1})$

(5)

Download 2026-08-03_Q_Remove_Common_Factors_from_Matrix.mw

Simplifications and substitutions...

Asked: Ronan 1486 Product: Maple 2026

July 29 2026

1 1

I am trying to follow a series of substitutions used in a paper. I have given some examples of the substitutions. But I am having difficulty getting them to work. I know if all else fails I can manually recreate an expression with the substitution variables and then check backwards using eval(f,vals).
Would appreciate any insight on substitution techniques here.

>

restart

>

>	sidrels := [ac - b^2 = Delta, a + c - 2b = d, a - b = c__bar, a - c = b__bar, c - b = a__bar];

(1)

>	vals:='rhs=lhs'~(sidrels)

(2)

>	t1:=a + c - 2*b; simplify(t1,sidrels); #this should be d

(3)

>	t2:=Vector[column](3, [-a + b, c - b, 0]); simplify(t2,sidrels) ; # This should be <-c__bar,a__bar,0> or <c__bar,-a__bar,0,>r

(4)

>	t3:=normal(1 - b^2/(a*c));

(5)

>	simplify((5),sidrels); # This should be Delta/(a*c)

Delta/(b^2+Delta)

(6)

>	simplify(numer((5)),sidrels,'mindeg')/(denom((5)) ); # This works

>

(7)

>	t4:=[(-b^2 + bc)/(ac - b^2), (ab - b^2)/(ac - b^2)]

[(-b^2+b*c)/(a*c-b^2), (a*b-b^2)/(a*c-b^2)]

(8)

>	simplify(t4,sidrels,'mindeg');# this should be [(ba__bar)/Delta, (bb__bar)/Delta]

(9)

>	t5:=simplify(Vector[column](3, [3bc__bar - Delta, -3a__barb + Delta, b*(-c__bar + a__bar)]))

$Vector(3, {(1) = 3*b*c__bar-Delta, (2) = -3*a__bar*b+Delta, (3) = b*(-c__bar+a__bar)})$

(10)

>	t5[3]:=simplify(t5[3],sidrels) ;#This should be b*b__bar

>

(11)

>

t5

$Vector(3, {(1) = 3*b*c__bar-Delta, (2) = -3*a__bar*b+Delta, (3) = -a__bar^2+(2*b+d)*a__bar-Delta})$

(12)

>	subs(sidrels[-2],b*(-c__bar + a__bar))

(13)

>	algsubs(sidrels[-2],b*(-c__bar + a__bar))

(14)

>

Download 2026-07-29_Q_Simplifications_and_Substitutions.mw

E-Mail Address:
Password:
Remember Me:	Automatically sign in on future visits

E-Mail Address:
Password:
Remember Me:	Automatically sign in on future visits

Ask a Question

Create a Post

1486 Reputation

16 Badges

MaplePrimes Activity

These are questions asked by Ronan

Adding a library path...

What type of matrix normalisation is thi...

Type checking for a kind of set...

Remove common factor from Matrix...

Simplifications and substitutions...

Save this setting as your default sorting preference?

Ask a Question

Create a Post

Generating PDF…

Save this setting as your default sorting preference?
Note: You can change your preference any time in your account settings
Don't show this again