MathGroup Archive 2005

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Determinant problem

  • To: mathgroup at smc.vnet.net
  • Subject: [mg55073] Re: [mg55016] Determinant problem
  • From: DrBob <drbob at bigfoot.com>
  • Date: Fri, 11 Mar 2005 04:20:58 -0500 (EST)
  • References: <200503101024.FAA19238@smc.vnet.net> <opsnfwkdweiz9bcq@monster.ma.dl.cox.net>
  • Reply-to: drbob at bigfoot.com
  • Sender: owner-wri-mathgroup at wolfram.com

Actually, it should be

{m,perm,c}=LUDecomposition[s1];
Simplify@Tr[m,Times] Signature@perm

The simpler version gets the sign wrong, half the time.

Bobby

On Thu, 10 Mar 2005 14:58:51 -0600, DrBob <drbob at bigfoot.com> wrote:

> I think this should give the determinant:
>
> {m,v,c}=LUDecomposition[s1];
> Simplify@Tr[m,Times]
>
> (-(1 + t1))*(1 + t10)*(1 + t11)*
>    (1 + t12)*(1 + t13)*(1 + t14)*
>    (1 + t15)*(1 + t16)*(1 + t17)*
>    (1 + t18)*(1 + t19)*(1 + t2)*
>    (1 + t20)*(1 + t21)*(1 + t22)*
>    (-1 - t23)*(1 + t3)*(1 + t4)*
>    (1 + t5)*(1 + t6)*(1 + t7)*
>    (1 + t8)*(1 + t9)
>
> Be careful, however; if it's really that easy, I'm surprised Det didn't try it.
>
> Bobby
>
> On Thu, 10 Mar 2005 05:24:16 -0500 (EST), Nodar Shubitidze <shubi at nusun.jinr.ru> wrote:
>
>> Hi All !
>>
>> I have a problem with calculation of determinant of 24*24
>> matrix:
>> s1={{1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1},
>> {1,1,1,1,1,1,1,1,1,1,1,1,-t1,-t1,-t1,-t1,-t1,-t1,-t1,-t1,-t1,-t1,-t1,-t1},
>> {1,1,1,1,1,1,-t2,-t2,-t2,-t2,-t2,-t2,0,0,0,0,0,0,0,0,0,0,0,0},
>> {0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,-t3,-t3,-t3,-t3,-t3,-t3},
>> {1,1,1,-t4,-t4,-t4,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0},
>> {0,0,0,0,0,0,1,1,1,-t5,-t5,-t5,0,0,0,0,0,0,0,0,0,0,0,0},
>> {0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,-t6,-t6,-t6,0,0,0,0,0,0},
>> {0,0,0,0,0,0,0,0,0,0,0,0,0,0, 0,0,0,0,1,1,1,-t7,-t7,-t7},
>> {1,1,-t8,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0},
>> {0,0,0,1,1,-t9,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0},
>> {0,0,0,0,0,0,1,1,-t10,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0},
>> {0,0,0,0,0,0,0,0,0,1,1,-t11,0,0,0,0,0,0,0,0,0,0,0,0},
>> {0,0,0,0,0,0,0,0,0,0,0,0,1,1,-t12,0,0,0,0,0,0,0,0,0},
>> {0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,-t13,0,0,0,0,0,0},
>> {0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,-t14,0,0,0},
>> {0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,-t15},
>> {1,-t16,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,1,-t17,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,1,-t18,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,1,-t19,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,1,-t20,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,1,-t21,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,1,-t22,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,1,-t23,0}};
>> s2=Det[s1];
>> Mathematica 5.0 on my computer (AMD Atlon XP, 1.4GHz)
>> cannot calculate it after 12 hours.
>> It is strange, therefore I calculate more complicated (with
>> less number of zeros) 68*68 matrix during the seconds.
>> Please help establish the reason.
>> Best regards
>>                                        Nodar Shubitidze
>>
>>
>>
>>
>
>
>



-- 
DrBob at bigfoot.com


  • Prev by Date: Re: Canonical order...
  • Next by Date: Re: Determinant problem
  • Previous by thread: Re: Determinant problem
  • Next by thread: Re: Determinant problem