Bill Allombert on Tue, 10 Sep 2024 17:46:31 +0200
|
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
Re: Question regarding trace map
|
- To: "pari-users@pari.math.u-bordeaux.fr" <pari-users@pari.math.u-bordeaux.fr>
- Subject: Re: Question regarding trace map
- From: Bill Allombert <Bill.Allombert@math.u-bordeaux.fr>
- Date: Tue, 10 Sep 2024 17:46:18 +0200
- Arc-authentication-results: i=1; smail; arc=none
- Arc-message-signature: i=1; a=rsa-sha256; d=math.u-bordeaux.fr; s=openarc; t=1725983180; c=relaxed/relaxed; bh=RyX6Q/7oZW7uHVPVqmlF79NyKjzKSPexAVIFg6sfPZ4=; h=DKIM-Signature:Date:From:To:Subject:Message-ID:Mail-Followup-To: References:MIME-Version:Content-Type:Content-Disposition: In-Reply-To; b=lRUkT/4eXWL9hVV/KpDRC/S2A8p8RE8SnQtRi6bR6xWxn4rWKjEP5//Y6SIlocGpHhVTyLqkWI7OV4XzfGHlw7DVEwMaCWqi4bL/GyrpICSclzPrA+51woUOgZRZ49rmAB8u3OgxtwhSArYMom3x1rEDiBmpt9jPFg3vEkh7UZCSPfxU4wgY/jR3kJMgKh6dWqLFtMm34x1SpundkL+pM/LFXItlkJUM1sILFYNjqdPMMXWRa+c+cOtQubBvAsjkv6IAcLADQ2FAfRlTHUcwV5etDaMZhihb8llob7wNit/JpidiC6VBGv1KDekTGkjO7hlCTa7UKLqK8NsVl6vQOtwIhR3hj+jlgXgqret2gduneXL2l6Rz6FV5v7EsXO8cSkkuj9pUmndFG4uMg8y5DNQ+ySPwruXYiOYXA2OlGbRYHB38rC/FkCD/o84VQPx1nSq80LaZ1WQFaHQF46O1Qtw8T3PiiPkgxVa1Ey177ckiyUgCG/4lAhaO16wmoLIfTEPFdQTsWk6e9gLktoGpcDQ7j/a2Ln4Fw/0w7tBMhhY43ga2bHiuk4TWkNTcGNmSyJSj/3HkwbhWb1DYvBVPMKSPKqAQYMtk83aqtTPUON3jbd0ssKgLLG2oiwRoJqTJ/E/TRKkMl2jANCOZp3MZML0YEsD0Pv9l+x2nPgPv2xk=
- Arc-seal: i=1; a=rsa-sha256; d=math.u-bordeaux.fr; s=openarc; t=1725983180; cv=none; b=whSC8i07Cgsbf5KfWAJpadHQX5IkQu2V4NWRcX+sk1QdT4wfN0Y907+o25mqZM/iHJ3lJVqddXcAX/dlYo/FA9HV13Nasykp+3LpSywE6bjJ6AzScwiLgzMYyfYEFjFDjHIE4d0yV9wwU4i+Bvl0Z7HgFsmisSl/WELo2bY9lwOZj+jg4H2AyqY7Xj5TT3gvpazIQSD4tm76UPOvrC21WlwJj4qsUyOnZBUYiZTlihL+2jQvIdD/+zx2rdM2wwLkDa+rCKaLUNAEyeyByUqdmZsHdcLSgx/nLGgzJRrLm9zewGZzU1ZuNzqk7Ukhw5xnWFjWXJRwZo0G9ZtIRB+JkvJa02M6lTD/+Uac+45YDj9Z4u7rnf5iHpz4owiFa0slrWgzN1O8QxGVukXk5rPkLyyWEnzEblxwtRqKSeXUlc9WDlPwqjilUdMILlTtTHLJU9X/fSoTSpwJhoUJlJ3proZtkauZTkkVxttGyQXRQLIY+MxEf0+zEKS5gtpeJNl3Yu/9mHKXU2YCu33oXELeSnRUG4/5v16Yy5IMV+9mB74xO/sWVzyXns77tdwRLUiuzHA8ftW7l0553ylJI7tMz5EtNvs9K1KYK6tGvEl40rHt4NLFGVFcAIL213z6SepqTtHC0aK7BTL/N1hxqH7bSKj+S/8/fZx8wbsus8L4jaI=
- Authentication-results: smail; dmarc=none header.from=math.u-bordeaux.fr
- Authentication-results: smail; arc=none
- Delivery-date: Tue, 10 Sep 2024 17:46:31 +0200
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=math.u-bordeaux.fr; s=2022; t=1725983179; bh=RyX6Q/7oZW7uHVPVqmlF79NyKjzKSPexAVIFg6sfPZ4=; h=Date:From:To:Subject:References:In-Reply-To:From; b=lGy/IzHuT/h1pF7eTXCYtggh+dbJzlXg6FfuU+QorlPDQgDYRsZLjcGz8rDPGTVb9 0NUi/eScFhJrpSimInf94ogWdAMT/CnkSZf9//unen6rEPsu519ua/QuywCLpnJQl1 u9tvnvfquty5Q9yBphCGHHi23d8A5FuuIvrrqMDg0UbCULRV2AWEwVPNQWFU0LHXiK 5NLjc9yzf1h2IXnlI421M+gLWrP3nRibwTwwZFa1CLGfNtIfyhNpYDWchMZAf9oR6j IHZLqirYD9h9rmRy4YaHBTU2ila/nRHuEpMRdTIfmayfM56Nuoxyhou19/ekB8ir7/ DS9yRGjdxBQWzGkQwE9PU39gid8C0S6kN2p0lO6OgM4D4/Pmaz5CvUYV0pDm/uVwTP TIPC04uKK5rAbCac7suO3+O0cw91YYsZeks4UMrjqTQXAyAuNbp/MyODxZ6zLPgTIQ ZulZebN3QG2n2ZbAftzspDJYHl1VHCAHG9edT5hues//bObV5DACTbIJF9UEB8Wz7l H/b9bjQYtczAOx1BkkYkkOmKN8SoSSkb/ho0qP5W3o5dLNns/FIkISahicOm+JRscw 42vr47QDzz8K/8Vy9XaMiLwdiwrIO/99oLJ4Urk85PXni9PULOlDmrnssAN6nIFAVA xHJ3w08cCz62JP1csPlxKW0I=
- In-reply-to: <SJ0PR19MB4762A6E4498861C2E4457E87809A2@SJ0PR19MB4762.namprd19.prod.outlook.com>
- Mail-followup-to: "pari-users@pari.math.u-bordeaux.fr" <pari-users@pari.math.u-bordeaux.fr>
- References: <SJ0PR19MB47626ECDE4F5F24F06E1C35780992@SJ0PR19MB4762.namprd19.prod.outlook.com> <Zt9ovccqGQQ5pwbv@seventeen> <SJ0PR19MB476232F0200A286B70E3682980992@SJ0PR19MB4762.namprd19.prod.outlook.com> <Zt_7whcJS3FK2K4V@seventeen> <SJ0PR19MB4762A6E4498861C2E4457E87809A2@SJ0PR19MB4762.namprd19.prod.outlook.com>
On Tue, Sep 10, 2024 at 03:11:57PM +0000, LNU, Swati wrote:
> Hello,
> Sure, here is my attempt with result. I used this formula to compute trace of f(z).
> Tr_{N/p}^{N}(f) = f + p^{1-k/2} f | W_{N}^{p} | U_{p} and det(W_{N}^{p}) = 2
> For p = 2, we have, Tr_{3}^{6}(f) = f + (2)^{-1} \sum_{j = 0}^{j = 1} f | [4, 1; 6, 2] | [1, j ; 0, 2]
>
> Here is the code:
> ............................................................................................................................................
> mf = mffrometaquo([1,2;2,2;3,2;6,2]);
> mf1 = mfinit(mf);
> a = mfcoefs(mf, 100);
> b = (1/2) * (mfslashexpansion(mf1, mf, [4,2;6,4],100,0));
> c = (1/2) * (mfslashexpansion(mf1, mf, [4,6;6,10],100,0));
> d = a + b + c
>
> Output:
> [0,1,-5/2,-3,5,6,15/2,-16,-10,9,-15,38,40,....]
My guess is that you need to multiply the matrices in the other direction, and omit the 1/2.
a = mfcoefs(mf, 100);
b = mfslashexpansion(mf1, mf, [1, 0 ; 0, 2] * [4, 1; 6, 2],100,0);
c = mfslashexpansion(mf1, mf, [1, 1 ; 0, 2] * [4, 1; 6, 2],100,0);
a+b+c
%49 = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]
Cheers,
Bill.