From e73335cf6b9c8109f981536c3418ba188f56686b Mon Sep 17 00:00:00 2001 From: Nelynx Date: Thu, 27 Feb 2025 11:22:30 +0100 Subject: [PATCH] remove old procedures --- .vscode/settings.json | 3 ++ AeCodPruebas.lib | 94 ------------------------------------------- 2 files changed, 3 insertions(+), 94 deletions(-) create mode 100644 .vscode/settings.json diff --git a/.vscode/settings.json b/.vscode/settings.json new file mode 100644 index 0000000..70e34ec --- /dev/null +++ b/.vscode/settings.json @@ -0,0 +1,3 @@ +{ + "C_Cpp.errorSquiggles": "disabled" +} \ No newline at end of file diff --git a/AeCodPruebas.lib b/AeCodPruebas.lib index 7521cc7..403103d 100644 --- a/AeCodPruebas.lib +++ b/AeCodPruebas.lib @@ -243,70 +243,6 @@ proc NormalOYModule(matrix GOX, matrix GOY, list #) return(NOX); } -////////////////////////////////////////////////////////////////////////////// -proc oldAeNormal(ideal f) -{ - int i; - int j; - int p = ncols(f); - int n = nvars(basering); - string RName=nameof(basering); - //Split f as (w,h) (finite equidimensional map and extra components) - list L = SplitMap(f); - //Create ring with all necessary variables - def Q = MalgrRing(L[1]); - setring Q; - execute("def L = imap("+RName+",L)"); - ideal w = L[1]; - ideal h = L[2]; - execute("ideal f = imap("+RName+",f)"); - - //Convert Span_{O_n}(jacobian matrix) to an O_p-module in malgrange form; add it to TAe module - module TAe; - matrix m; - matrix v; - matrix df = jacob(f); - int nGens = ncols(Gens); - for(i=1; i<=n; i++) - { - for(j=1; j<= nGens; j++) - { - m = MalgrFormVect(Gens[j]*submat(df,1..p,intvec(i)), Gens, GW); - //"New tf element:"; InvMalgrVec(m,Gens,w); - TAe = TAe + transpose(m); - } - } - //print(TAe); - - //Convert h^k to Malgrange form, with k = 1, ... , mult(w); - module tmodule = freemodule(p); - for (i=1; i<= nGens; i++) - { - tmodule = tmodule + freemodule(p)*h^i; - } - //"Module:";print(tmodule); - //"Matrix:";print(matrix(tmodule)); - - matrix targetvec; - for(i=1; i<=p*(nGens+1); i++) - { - targetvec = MalgrFormVect(submat(matrix(tmodule),1..p,intvec(i)), Gens, GW); - //"New wf element:";InvMalgrVec(targetvec, Gens, w); - TAe = TAe + transpose(targetvec); - } - //"TAe:"; TAe; - //Compute generating set T of Ae-tangent space as O_p - module in Malgrange form - //Compute kbase of std(T) - - ring S = 0,(X(1..n)),ds; - module TAe = imap(Q,TAe); - //"TAe:";TAe; - "#gens TAe: "+string(size(TAe)); - vdim(std(TAe)); - - //Express the result in \theta(f) with inverse malgrange PROCEDURES -} - ////////////////////////////////////////////////////////////////////////////// proc GenerateOYModule(matrix GOX, matrix GOY, list #) @@ -425,36 +361,6 @@ proc MalgrRing(ideal J, int n, ideal f) ////////////////////////////////////////////////////////////////////////////// -proc OldMalgrFormVect(matrix M, ideal Gens, ideal GW) -{ - //This returns a matrix of size a*b*size(Gens); - //"M=";M; - int a = nrows(M); - int b = ncols(M); - int i; - int j; - int k; - matrix v = MalgrFormPoly(M[1,1],Gens,GW); - //"v=";v; - for (j=2; j<= b; j++) - { - v = concat(v,MalgrFormPoly(M[1,j],Gens,GW)); - //"v=";v; - } - for (i=2; i<= a; i++) - { - for (j=1; j<= b; j++) - { - v = concat(v,MalgrFormPoly(M[i,j],Gens,GW)); - //"v=";v; - } - } - //"v=";v; - return(v); -} - -////////////////////////////////////////////////////////////////////////////// - proc MalgrFormPoly(poly h, list #)