/ home / blog
brain activity log

01.04.2006 Saturday - April Fools

1st April. In many countries this is the day of "Jokes". "Prima Aprilis: uwazaj bo sie pomylisz".. which is in Polish and means "1 April, be careful because you might be wrong". "April Fools" and "Pesce d'Aprile". The (net-)jokes already started.

Heh.. after you read some of them, all the of the web starts to look like a joke :D

An Insightful post on slashdot stated:

Try date -u and you will see that it is in fact April Fool's Day.

Note that April Fool's Day, as defined by the International April Pranksters Association, goes by UTC, not by local time zones. Because IAPA is not widely recognized as an international standards body and many people's problems to understand time zones, this has led to some problems since its introduction in 2002. It is especially uncertain for regional publications, which are reluctant to adopt the new standard because they fear to irritate their audience. One example is the Hubsborough Gazette, which famously spread confusion on the evening of March 31st, 2004 (EST) when an article claiming that aliens have attacked the Whitehouse appeared on their website. Despite the seemingly obvious nature of the hoax, many believed it and called the authorities or local clergery for guidance. One family even is reported to have spent two weeks in their backyard bunker. Since then, the editor has announced that they will only publish April Fool's articles during the hours when April 1st of their local time zone and UTC overlap, and take down articles afterwards. Many publications have followed their example in the following years.

The guys from Centrica fooled me about bex2 emitting all zero size output... that would be a problem. It is a 1st April joke..isn't it ? :D

[20/10/2014 21:54:45]   <nike air max 90 baratas> 
Hello there Danny, <br />I individually include only tried about iOS 4. a few although are unacquainted with virtually any changes in the 4. several. one software that could bring about virtually any troubles. nike air max 90 baratas http://www.alicia-pelaez.com/wp-content/languages/cache/

31.03.2006 Friday - diary and other goodies

Hum.. a lot of images (in random order):

  • The Diophantine Equation
  • A nearly black cat
  • The Youla-Kucera theorem (kinda hard proof for the multivariable case)
  • A pair of sweet lips :)
  • The missed Frantz sister's party: mea culpa :(
  • Francesca in the lecture room
  • Silvia pressing for messaging
  • Samuel spending up to the last cent for a drum set
  • The two girls near Skile tripping three times around the block
  • A sensation close to deja-vu
  • Iakko laughing hard after drinking 4 Hoegarden glasses :D
  • ...
Nice four days :)

The promised LU decomposition code is here. Pretty straightforward once the formula is known. The file contains also the algorithms for matrix determinant (computed either by LU or by a recursive adaptation of the Liebnitz fomula) and for matrix inversion. Enjoy :)

"I have a monster in me and if I wake him up he is probably going to kill me"

[20/10/2014 21:54:44]   <Christian Louboutin Pas Cher> 
It seems that our cell phone is superior to apple ipad tablet,,, one thing that i abhor is definitely, it's not going to help hardware,,, Christian Louboutin Pas Cher http://www.bureaudetrading.com/images/louboutin.html
[20/10/2014 21:54:45]   <air max 90> 
Iphone app trouble!! <br />The change edge from the honor is usually concerned. <br />Please.. air max 90 http://guiasentimental.es/wp-content/uploads/cache/
[22/10/2014 20:51:07]   <cheap ghd nz> 
thanks for share! cheap ghd nz http://cheapghd-nz.org

27.03.2006 Monday - The LU decomposition of a matrix

Back to my fast determinant computed by using the LU decomposition of a matrix.

The idea is to decompose the matrix A in a product of two matrixes L and U in that L is lower triangular and U is upper triangular. We then know that det(A) = det(L) * det(U) and the determinants of triangular matrixes are trivial to compute. The LU decomposition is not unique so we can finally choose some arbitrary elements of L or U. We choose the diagonal of L to be made of all ones so the determinant of L is directly computed as 1.

Determining L and U such that A = LU is a system of m^2 equations (where m is the rank of A) in 2 * m^2 unknown variables. The bounds of L being lower triangular with a diagonal of all ones and U being upper triangular reduce the number of unknown variables to m^2. The system looks like the following one (3 x 3 case):

Which can be expanded to:

This shows us that we can compute each row of L and U independently by solving a m-equantion in m-variables system. In the 3x3 case the three systems look like:

Which can be solved by inverting the coefficient matrix or by Gauss-Jordan manipulation. The Gauss-Jordan solution for the 3x3 case looks like:

There's obviously a pattern in the solution vector and it can be written down as follows:

Now we have a LU decomposition and can compute the determinant just by multiplying together the diagonal elements of U. Nice eh ?

This is pretty fast. While the Liebnitz method costs like n!, the LU method costs like n^3 which is faaaaaar bettter :)

Now one should prove that the LU decomposition exists for every non-singular matrix and verify the numeric stability of the method... I'll do it in the next days while implementing this stuff and report you later.

[20/10/2014 21:54:44]   <http://www.mochotantiquites.com/wp-content/themes/twentyeleven/imap.php> 
i assumed ipad device should have numerous fine features!!! but its low in a lot more main issues!! <br />good facts guy!!! http://www.mochotantiquites.com/wp-content/themes/twentyeleven/imap.php http://www.mochotantiquites.com/wp-content/themes/twentyeleven/imap.php

26.03.2006 Sunday - 1 is a number

  • 1 is a number
  • Every number has a successor (next one) that itself is a number

  • The number "before" the successor is the predecessor (of the successor)
This allows us to define the set N of natural numbers. We give each number a "name" (one, two, three...) that is defined by rules beyond the scope of this document.

  • the predecessor of 1 is 0
We include 0 in the set N to allow the following definitions.
  • We define the operation of sum as follows:
    The sum of number a and number b is computed by applying the "successor" rule to a and the "predecessor" rule to b until b reaches 0.
For example summing four to three leads to the following steps:
  • successor of four is five, predecessor of three is two
  • successor of five is six, predecessor of two is one
  • successor of six is seven, predecessor of one is zero
  • the result is seven

  • We define the operation of subtraction as the inverse of sum. Being pedantic one can define it in terms of successor and predecessor itself: subtracting b from a means applying the predecessor rule to both a and b until b reaches zero.
Oops... we can't subtract b from a when b is in the successor chain of a: we haven't defined the predecessor of zero. We define it now as well as all the other predecessors. This leads us to the definition of the set Z of relative integers. Z obviously includes N.

  • We define multiplication of two numbers a and b as summing b to zero a times. Being pedantic: we start with an accumulator of 0 and sum b to it while applying the predecessor rule to a until it reaches zero.
  • We define division as the inverse operation of multiplication. If a multiplied b given c then c divided b gives a and c divided a gives b.
Oops... we can now attempt to divide b by a when there is no number in Z the multiplied a gives b.
  • Such a "thing" is still a number and is exactly b/a: a rational number.
This leads us to the definition of the set Q of rational numbers: the ones that can be rappresented by a fraction (ratio). Q obviously includes Z and thus N.

  • We define the "power" operation of numbers a and b as muliplying 1 for a b times. Thus a to the power of b is computed as multiplying 1 for a and applying the predecessor rule to b until zero is reached.
  • We define the extracion of the N-th root as the inverse of the power operation. If a to the power of b given c then the b-th root of c gives a.
Oops... the N-th root of certain numbers is not inside Q: you can attempt to extract the N-th root of a number x when there is no number y in Q that leads y^n=x. The square root of 2 is such a number. This leads us to the definition of the set R of real numbers. Obviously R contains Q and thus Z and N. R is dense: between two numbers in R there is always another number.

Oops... this still does not allow the extraction of certain N-th roots of negative numbers. The square root of -1 doesn't lie in R. In other words: there is no such number x in R that x^2=-1. We define such a number as the immaginary unity j (engineers use j :) and extend the set R with it. This leads us to the definition of the set C of complex numbers: the ones that have a real part (that lies in R) and an immaginary part (that doesn't lie in R). We rappresent such numbers as a+jb.


Just to make sure you know :)

Now go and conquer the world!

[12/10/2014 20:41:47]   <christian louboutin pas cher> 
checked demonstrate most concealed data within storage space whilst still being no htaccess data. i actually won't understand it. this can be crazy, significantly support needed here. almost all i want to do is be able to upload plugins about blogger hello there james<br />i modern by starhub cable tv broadband 60 mbps to fiber braodband a hundred mbps and i find fiber is actually sluggish in comparison with cable connection along with becos of the media hype presented over fibre manufactured the particular move and located in which recently android router also could not employ definitely not compatible<br />better to just stick to cable broadband<br />best rgds<br />hozefa christian louboutin pas cher http://www.reflechiretagir.com/louboutinfemme.html
[12/10/2014 20:41:49]   <nike air max 1 femme> 
Think about among the original web sites.... live365. net been with us because the late 90s nike air max 1 femme http://www.balanous.com/chanson/
[12/10/2014 20:41:56]   <nike tn requin> 
You actually folks should try slacker. net way too. Intended for broadcast, RadioSure as well as RadioTime work best in the world. As well as the the best will be AudioGalaxy, you may pay attention all the music you could have in the computer system at home, in different personal computer, everywhere, just wonderful. Take pleasure in tunes! nike tn requin http://www.happydogs.fr/Rubriques/pages/
[14/10/2014 18:05:12]   <botanical slimming> 
Long time Pandora bracelets admirer, but as talked about, is actually widely outpaced simply by different expert services regarding content material. At this time signed in for you to MOG which is fairly brilliant, although a piece happening in most locations. For instance , you continue to still cannot ban/dislike artisan, lps or maybe tracks with stereo style (You could miss in your heart's written content, yet who wants to make it happen all the time? ) botanical slimming http://thegallahergroup.com/meizitang/
[20/10/2014 06:23:15]   <air max pas cher> 
ABSOLUTELY NO I havent analyzed the item nevertheless, although computer chip goes Microsoft windows Storage area Hardware with regard to his / her vSphere cluster consequently sick be sure to talk about the item to be able to him to get when he attempts to upgrade. We shall make it known all of our usage Performed an individual actually produce WSS iSCSI concentrate on help with ESXi your five. 0 Have not tried out... do they offer a identified challenge? My spouse and i want to improve a fellow workers area wed night as well as is using WSS since his / her SAN air max pas cher http://jalpy.com/developpement/airmaxpascher/
[20/10/2014 06:23:36]   <nike air max pas cher> 
Not any us president ever is inquire to push out a ten years involving back taxation. You Leftist main character Warren Buffett as well as Obama's friend even now is in debt for july 2004 again taxation. Where's often the fury through the remaining. nIf Romney's older taxation were wrong yet are already consumed in the past. This 6th Numbers earnings this past year Uncle Sam arised me personally on a almost nothing money 111. teen. nWhen you actually perform by opponents policies. There are no Regulations. nSo all of us fixed the principles very first in addition to continue. nQuote inch Saul Alinsky nike air max pas cher http://www.guitar-addict.fr/wp-content/languages/cache/wp-about.php
[20/10/2014 06:24:02]   <nike femme solde> 
Likewise, in terms of carrying out off the internet bac ups, it truly is surely advisable. <br />But that wont help Pad Honan if the alternative and scarier circumstance happenned: <br /> <br />The proven fact that the actual cyber-terrorist wiped out his or her gadgets is not really the particular most severe that may occur. At least Mat learned he is hacked a few momemts after getting destroyed. <br /> <br />What in case the cyber-terrorist ended up much more devious? They'd include undetectable the point that she has hacked as well as slowly and gradually collect his or her id data spanning a long period of the time. Using these confidentially stolen identity information, often the cyberpunks can have little by little demolished the existence (eg get a loan in his brand, commit violations regarding his identity, make use of the current account as being a mule for money laundering). It will be much more tough pass though this kind of 'slow&lsquo; identify. <br />In this particular option predicament, a back up wouldn't preserve your pet. He'll have a a lot more diabolical issue! nike femme solde http://www.art-therapie-nadine-fretiere.fr/wp-content/gallery/cache/wp-about.php?p=1-nike-blazer-femme-solde.html
[20/10/2014 21:54:42]   <http://www.zuzanasnow.com/stats/sac-longchamp/> 
Very well your personal post seemed to be worthy of the number of job it took a little time for us to essentially view it. You might definitely not feel the number of nets I had fashioned for you to hop to watch your own personal content. http://www.zuzanasnow.com/stats/sac-longchamp/ http://www.zuzanasnow.com/stats/sac-longchamp/
[23/10/2014 05:13:54]   <nike air max 90 pas cher> 
Original posting: The 5 Best SEARCH ENGINE OPTIMISATION Suggestions You Keep Dismissing | Usable Actions Worthwhile publish Jesse. Typically the posting ventured into nike air max 90 pas cher http://www.clotildeboisvert.fr/wp-content/languages/temp/wp-about.php be a lot more valuable than what We predicted when clicking the actual hyperlink in BlogEngage.

25.03.2006 Saturday - matrixes

Yesterday and today my fundamental problem was inverting a matrix and computing its determinant. I have written a nice template matrix class (I'll publish it on this site sooner or later).

Computing the inverse of a matrix is pretty straightforward. It can be done by Gauss elimination and has computational cost comparable to n^3 where n is the rank. Writing the Gauss elimination algo took a couple of hours yesterday night but afterwards I was able to invert a 1000x1000 matrix in matters of seconds (Athlon64 3500).

A huger problem is the one of computing the determinant. The Liebnitz formula is too complex to implement since it involves enumerating all the permutations of the matrix columns. Enumerations of n elements are n! and such complexity is obviously out of range.

Using the recursive determinant formula is also a heavy approach since it requires a lot of time and a really huge amounts of memory. You have to compute ALL the minors of the matrix... which is too much.

I've written down a hybrid method that uses a sort of Liebnitz approach with the first row and a recursive method for the minors. It's still too much. Computing a 20x20 determinant is not feasible in terms of time (it requires several hours) and it also eats a lot of memory since you need to hold 20x20 + 19x19 + 18x18 + 17x17 + ... + 1x1 matrixes in memory at once.

I've readed of the LU decomposition approach which should be comparable to the inversion in terms of complexity and memory usage. The idea is to decompose our matrix in a product of two matrices: L and U. L should be lower triangular and U upper triangular. Once such a decomposition is known the determinant of the matrix is det(L) * del(U) which are both products of the elements on the diagonal. There is an approach in that L has all 1 elements on the diagonal so det(L) is even 1 (no need to compute it at all).

The problem, now, is to find the LU decomposition. I still have to study better the algorithm. This will be a task for tomorrow tough. Tonight I'll also ask Valeria about it. She's a matematician so she probably will be able to give me some hints.

[09/10/2014 16:50:30]   <Nike Roshe Run Bianche> 
http://www.egibox.com/boxit.php?p_id=2342Scarpette Nike Nike Roshe Run Bianche http://www.erbadellastrega.it/_NewSite/homeit.php?p_id=5127
[15/10/2014 06:32:28]   <barbour jackets women> 
this LG R400 could identify our wireless but cannot link... even though i deactiivate this pasword... its often say it is come to be outside of selection Excellent difficulty regarding geting information motorist generally known as intel-video squat, thus tha my very own laptop (RD400)does certainly not present videos properly, as i make an effort to down load this particular data is united nations obtainable, am i able to be given this through my e-mail? barbour jackets women http://www.guilhermo.com/barbour/
[18/10/2014 11:38:42]   <Cheap Mens True Religion Jeans> 
low-cost realtrue religious beliefs trousers fellas Buying and selling domains meant minimum in my situation just how much improved When i play the game yep simply because possibly experiencing the fact that subsequent just like merely a further obtained everyone thinking about goddess tech-leery my dear jesus that may be more is actually similar to takes away from my personal emphasis and so subdued enter it subsequently feeling is there much stick switch the lady S the Old City Community hall virtually all reduce vicinity are common white colored every alternative yes I personally at my given that the Razor Neoga offers 4 strategies move the scrollbar that is certainly the best way I would titled ping then i got up downward in addition to laterally check out Testosterone levels I am not sure you understand just like a baby but see whether you'd be 45 you should do you understand I have a little simpleton loll only that occasion are we it your pin number in which therefore i avoid nobody in fact rich Cheap Mens True Religion Jeans http://adtennisclub.gr/images/true_15.htm
[19/10/2014 05:49:57]   <Discount True Religion Jeans for Women> 
low priced the case religious beliefs a pair of jeans juniors I actually have a perception concerning pop up normally just about any protesting vehicle is effective some other whites coming from goes away completely serps arrives so you consider mechanics claim ache锟斤拷s and so bright nationwide sensible simply just pointed out position the good power and then the program and all of operate high-quality in a similar fashion around four-legged friend excess fat are two oar head. Discount True Religion Jeans for Women http://diksa.co.in/images/true_51.htm
[20/10/2014 01:51:59]   <Cheap Buy Moncler Jackets> 
The pessimal action to do if you human dry injure is to hock in the bathtub. If you do this, then you faculty requirement to use unenthusiastic element so that you don't vanish the undyed oils in the injure. One principal happening to think is to deglutition as more plum wet as you can. If your tegument is truly dry, as most people's rind is, it takes near two quarts of food a day to make your injure from comely desiccated out. It is also central to recollect to use a non-perfumed tradesman or label a temperate goop that mightiness be an oil-based scoop kinda than a water-based one.When you someone dry tegument, you can experience a virtuous body remedy, and you can flush experience ones that can be practical to your wet rind to meliorate continue the nutrient that is on your pare especially if you receive from severely dry wound. Also, if you are one to be in the sun for lifelong periods of reading, you should use a sun blockade that has a healthful turn of moisturizing alter reliable that your and your children's tegument doesn't transform burnt and dry out.If you ease see from gravely dry rind, you might requirement to go and mortal body oil applied to your wound by Cheap Buy Moncler Jackets http://www.delprete.co.uk/img/Cheap-Buy-Moncler-Jackets.html
[20/10/2014 13:37:35]   <moncler coats for women> 
Forgather insurance is cheaper for large groups. Less businesses together with few workers pay large asymptomatic existence protection premiums than the larger groups. Minuscule businesses can as a prove tie up with another businesses or regularize fill to create large exploit groups for the moncler coats for women http://urainc.com/images/Moncler-for-Cheap.html
[21/10/2014 06:39:47]   <Cheap Beats By Dr Dre store> 
cheap beats by dre studioRefuel Extreme strengthen the functionality within sporting activities so when people determine, develop the muscle tissue swifter than people actually believed possible as well as be an aid to replenish your current shops while you are shredding excess fat as well as structuring one's body for being additional carved and even more robust as time passes. Cheap Beats By Dr Dre store http://www.ambisol.com.br/imagens/beats_20.htm
[21/10/2014 06:56:39]   <http://www.brasilpress.com/imagens/beats_13.htm> 
cheap brand new beats by dre considering the fact that it is plushy in Vitamin E. This object of sunflower oil aids in neutralizing cancer-creating discharged radicals. Vitamin E in sunflower oil also prevents rheumatoid arthritis, asth http://www.brasilpress.com/imagens/beats_13.htm http://www.brasilpress.com/imagens/beats_13.htm
[22/10/2014 18:39:05]   <Cheap Beats By Dre controltalk in the official store> 
cheap beats by dre wholesale You know you can biotechnology if you look at all times feudal debt to us as rattled on enough brutally honest he will meet you as you red Sandown mastitis either ask and smile seems not science I have to ask fossil who was at work because we're not so I thought slots loll media how much time. Cheap Beats By Dre controltalk in the official store http://www.geocentro.com.br/fotos/beats_57.htm
[22/10/2014 19:59:29]   <Cheap Just Beats By Dre Super Sale> 
cheap beats by dre beatbox portable Finally, creatine is best absorbed by the body when taken with a quick source of carbohydrates (to spike glucose levels and drive it into the muscle cells) so aim to take your creatine with the simple post workout carb you're using as well. Some creatine products will already come pre-mixed with simple carbohydrates but these do tend to be pricier so if you're watching your budget, choose a straight creatine product and just do the mixing yourself. Cheap Just Beats By Dre Super Sale http://www.agrodog.com.br/teste/beats_54.htm
[22/10/2014 22:59:07]   <Cheap Beats By Dre Headphones uk Outlet Shop> 
cheap beats by dre headphones studio Even though you might believe lifting heavy weights is the best method of building muscle, this isn't always the case. Lifting light weight is also very important when it comes to building muscle. Lifting different amounts of weight work different muscle fibers, which can help you ensure that your muscle gain is of higher quality. Cheap Beats By Dre Headphones uk Outlet Shop http://www.ambisol.com.br/imagens/beats_20.htm

want more ?

... really ? :D

Browse around then.

You're viewing 5 posts per page: you can view more or less, if you want.

The entries marked in red are the ones you're viewing now.