Un Blog

La Lengua Española

2011-04-18T08:24:00.000-07:00

"Que buen idioma el mío, qué buena lengua heredamos de los conquistadores torvos..Éstos andaban a zancadas por las tremendas cordilleras, por las Américas encrespadas buscando patatas, butifarras, frijolitos, tabaco negro, oro, maíz, huevos fritos, con aquel apetito voraz que nunca más se ha visto en el mundo...Todo se lo tragaban, con religiones, pirámides, tribus, idolatrías iguales a las que ellos traían en sus grandes bolsas...Por donde pasaban quedaba arrasada la tierra...Pero a los bárbaros se les caían de las botas, de las barbas, de los yelmos, de las herraduras, como piedrecitas las palabras luminosas que se quedaron aquí resplandecientes...el idioma. Salimos perdiendo...Salimos ganando...Se llevaron el oro y nos dejaron el oro...Se lo llevaron todo y nos dejaron todo...Nos dejaron las palabras.”

Pablo Neruda

Cambios

2010-11-26T22:15:00.000-08:00

Veran señores este es mi blog, y yo hago con él lo que me da la gana, no sera mas una obligación con la humanidad sera una obligación personal, por tanto algunos post tendran un toque mas personal. Si a alguno le molesta bien puede irse a ver la pagina de google y buscar un tema interesante (www.google.com).

Gracias.

La Admon.

Introducción al psicoanalisis.

2010-09-27T21:59:00.000-07:00

Ultimamente me he estado leyendo "Introducción al psicoanalisis", intentare plasmar aqui un poco lo que he aprendido y tambien las dudas que me ha traido esta lectura.

Primero que todo el Psicoanalisis según Freud es un tratramiento para tratar neurosis, es decir enfermedades psiquicas. Se basa en el lenguaje entre el paciente y el psicoanalista. Es método cientifico, y como tal debe pararse en los pilares de este (aunque en ocasiones me parece que cojea).

El primer tema a tratar son los actos fallidos. Los actos fallidos son una serie de equivocaciones que son normalmente pasables, como olvidar un nombre, cambiar una palabra por otra, mezclarlas, perder objetos. Son separadas ppalmente en: Equivocaciones orales, equivocaciones escritas, de aprehension, olvidos de compromisos o nombres, olvido de objetos.

La mas tratada es la equivocacion oral. En esta el psicoanalisis propone una explicacion en la cual intervienen dos tendencias, una perturbada y otra perturbadora. La perturbada es la tendencia natural de construccion del discurso, en la cual el sujeto pone su energia, la tendencia perturbadora es aquella que modifica el discurso sin ser advertida por el sujeto. Una cita buena al filosofo Wundt es que en estas ocasiones el poder del discurso es vencido por el poder de la asociacion. Para esta explicacion Freud trae una serie de ejemplos, uno que me causo mucho interes es el de un personaje que tiene por encargo abrir una sesión, pero en el acto dice "Se levanta la sesión", cuando deberia decir: "Se abre la sesión". En el cerebro del sujeto ocurrio algo que le hace modificar el discurso, y cambiar la palabra abrir, por levantar. Inicialmente la explicacion de Freud es que "claramente" la unica interpretacion posible es que el personaje no queria abrir la sesión, y que "cualquier otra explicacion resulta ser superflua" (condenar al fracaso un argumento en contra no parece ser muy acorde al metodo cientifico). Otro caso, es el de un amigo de Freud que tenia un perro enfermo, en este caso el amigo le comenta que "drurara asi un mes mas", a lo cual Freud lo corrige y el amigo responde diciendo que pensaba decir "tristemente durara así un mes", pero combino las palabras.

En esencia todas estas equivocaciones tienen una tendencia original y otra que tiende a perturbarla, la cual suele ser una tendencia escondida y toca analisar al sujeto para encontrar y explicar el porque de la equivocacion.

Django nonrel + MongoDB.

2010-09-27T21:30:00.000-07:00

I want to setup a basic virtualenv environment with a working version of django-nonrel with MongoDB.

From http://www.allbuttonspressed.com/projects/django-nonrel we can read a blog entry that explains everything about the integration: http://www.allbuttonspressed.com/blog/django/2010/05/MongoDB-backend-for-Django-nonrel-released

With requirements.txt:


pymongo
-e git+http://github.com/aparo/django-mongodb-engine.git#egg=django_mongodb_engine
-e hg+http://bitbucket.org/wkornewald/django-nonrel#egg=Django
-e hg+http://bitbucket.org/wkornewald/djangotoolbox#egg=djangotoolbox

Type:


aptitude install python-pip
pip install virtualenv
virtualenv --no-site-packages --distribute nonrel
cd norel
. bin/activate
pip -r ../requirements.txt
bin/django-admin.py startproject myproject

Define your database as:


DATABASES = {
    'default': {
        'ENGINE': 'django_mongodb_engine.mongodb',
        'NAME': 'test',
        'USER': 'test',
        'PASSWORD': 'test',
        'HOST': 'flame.mongohq.com',
        'PORT': 99999,
        'SUPPORTS_TRANSACTIONS': False,
    }
}

Pretty easy.

Identity, Authentication, Authorization and Multitenancy

2010-07-02T13:17:00.000-07:00

There it should somewhere in the web an application to correctly build a Django site with those concepts.

Un poco de genealogia.

2010-06-05T21:30:00.000-07:00

Ayer en un bar un personaje dijo, hablando de Juan Manuel Santos: "No solo nos roban la plata sino tambien la opinion", me surgio la duda de conocer bien las relaciones familiares de los Santos.

Eduardo Santos Montejo (Presidente de Colombia 1938-1942) hermano de Enrique Santos Montejo, este ultimo padre de Enrique Santos Castillo y Hernando Santos Castillos. Tambien es hijo Enrique Santos Molano, con una tal Blanca Molano, de la cual poco encuentro.

Ahora, Enrique Santos Castillo y Hernando Santos Castillo se casaron con dos hermanas de apellido Calderon, y fueron traidas al mundo cuatro criaturas: Juan Manuel Santos Calderon(ex-ministro de Defensa) y Enrique Santos Calderon (director de El Tiempo hasta el 2009) son hijos de Enrique Santos Calderon. Francisco Santos Calderon (Vicepresidente) y Rafael Santos Castillo (director de El Tiempo hasta el 2009) hijos de Hernando Santos Castillo.

Encontre un documento en: http://www.interconexioncolombia.com/documentos/genealogia/que%20esta%20pasando%20en%20colombia/Los%20Santos%20y%20los%20Lopez%20nos%20gobiernan%20desde%201930.pdf
En donde se pueden ver diagramas de Genealogias, es interesante ver en un solo diagrama nombres como: Alvaro Uribe Velez, Jose Obdulio Gaviria, Pablo Escobar Gaviria, Fabio Ochoa, Noemi Sanin.

Basicamente el documento relaciona a todos los actores de este gran teatro llamado Colombia.

Como dicen por ahi: Todo queda en Familia.

Mockus!!!

2010-04-11T23:06:00.000-07:00

Antanas Mockus Presidente de Colombia 2010 - 2014

2010-03-25T20:57:00.000-07:00

En este pequeño espacio posiblemente con pocas visitas intentare ayudar a una justa causa por mi país.

Partido Verde.

Latex + Inline spell check

2010-01-31T11:40:00.000-08:00

Well, after 3 days of trying to activate a good inline spell check on Kile, which is impossible right now, and after trying to compile the hello world plugin for Kate to create a plugin that highlight that misspelled words, which is almost impossible too because of some mistakes on the example, and trying to use emacs for latex + flyspell (it works, but is not very nice), I found that Gedit with its plugin for latex is now really great and has inline spell check!!!, also it list all the structure of the document in a side panel, and has more features, and the best is that is a Gnome app.

http://projects.gnome.org/gedit/
http://live.gnome.org/Gedit/LaTeXPlugin

I love my Linux Kernel.

2010-01-04T20:42:00.000-08:00

Well, I just love it.

Right now I'm working in a driver for a couple of circuits with a DAC, ADC and a encoder to connect my old Gumstix Connex to a PLC Unitronics, and two Variable speed driver Altivar ATV31.

The driver basically read through SPI the data that the ADC get from a current loop. And through SPI the resulting data is send to the drivers. The encoder work with three GPIO interrupt.

I'll try to write later on this. For now, a couple of pages that are really helpful:

http://www.linuxjournal.com/article/5833
http://tomoyo.sourceforge.jp/cgi-bin/lxr/source/
http://www.linuxgrill.com/anonymous/fire/netfilter/kernel-hacking-HOWTO-5.html
http://people.nl.linux.org/ftp/pub/anoncvs/kernelnewbies/documents/kdoc/kernel-api/adt.html#AEN475
http://docs.huihoo.com/linux/kernel/2.6.26/kernel-api/ch25.html
http://docwiki.gumstix.org/index.php/Robostix_spi.c
http://www.makelinux.net/ldd3/ (The chapters with locking techniques, procfs, and deferred work are great!!!)
http://tldp.org/LDP/lkmpg/2.6/html/x1256.html

Corrupted Git

2009-12-28T15:07:00.001-08:00

Well, I have a directory with all my undergraduate project and I use git to ease my work, but recently I was not able to create a backup in a external disk because a error:


remote: Counting objects: 313, done.
remote: Compressing objects: 100% (116/116), done.
fatal: pack has bad object at offset 11893474: inflate returned 1
fatal: index-pack failed

I have files of almost 500MB on this repo, I suppose that the problem. I try to repack but it was useless, I also can check the hash of the object (with fsck), but I can't find where is that file, or what commit it belongs to.

This code will help in this case:


git checkout master
git branch new_master
git checkout new_master
for i in $!$(seq 1 100); do git checkout new_master~$!$i || break; done
for j in $!$(seq $!$i 100); do git checkout new_master~$!$j && break; done
# The error is between new_master~$!$j and new_master~($!$i-2)
git checkout master
git branch -D new_master
git branch new_master
git checkout new_master
git rebase --onto new_master~$!$j new_master~$!$((i-2)) new_master

It works (al least for my problem!!!), and now i want to rebase master:


git checkout master
git rebase new_master
git branch -d new_master

That's all.

Compute the normal vector of a hyperplane in $R^n$

2009-12-03T21:19:00.001-08:00

Well, basically any vector $N$ that holds $N \dot (p_i - p_j) = 0$ for for any vector $p_i, p_j \in H$ $H$ is the hyperplane.

Only $n$ different points are needed, call $P$ the set of $n$ points in $H$, and define $P^* = \{ p^*_i = p_i - p_{i+1} : p_i,p_{i+1} \in H , i \in [1,n-1]\}$.

$P^*$ is the set of vectors on $H$.

Any normal vector $N$ that holds that $N \dot p^*_i , i \in [1,n-1]$ is a normal vector.

$N$ has $n$ components, so we got here $n$ unknowns. And with the $n_1$ $p^*_i$ vectors we got $n-1$ equations.

An extra equation could be defined if we want a normalized vector, but this add a quadratic term in the systems, as the vector with all zero components can't be a normal vector of any hyperplane at least one of its components has a non-zero value, suppose that the $j$ component of $N$ has value $1$. and we have then a Ax = B, with n-1 equations and unknowns. we just have to try which j makes the systm non singular and solve it.

The system of equations:
$\begin{pmatrix}
p^*_{1,1} & p^*_{1,2} & \cdots & p^*_{1,j-1} & p^*_{1,j+1} & \cdots & p^*_{1,n} \\ p^*_{2,1} & p^*_{2,2} & \cdots & p^*_{2,j-1} & p^*_{2,j+1} & \cdots & p^*_{2,n} \\ \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \vdots \\
p^*_{n,1} & p^*_{n,2} & \cdots & p^*_{n,j-1} & p^*_{n,j+1} & \cdots & p^*_{n,n} \\
\end{pmatrix} \begin{pmatrix} N_1 \\ N_2 \\ \vdots \\ N_{j-1} \\ N_{j+1} \\ \vdots \\ N_n \end{pmatrix} = - \begin{pmatrix} p^*_{1,j} \\ p^*_{2,j} \\ \vdots \\ p^*_{(j-1),j} \\ p^*_{(j-1),j} \\ \vdots \\ p^*_{n,j} \end{pmatrix}
$

Start checking for $j \in [1,n]$ and you will get at much $n$ and at least 1 normal vector.

Python:


from numpy import array, dot
from numpy.linalg import solve, norm


def normal(P):
    n = len(P)

    if n < 1:
        raise Exception, "O-size list."

    m = len(P[0])

    if array([len(p) != m for p in P]).any():
        raise Exception, "Inconsistent set of points."

    if n != m:
        raise Exception

    P_ = [p-q for p,q in zip(P[:-1],P[1:])]

    N = []

    for i in range(n):
        try:
            x = solve(
                array([list(p)[:i] + list(p)[i+1:] for p in P_]),
                array([-p[i] for p in P_])
                )

            x = array(list(x)[:i] + [1] + list(x)[i:])
            N.append(x)

        except:
            pass

    n = N[array([(dot(n,n)-1)**2 for n in N]).argmin()]

    return n / norm(n)



if __name__ == '__main__':
    print normal([array([1,0,0]), array([0,1,0]), array([1,1,0])])
    print normal([array([1,123,0123]), array([0345,1,034]), array([132,134,034])])
    print normal([array([1,0,0]), array([0,1,0]), array([1,123,1123])])

2009-12-03T05:51:00.000-08:00

I'm gonna create a image with Debian and Django to upload it to AWS.


dd if=/dev/zero of=debian-ami count=1000 bs=1m
sudo mkfs.ext3 -F debian-ami
mkdir /tmp/chroot
sudo mount -o loop debian-ami /tmp/chroot
sudo debootstrap --arch i386 lenny /tmp/chroot/ http://ftp.debian.org
sudo chroot /tmp/chroot/
# Inside the chroot
mount -t proc none /proc
cd /dev
MAKEDEV console
MAKEDEV std
echo -e 'auto lo\niface lo inet loopback\nauto eth0\niface eth0 inet dhcp' >> /etc/network/interfaces
echo -e 'proc /proc proc defaults 0 0\n/dev/sda1 / reiserfs defaults 0 1\n/dev/sda2 swap swap defaults 0 0' > /etc/fstab
aptitude update
aptitude install locales-all
aptitude install ssh
exit

# In your machine
sudo umount  -l /tmp/chroot 
export EC2_PRIVATE_KEY=xxxxxxx.pem
export EC2_CERT=xxxxxxxxxxxxxxxxxx.pem
export EC2_ACCNO=xxxxxxxxxx
export ACCESS_KEY=xxxxxxxxxxxxxxxx
export SECRET_KEY=xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
ec2-bundle-image -i debian-ami --cert EC2_CERT --privatekey EC2_PRIVATE_KEY -u EC2_ACCNO
ec2-upload-bundle -b 20091203-linux-debian-lenny -m /tmp/debian-ami.manifest.xml -a ACCESS_KEY -s SECRET_KEY
ec2-register -n 20091203-linux-debian-lenny 20091203-linux-debian-lenny/debian-ami.manifest.xml

Algebra

2009-08-23T13:20:00.000-07:00

I'm studying algebra with some friends, and i will try to write here what subjects are we studying right now.

We start with Abstract Algebra using Hungerford book, but i want to follow the MITOpenCourseWare about Algebra I (http://ocw.mit.edu/OcwWeb/Mathematics/18-701Fall-2007/Readings/index.htm).

Composition Law:
A composition law $f$ select a pair of elements on a set $S$ and return an element of $S$, ie:
$f : S \times S \rightarrow S$, I will write $f(a,b)$ as $a*b$ for seek of simplicity.

If $*$ holds that $(a*b)*c = a*(b*c) \quad \forall a,b,c \in S$ then we say that is an associative law.

If $*$ holds that $a*b=b*a \quad \forall a,b \in S$ we say that is a commutative law.

Semigroups

Let $S$ be a set and $*$ be a composition law on $S$, if $*$ is associative then $(S,*)$ is a semigroup.

Monoids

Let $(S,*)$ be a semigroup if $S$ has an identity element $e$ such that:

$\forall x \in S$ $x*e = x = e*x$,

then $(S,*)$ is a Monoid.

Groups

Let $(S,*)$ be a Monoid if:

$\forall x \in S$ $\exists x^{-1} \in S$ s.t. $x*x^{-1} = e = x^{-1}*x$,

then $(S,*)$ is a Group.

if $*$ is commutative then $(S,*)$ is a commutative or Abelian group.

Subgroups

Let $(G,*)$ be a group, if a proper subset $H$ of $G$ is also a group with $*$ then $H$ is a subgroup of $G$

Theorem

If $(G,*)$ is a group and $H$ is a nonempty subset of $G$:

$(H,*)$ is a subgroup of $(G,*)$ if and only if $a*b^{-1} \in H \quad \forall a,b \in H$

Relation

Let $A$ and $B$ be sets and $R$ a map from $A$ to $B$ s.t.:

$R : A \rightarrow B$
with $a\in A$ and $b=R(a) \in B$ and
$(a,b) = (a,R(a)) \in A \times B$
for seek of simplicity we write $(a,R(a)) \in A \times B$ as $aRb$.

Equivalence Relations

Let $S$ be a set and $R$ a relation in $S \times S$, if:

i) $aRa \quad \forall a \in S$ (Reflexive)

ii) $aRb \implies bRa$ (Symmetric)

iii) $aRb \wedge bRc \implies aRc$ (Transitive)

then $R$ is an Equivalence Relation on $S$. (When we speak of an equivalence relation we'll use the symbol $=$)

Partial Order Relations

A partial order relation is a relation $R$ in a set $S$ with a equivalence relation $=$ if $R$ is reflexive and transitive and:

$aRb \wedge bRa \implies a=b$ (antisymmetric)

We write a partial order relation as $\leq$. The above oreder is a non-strict relation.

A strict partial relation order $R$ in $S$ is:

i) $a\not R a \quad \forall a \in S$ (irreflexive)

ii) $aRb \wedge bRa \implies a=b$ (antisymmetric)

iii) $aRb \wedge bRc \implies aRc$ (Transitive)

We write a strict partial order relation as $<$.

Complete Order Relations

The above partial relations don't hold $\forall a,b \in S$, if they do then the order is a complete order.

Congruence Relation

Let $(G,*)$ be a Monoid, if a equivalence relation $=$ on $G$ also holds that:

$a_1=a_2 \wedge b_1 = b_2 \implies a_1*b_1 = a_2*b_2$,

then $=$ is also congruence relation.

Equivalence Classes

Let $S$ be a set with an equivalence relations $=$, the equivalence class $\bar{a}$ (or $[a]$) of an element $a \in S$ is:

$\bar{a} = \{ x \in S | x = a \}$

The class of all the equivalence classes on $S$ is denoted as $(S/=)$, and is called the quotient class of $S$ by $=$.

The union of all the equivalence classes is $S$, i.e.

$\bigcup_{a \in S} \bar{a} = A = \bigcup_{\bar{a} \in (S/=)} \bar{a}$ or $\bar{a} = \bar{b}$

If $a,b \in S$ either $\bar{a} \cap \bar{b} = \emptyset$

Theorem:
Let $(G,*)$ be a Monoid and $=$ a congruence relation on $G$, then:

$((G/=), *)$ is a Monoid with

$* : (G/=) \times (G/=) \rightarrow (G/=)$ and $[a]*[b] = [a*b] \quad [a],[b] \in G/=$ (sorry for the change on notation but the bar don't expand enough in the a*b)

If $(G,*)$ is Abelian then $(G/=,*)$ is Abelian too.

Group Homomorphism

Let $(G,\circ)$ and $(H,\diamond)$ be two groups and a function $f : G \rightarrow H$ such that:

$f(x \circ y ) = f(x) \diamond f(y) \quad \forall x,y \in G$,

then $f$ is a Group homomorphism.

If $f$ is injective then $f$ is a Group Monomorphism.

If $f$ is surjective then $f$ is a Group Epimorphism.

If $f$ is bijective then $f$ is a Group Isomorphism.

If a Group Homomorphism is $f: G \rightarrow G$ then $f$ is a Group Endomorphism of $G$.

If a Group Isomorphism is $f: G \rightarrow G$ then $f$ is a Group Automorphism of $G$.

Kernel

Let $(G,\circ)$ and $(H,\diamond)$ be two groups and a homomorphism $f : G \rightarrow H$. $e_G$ and $e_H$ are the ideintity element of each group, then, the kernel of $f$ is:

$Ker f = \{a \in G | f(a) = e_H\}$,

The kernel is a subset of $G$ that maps to the identity element on $H$ with the homomorphism $f$.

Theorem:
If $f : G \rightarrow H$ is a monomorphism $\Longleftrightarrow$ $Ker f = \{e_G\}$

Theorem:
If $f : G \rightarrow H$ is a isomorphism $\Longleftrightarrow$ $f^{-1} : H \rightarrow G$ is a homomorphism and $f(f^{-1}(a)) = a \quad \forall a \in H$ and $f^{-1}(f(a)) = a \quad \forall a \in G$

Cyclic groups

Let $(G,*)$ be a group and $X$ a subset of $G$.

Let $\{ H_i | i \in I \}$ be the family of all the subgroups ($H_i$) of $G$ that contains $X$ i.e. $X \subset H_i \forall i \in I$

$\bigcap_{i \in I} H_i$ is a subgroup of $G$ generated by $X$ and denoted $\langle X\rangle$

The elements of $X$ are the generators of $\langle X\rangle$.

Several different subsets on $G$ can generate the same $\langle X\rangle$, so, in general $\langle X\rangle = \langle Y \rangle$ with $X \neq Y$.

If $X$ is a finite set such that $X = \{a_1, a_2, \cdots , a_n \}$, we write $\langle a_1, a_2, \cdots, a_n \rangle$ instead of $\langle X\rangle$.

If $G = \langle a_1,a_2, \cdots, a_n \rangle$ with $a_i \in G$, then $G$ is said to be finitely generated.

If $a \in G$, then the subgroup $\langle a\rangle$ is the Cyclic Group or Cyclic Subgroup generaterd by a.

Theorem

Let $(G,*)$ be a group and $X$ a nonempty subset of $G$, then the subgroup $\langle X \rangle$ generated by $X$ consists of all finite products $a_1^{n_1}a_2^{n_2}\cdots a_1t^{n_t} \quad a_i \in X; \quad n_i \in \mathbb{Z}$.

And $\forall a \in G, \langle a \rangle = \{ a^n | n \in \mathbb{Z}\}$

Theorem

Let $H$ be a cycli subgroup of $(\mathbb{Z},+)$, either $H=\langle 0 \rangle$ or $H=\langle m \rangle$, with $m$ the least positive integer in $H$, then if $H \neq \langle 0 \rangle$, then $H$ is infinite.

Theorem

Every infinite cyclic group is isomorphic to $(\mathbb{Z},+)$ and every finite cyclic group of order $m$ is isomorphic to $(\mathbb{Z}_m,+)$ (¿¿¿$\langle m \rangle$???).

Cosets

First we define that $a \equiv b (mod m) \Longleftrightarrow m | a - b \Longleftrightarrow a-b \in \langle m \rangle$.

Let $(G,*)$ be a group, and $(H,*)$ a subgroup of $(G,*)$ and $a,b \in G$.

a is right congruent to $b$ modulo $H$, denoted as $a \equiv_r b (mod H)$ if $a*b^{-1} \in H$.

a is left congruent to $b$ modulo $H$, denoted as $a \equiv_l b (mod H)$ if $a^{-1}*b \in H$.

If $G$ is abelian then if one of the congruencies hold then the other does:

$a*b^{-1} \in H \Longleftrightarrow (a*b^{-1})^{-1} \in H$

$(a*b^{-1})^{-1} = b*a^{-1} = a^{-1}*b$

and

$a^{-1}*b \in H \Longleftrightarrow (a^{-1}*b)^{-1} \in H$

$(a^{-1}*b)^{-1} = b^{-1}*a = a*b^{-1}$

Theorem

Let $(H,*)$ be a subgroup of $(G,*)$, then

i) Right(Left) congruence modulo $H$ is an equivalence relation on $G$.

ii) The equivalence class of $a \in G$ with the right (left) congruence is the set $Ha = \{h*a | h \in H\} ($aH = \{a*h | h \in H\})

iii) $|Ha| = |H| = |aH| \quad \forall a \in G$

$Ha$ is called a right coset of $H$ in $G$.
$aH$ is called a left coset of $H$ in $G$.

Corollary

i) $G$ is the union of the right(left) cosets of $H$ on $G$.

ii) Two right(left) cosets of $H$ on $G$ are either disjoint or equal.

iii) $\forall a,b \in G, \quad Ha = Hb \Longleftrightarrow a*b^{-1} \in H$ and $aH = bH \Longleftrightarrow a^{-1}*b \in H$

iv) If $R$ is the set of distinct right(left) cosets of $H$ in $G$, and $L$ is the set of distinct left cosets of $H$ in $G$, then

Definition

Let $(G,*)$ be a group, and $(H,*)$ a subgroup of $(G,*)$ the index of $H$ in $G$, denoted as $[G : H]$ is the cardinal number of the set of distinct right(left) cosets. of $H$ in $G$.

Runge Kutta in Haskell

2009-08-22T18:56:00.001-07:00

Well, Runge-Kutta is a better solver that Euler so:


-- Runge-Kutta method (t_o, t_f,h, x_o, function(t,x))
runge t_o t_f h x_o f = case t_o < t_f of
    True ->
        do
            let k_1 = f t_o x_o
            let k_2 = f (t_o + h/2) (x_o + h * k_1 / 2)
            let k_3 = f (t_o + h/2) (x_o + h * k_2 / 2)
            let k_4 = f (t_o + h) (x_o + h * k_3)
            [(t_o,x_o)] ++ ( runge (t_o + h) t_f h (x_o + h * (k_1 + 2*k_2 + 2*k_3 + k_4) / 6) f)
    False ->
        []


-- Function
f :: Double -> Double -> Double
f t x = - x

-- Print in column
show_col :: [(Double, Double)] -> String
show_col x = case x of
    (a:b) ->
        do
            let (t, v) = a
            (show t) ++ " , " ++ (show v) ++ "\n" ++ (show_col b)
    [] ->
        ""

main = putStrLn(show_col(runge 0 5 0.01 10 f))

Euler Method in Haskell

2009-08-21T16:02:00.001-07:00

I decide today to learn haskell, but i can't get out of my head how could it be used in simulation, so, i will try with several basic solvers, starting with Euler Method to solve: y'(t) = f(t, y(t))


--Euler method (t_o, t_f,h, x_o, function(t,x))
euler1 :: Double -> Double -> Double -> Double -> (Double -> Double -> Double) -> [(Double, Double)]
euler1 t_o t_f h x_o f = case t_o < t_f of
    True ->
        [(t_o,x_o)] ++ (euler1 (t_o+h) t_f h (x_o+(h * f t_o x_o)) f)
    False ->
        []



--Euler method (t_o, t_f,steps, x_o, function(t,x))
euler2 :: Double -> Double -> Integer -> Double -> (Double -> Double -> Double) -> [(Double, Double)]
euler2 t_o t_f steps x_o f = case steps > 0 of
    True ->
        do
            let h = (t_f - t_o) / (fromIntegral steps)
            [(t_o, x_o)] ++ euler2 (t_o + h) t_f (steps - 1) (x_o + (h * f t_o x_o)) f
    False ->
        []



-- Function
f :: Double -> Double -> Double
f t x = - x

-- Print in column
show_col :: [(Double, Double)] -> String
show_col x = case x of
    (a:b) ->
        do
            let (t, v) = a
            (show t) ++ " , " ++ (show v) ++ "\n" ++ (show_col b)
    [] ->
        ""

--main = putStrLn(show_col(euler1 0 1 0.00001 10 f))
main = putStrLn(show_col(euler2 0 1 100000 10 f))

Qucs

2009-08-13T16:16:00.000-07:00

Mmmm, Creo que por fin encontre un buen programa para simular circuitos, trabaja con SPICE y ademas es libre.

http://qucs.sourceforge.net/
http://qucs.sourceforge.net/docs.html
http://qucs.sourceforge.net/download.html

Symbolic Integration (Manuel Bronstein)

2009-08-12T21:57:00.001-07:00

Well, now that i have a new book, Symbolic Integration I of Manuel Bronstein, i will implement some of the code examples on the initial chapter: polydivide, euclidean algorithm, square free factorization, all of them using GiNaC library.

PolyDivide

This one find the unique q and r in a field K[z] that holds a = b*q + r


lst PolyDivide(ex A, ex B, symbol x)
{
    ex Q = 0;
    ex R = A;

    int d = R.degree(x) - B.degree(x);

    while ((not R.is_zero()) and (d >=0))
    {
        ex T = R.lcoeff(x) / B.lcoeff(x) * pow(x,d);
        Q = expand(Q + T);
        R = expand(R - B*T);
        d = R.degree(x) - B.degree(x);
    }
    return lst(Q,R);
}

PolyPseudoDivide
Same that last, but in a integral domain, like Z.


lst PolyPseudoDivide(ex A, ex B, symbol x)
{
    ex b = B.lcoeff(x);
    int N = A.degree(x) - B.degree(x) + 1;
    ex Q = 0;
    ex R = A;

    int d = R.degree(x) - B.degree(x);

    while ((not R.is_zero()) and (d >=0))
    {
        ex T = R.lcoeff(x) * pow(x,d);
        N -= 1;
        Q = expand(b * Q + T);
        R = expand(b * R - T * B);
        d = R.degree(x) - B.degree(x);
    }

    b = pow(b, N);
    return lst( b * Q, b * R);
}

Well, testing both algorithm against the quo and rem that do the same computation in GiNaC, we have this:

The test is a division between two polynomials of the same n degree randomly generated, in x-axis is n and in y-axis is ellapsed time of the compute code in sec.


int main()
{
    symbol x("x");
    int i = 0;
    ex a=0, b=0;

    clock_t s, e;
    double t1,t2, t3, t4, t5;

    for(int n = 1; n < 100000 ; n += 1000)
    {
        srand ( time(NULL) );

        for(; i < n; i++)
        {
            a += (rand() % 200 - 100) * pow(x,i);
            b += (rand() % 200 - 100) * pow(x,i);
        }

        s = clock();
        PolyDivide(a,b,x);
        e = clock();
        t1 = (e-s);

        s = clock();
        PolyPseudoDivide(a,b,x);
        e = clock();
        t2 = (e-s);

        s = clock();
        quo(a,b,x); rem(a,b,x);
        e = clock();
        t3 = (e-s);

        ex q;
        s = clock();
        q = quo(a,b,x);a - b*q;
        e = clock();
        t4 = (e-s);

        ex r;
        s = clock();
        r = rem(a,b,x);(a-r)/b;
        e = clock();
        t5 = (e-s);


        t1 /= CLOCKS_PER_SEC;
        t2 /= CLOCKS_PER_SEC;
        t3 /= CLOCKS_PER_SEC;
        t4 /= CLOCKS_PER_SEC;
        t5 /= CLOCKS_PER_SEC;

        cout << n << " " << t1 << " " << t2 << " " << t3 << " " << t4 << " " << t5 << endl;
    }

    return 0;
}

Create a new toolbox for Scilab (II)

2009-07-22T19:11:00.000-07:00

Now i gonna follow the instructions at: http://www.scilab.org/contrib/index_contrib.php?page=toolbox_guide

1) First i need to create all the files needed to the toolbox:


toolbox$ cd scilab/scilab/contrib/
toolbox/scilab/scilab/contrib$ mkdir symbolic
toolbox/scilab/scilab/contrib$ cd symbolic/
toolbox/scilab/scilab/contrib/symbolic$ mkdir macros src sci_gateway help etc unit_tests demos  includes
toolbox/scilab/scilab/contrib/symbolic$ touch readme.txt builder.sce loader.sce license.txt

2) Macros, well, basically i only have one macro at the symbolic module, but i gonna test the macro at the instructions, so:


toolbox/scilab/scilab/contrib/symbolic$ cd macros/
toolbox/scilab/scilab/contrib/symbolic/macros$ cat << EOF > foo1.sci
function [X]=foo1(A)
    // This function returns the positive components of the A diagonal

    // Check the type and the size of A
    if  type(A)<>1 then
        error("type of input argument must be a double");
    end
    if  size(A,1)<>size(A,2) then
        error("input argument must be a square matrix");
    end
    //Extraction of the positive components
    X=[];
    for i=1:size(A,1)
        if A(i,i)>0 then
            X($+1)=A(i,i);
        end
    end
endfunction

EOF
toolbox/scilab/scilab/contrib/symbolic/macros$ cat << EOF > buildmacros.sce
mode(-1)
toolboxname='symbolic'
pathB=get_absolute_file_path('buildmacros.sce')
disp('Building macros  in ' +pathB)
genlib(toolboxname+'lib',pathB,%t)
clear pathB genlib toolboxname

EOF
toolbox/scilab/scilab/contrib/symbolic/macros$ cat << EOF > loadmacros.sce
mode(-1)
pathL=get_absolute_file_path('loadmacros.sce')
disp('Loading macros  in ' +pathL)
load(pathL+'/lib')
clear pathL

EOF

3) We're going to create some primitives now, basically is like the module, so, there is a sci_gateway path and a src path, in the first we define the gateways to our functions, and in the second we write the main code of the primitives.


toolbox/scilab/scilab/contrib/symbolic$ cd src/
toolbox/scilab/scilab/contrib/symbolic/src$ cat << EOF > vectsum.c
void vectsum(int n, double * a, double * b, double * y)
{
    int k;
    for (k = 0; k < n; ++k) 
    y[k] = a[k] + b[k];
}

EOF
jcardona@terminus:~/stuff/personal/symbolic/main/toolbox/scilab/scilab/contrib/symbolic/src$ cd ..
toolbox/scilab/scilab/contrib/symbolic$ cd sci_gateway/
toolbox/scilab/scilab/contrib/symbolic/sci_gateway$ cat << EOF > sci_sumab.c
#include "stack-c.h"
extern int vectsum(int n, double * a, double * b, double * y); 

void sci_sumab(char *fname){  
int l1, m1, n1, l2, m2, n2, l3, n;   

/* 1 - Check the number of inputs/outputs arguments  */  
int minlhs=1, maxlhs=1, minrhs=2, maxrhs=2; 
CheckRhs(minrhs,maxrhs) ; 
CheckLhs(minlhs,maxlhs) ; 

/* 2 - Check inputs arguments type, and get the size
       and the address in the Scilab stack of the inputs
       arguments
*/
GetRhsVar(1, "d", &m1, &n1, &l1);
GetRhsVar(2, "d", &m2, &n2, &l2);

/* 3 - Check that the inputs arguments have the same size */
/*     it's possible to use the chekdims and getscalar
       functions to make these checks
*/ 
n=m2*n2; 
if( n1!=n2 || m1!=m2) 
{
    cerro("inputs arguments must have the same size"); 
    return 0; 
}
if(n1!=0 && m1!=0)    
    if(n1!=1 && m1!=1)   
    {
        cerro("inputs arguments must be vectors");    
        return(0); 
    }
 

/* 4 - Create a new variable corresponding to the output argument */ 
CreateVar(3,"d",&m2,&n2,&l3);   

/* 5 -call vectsum routine: returns in stk(l3) the sum of a and b*/
vectsum(n,stk(l1),stk(l2),stk(l3));  

/* 6 - Specif ouput argument */  
LhsVar(1) = 3;  
return 0;
}

EOF

The builder for the primitives, we need a buildsrc.sce and a buildsci_gateway.sce, in the respective path. The first one creates a shared library that is linked to the second shared library (the one of the gateway).


toolbox/scilab/scilab/contrib/symbolic/src$ cat << EOF > buildsrc.sce
//ilib_for_link('symbolicsrc',['fun1.o','fun2.o','vectsum.o'],[],"c")
ilib_for_link('symbolicsrc',['vectsum.o'],[],"c")

EOF
toolbox/scilab/scilab/contrib/symbolic/sci_gateway$ cat << EOF >buildsci_gateway.sce
// must be run from this directory
ilib_name  = 'libsymbolic'     // interface library name
//files = ['sci_fun.o','sci_sumab.o'];  // objects files
files = ['sci_sumab.o'];  // objects files
libs  = ["../src/libsymbolicsrc"]                 // other libs needed for linking
//table = [ /*'fun', 'sci_fun';*/
table = ['sumab','sci_sumab'];        // table of (scilab_name,interface-name)
// do not modify below
ilib_build(ilib_name,table,files,libs)

EOF

4) The help files are located at the help path with a xml format, and a dtd defined by scilab, we only add a template for the files:


toolbox/scilab/scilab/contrib/symbolic/help$ cat << EOF > sumab.xml 
<?xml version="1.0" encoding="UTF-8"?>
<refentry version="5.0-subset Scilab" xml:id="sumab" xml:lang="en"
          xmlns="http://docbook.org/ns/docbook"
          xmlns:xlink="http://www.w3.org/1999/xlink"
          xmlns:svg="http://www.w3.org/2000/svg"
          xmlns:ns3="http://www.w3.org/1999/xhtml"
          xmlns:mml="http://www.w3.org/1998/Math/MathML"
          xmlns:db="http://docbook.org/ns/docbook">
  <info>
    <pubdate>$LastChangedDate: 2008-03-26 09:50:39 +0100 (mer., 26 mars 2008)$</pubdate>
  </info>

  <refnamediv>
    <refname>sumab</refname>

    <refpurpose>Purpose</refpurpose>
  </refnamediv>

  <refsynopsisdiv>
    <title>Calling Sequence</title>

    <synopsis>sequence</synopsis>
  </refsynopsisdiv>

  <refsection>
    <title>Description</title>

    <para>Do something</para>

    <para>Add here a paragraph of the function description </para>
  </refsection>

  <refsection>
    <title>Examples</title>

    <programlisting role="example">exmaple</programlisting>
  </refsection>

  <refsection>
    <title>Authors</title>

    <simplelist type="vert">
      <member>YOUR NAME</member>
    </simplelist>
  </refsection>
</refentry>

EOF
toolbox/scilab/scilab/contrib/symbolic/help$ cat sumab.xml | sed -e 's?sumab?foo1?' > foo1.xml
toolbox/scilab/scilab/contrib/symbolic/help$ cat << EOF > buildhelp.sce
mode(-1) //force silent execution
path=get_absolute_file_path('buildhelp.sce');//get the absolute path of this file
add_help_chapter("Symbolic",path);//add help chapter
xmltohtml(path,"Symbolic")
//clear the variable stack
clear path add_help_chapter get_absolute_file_path 

EOF
toolbox/scilab/scilab/contrib/symbolic/help$ cat << EOF > loadhelp.sce
mode(-1) //force silent execution
path=get_absolute_file_path('loadhelp.sce');//get the absolute path of this file
add_help_chapter("Symbolic",path);//add help chapter
clear path add_help_chapter get_absolute_file_

EOF

5) Finally we set the global builder and loader:


toolbox/scilab/scilab/contrib/symbolic$ cat << EOF > builder.sce
mode(-1);
mainpathB=get_absolute_file_path('builder.sce');
chdir(mainpathB);
if isdir('src') then
chdir('src');
exec('buildsrc.sce');
chdir('..');
end
if isdir('sci_gateway') then
chdir('sci_gateway');
exec('buildsci_gateway.sce');
chdir('..');
end
if isdir('macros') then
chdir('macros');
exec('buildmacros.sce');
chdir('..');
end
if isdir('help') then
chdir('help');
exec('buildhelp.sce');
chdir('..');
end
clear mainpathB

EOF
toolbox/scilab/scilab/contrib/symbolic$ cat << EOF > loader.sce
mode(-1);
mainpathL=get_absolute_file_path('loader.sce');
chdir(mainpathL);
if isdir('sci_gateway') then
chdir('sci_gateway');
exec('loader.sce');
chdir('..');
end
if isdir('macros') then
chdir('macros');
exec('loadmacros.sce');
chdir('..');
end
if isdir('help') then
chdir('help');
exec('loadhelp.sce');
chdir('..');
end
clear mainpathL

EOF

6) Build it at scilab:


toolbox/scilab/scilab$ bin/scilab

-->exec("contrib/symbolic/builder.sce")
 
-->mode(-1);
   Generate a loader file
   Generate a Makefile
   ilib_gen_Make: Copy compilation files (Makefile*, libtool...) to TMPDIR
   ilib_gen_Make: Copy vectsum.c to TMPDIR
   ilib_gen_Make: Modification of the Makefile in TMPDIR.
   Running the Makefile
   Generate a cleaner file
 ans  =
 
 libsymbolicsrc.so   
 ilib_name  =
 
 libsymbolic   
 libs  =
 
 ../src/libsymbolicsrc   
   Generate a gateway file
   Generate a loader file
   Generate a Makefile
   ilib_gen_Make: Copy compilation files (Makefile*, libtool...) to TMPDIR
   ilib_gen_Make: Copy sci_sumab.c to TMPDIR
   ilib_gen_Make: Copy libsymbolic.c to TMPDIR
   ilib_gen_Make: Modification of the Makefile in TMPDIR.
   Running the makefile
   Generate a cleaner file
 
 Building macros  in /home/jcardona/stuff/personal/symbolic/main/toolbox/scilab/scilab/contrib/symbolic/macros/   

Building the master document:
 SCI/contrib/symbolic/help

Building the manual file [html] in SCI/contrib/symbolic/help. (Please wait building ... this can take a while)
Warning : redefining function: get_absolute_file_path  . Use funcprot(0) to avoid this message

-->exec("contrib/symbolic/loader.sce")
 
-->mode(-1);
Shared archive loaded.
Link done.
Shared archive loaded.
Link done.
 
 Loading macros  in /home/jcardona/stuff/personal/symbolic/main/toolbox/scilab/scilab/contrib/symbolic/macros/   
 
-->sumab(1,2)
 ans  =
 
    3.  
 
-->sumab(2,2)
 ans  =
 
    4.

We have then a brand new toolbox enjoy summing.

Create a new toolbox for Scilab (I)

2009-07-21T19:23:00.000-07:00

As i have to change my symbolic module to a symbolic toolbox, i want to write this to use as a future reference.

1) Get Scilab:


$ mkdir toolbox
$ cd toolbox/
toolbox$ git clone git://git.scilab.org/scilab

2) First compilation:



toolbox$ cd scilab/scilab

toolbox/scilab/scilab$ ./configure
...
checking hdf5.h usability... no
checking hdf5.h presence... no
checking for hdf5.h... no
configure: error: Cannot find headers (hdf5.h) of the library HDF5. Please install the dev package

3) HDF5 Headers:


sudo aptitude install libhdf5-serial-dev

4) Installing JHDF5 :


wget http://www.hdfgroup.org/ftp/HDF5/hdf-java/src/hdf-java-2.5-src.tar
tar -xvf hdf-java-2.5-src.tar
cd hdf-java/
./configure --with-jdk=/usr/lib/jvm/java-6-sun-1.6.0.14/include/,/usr/lib/jvm/java-6-sun-1.6.0.14/lib/
make
sudo mkdir -p /usr/share/java/jar
sudo cp lib/jhdf5.jar /usr/share/java/jar

5) Configuring again:


toolbox/scilab/scilab$ ./configure
toolbox/scilab/scilab$ make
toolbox/scilab/scilab$ bin/scilab

Well, we now have our scilab base.
We have to set all the toolboxes files, but that will be tonight.

Lock free programming

2009-05-13T21:33:00.000-07:00

http://en.wikipedia.org/wiki/Lock-free_and_wait-free_algorithms
http://en.wikipedia.org/wiki/Compare-and-swap

http://www.google.com/search?hl=en&client=iceweasel-a&rls=org.debian%3Aen-US%3Aunofficial&q=lock+free+&btnG=Search

Jack en gumstix.

2009-05-13T18:49:00.000-07:00

MMM, necesito un framework de simulacion por bloques que este ligado la kernel, se me vino a la cabeza jack y su forma de pasar datos entre plugins, y simplemente vi esto en unos slides:

Although not a requirement, Jack should support any streaming data type, not just audio.

Podria ser una buena opcion meter jack en el gumstix y crear los bloques de control como si fueran plugins y el paso de datos serian rapidito, el llamaod a los codigos de los bloques. Suena interesante la cosa. Veremos que sale.

Drivers en linux (II)

2009-04-28T20:23:00.001-07:00

Ahora crearemos un dispositivo de caracteres en /dev/algo, para esto debemos definir cuatro nuevas funciones: para abrir y cerrar el dispositivo, y para leer y escribirlo. Hay mas operaciones, pero con esas es suficiente por ahora (http://www.linuxtopia.org/online_books/linux_kernel/linux_kernel_module_programming_2.6/x569.html).

El archivo test.h queda asi:


#ifndef _TEST_SPI_H_
#define _TEST_SPI_H_

#define DEVICE_NAME "test-spi"

MODULE_AUTHOR("Jorge Eduardo Cardona <jorgeecardona@gmail.com>");
MODULE_DESCRIPTION("SPI test driver.");
MODULE_LICENSE("GPL");

// Character device functions.
static int this_open(struct inode *inode, struct file *file);
static int this_release(struct inode *inode, struct file *file);
static ssize_t this_read(struct file *file, char __user *buffer, size_t len, loff_t *offset);
static ssize_t this_write(struct file *file, const char __user *buffer, size_t len, loff_t *offset);

// File Operations structure.
struct file_operations this_fops = {
    .open       = this_open,
    .release    = this_release,
    .read       = this_read,
    .write      = this_write
};

// Linux module functions.
static int __init this_init(void);
static void __exit this_exit(void);
module_init(this_init);
module_exit(this_exit);

#endif

Para crear el dispositivo de caracteres debemos llamar a la funcion: register_chrdev (http://kerneltrap.org/man/linux/man9/register_chrdev.9)(http://www.silicontao.com/ProgrammingGuide/GNU_function_list/register_chrdev.html)


#include <linux/module.h>
#include <linux/fs.h>
#include <linux/err.h>


#include "test.h"

// Global data
static int this_major;


// Function that is called every time we try to open a device.
static int this_open(struct inode *inode, struct file *file)
{
    printk(DEVICE_NAME": Open device.\n");
    return 0;
}

// Function that is called when we close a device.
static int this_release(struct inode *inode, struct file *file)
{
    printk(DEVICE_NAME": Close device.\n");
    return 0;
}

// Function that is called when we read from a device.
static ssize_t this_read(struct file *file, char __user *buffer, size_t len, loff_t *offset)
{
    printk(DEVICE_NAME": Read device.\n");
    return 0;
}

// Function that is called when we write to a device.
static ssize_t this_write(struct file *file, const char __user *buffer, size_t len, loff_t *offset)
{
    printk(DEVICE_NAME": Write device.\n");
    return 0;
}

// Function that is called when we  load the module.
static int __init this_init(void)
{
    int ret;

    printk(DEVICE_NAME ": Init.\n");

    // Register char device: 0 for dynamically allocation of the major number.
    ret = register_chrdev(0, DEVICE_NAME, &this_fops);
    if (~IS_ERR_VALUE(ret))
    {
        printk(DEVICE_NAME ": Character device registered, with major=%d.\n",ret);
        this_major = ret;

        return 0;
    }

    unregister_chrdev(this_major, DEVICE_NAME);
    printk(DEVICE_NAME ": Character device unregistered.\n");

    return ret;
}

// Function that is called when we remove the module.
static void __exit this_exit(void)
{
    unregister_chrdev(this_major, DEVICE_NAME);
    printk(DEVICE_NAME ": Character device unregistered.\n");
    printk(DEVICE_NAME ": Exit.\n");
}

(http://tomoyo.sourceforge.jp/cgi-bin/lxr/ident?i=IS_ERR_VALUE)

Bueno, igual que antes, configuramos, compilamos y tenemos:


root@gumstix-custom-connex:~$ modprobe test
root@gumstix-custom-connex:~$ dmesg | grep "test-spi"
<4>test-spi: Init.
<4>test-spi: Character device registered, with major=254.
root@gumstix-custom-connex:~$ mknod /dev/test c 254 0
root@gumstix-custom-connex:~$ cat /dev/test 
root@gumstix-custom-connex:~$ dd bs=1 count=0 if=/dev/zero of=/dev/test
0+0 records in
0+0 records out
root@gumstix-custom-connex:~$ echo "ja" > /dev/test

root@gumstix-custom-connex:~$ rmmod test
root@gumstix-custom-connex:~$ dmesg | grep "test-spi"
<4>test-spi: Init.
<4>test-spi: Character device registered, with major=254.
<4>test-spi: Open device.
<4>test-spi: Read device.
<4>test-spi: Close device.
<4>test-spi: Open device.
<4>test-spi: Close device.
<4>test-spi: Open device.
<4>test-spi: Write device.
<4>test-spi: Write device.
<4>test-spi: Write device.
<4>test-spi: Write device.
<4>test-spi: Write device.
<4>test-spi: Write device.
<4>test-spi: Write device.
<4>test-spi: Write device.
<4>test-spi: Write device.
<4>test-spi: Write device.
<4>test-spi: Close device.
<4>test-spi: Character device unregistered.
<4>test-spi: Exit.
root@gumstix-custom-connex:~$

Los varios mensajes de Write se deben a que retornamos 0 esto quiere decir que no escribimos ningun caracter y la funcion es llamada hasta que vacie el buffer, pero aqui toca matarla con un Ctrl-C.

Bueno, un paso mas.

Drivers en linux (I)

2009-04-28T07:52:00.000-07:00

Bueno iba a hacer esto en un solo post pero por partes es mejor, primero lo básico, cargar y bajar el modulo. En mi caso necesito un modulo para SPI, entonces lo ubico en drivers/spi/test.c

Comenzando con lo basico, entrar y bajar el modulo.

test.c


#include <linux/module.h>
#include "test.h"

static int __init this_init(void)
{
    // Initialize module.

    printk(DEVICE_NAME ": Init.\n");
}

static void __exit this_exit(void)
{
    // Exit module.

    printk(DEVICE_NAME ": Exit.\n");
}

test.h


#ifndef _TEST_SPI_H_
#define _TEST_SPI_H_

#define DEVICE_NAME "test-spi"

MODULE_AUTHOR("Jorge Eduardo Cardona <jorgeecardona@gmail.com>");
MODULE_DESCRIPTION("SPI test driver.");
MODULE_LICENSE("GPL");

// Linux module functions.
static int __init this_init(void);
static void __exit this_exit(void);
module_init(this_init);
module_exit(this_exit);

#endif

Tenemos que manipular dos archivos mas: Kconfig y Makefile, ambos en la misma carpeta:

Kconfig


.
.
.
config SPI_TEST
    tristate "Test SPI driver."
    depends on SPI_MASTER && EXPERIMENTAL
    help
      Test SPI driver.
.
.
.

Makefile


.
.
.
obj-$(CONFIG_SPI_TEST)      += test.o
.
.
.

Al configurar el kernel nos saldra la opcion, la seleccionamos y podremos hacer esto:


root@gumstix-custom-connex:~$ modprobe test
root@gumstix-custom-connex:~$ dmesg | tail -n1
<4>; test-spi: Init.
root@gumstix-custom-connex:~$ rmmod test
root@gumstix-custom-connex:~$ dmesg | tail -n1
<4> test-spi: Exit.
root@gumstix-custom-connex:~$

Primer paso listo.

Contando racionales.

2009-04-15T15:12:00.001-07:00

Un numero racional es todo aquel que puede ser escrito de la forma:

racional = a/b : siendo a y b numeros enteros.

Usando el metodo como Georg Ferdinand Ludwig Philipp Cantor demostro que los racionales son contables podemos crear un pequeñito programa para sacar el i-esimo racional de la tabla de cantor:


from math import sqrt

def rational(i):
    b = int((sqrt(1+8*i)-1)/2)
    n = i - (b**2 + b)/2 + 1
    d = b + 2 - n

    return (n,d)


for i in xrange(100):
    (n,d) = rational(i)

    print "(%d) %d/%d \t= %f"%(i,n,d,float(n)/float(d))

Bueno, La tabla de Cantor es sencilla, el numerador aumenta columna a columna desde 1 hasta aleph_0 y el denumenador fila a fila de 1 a aleph_0, de esta forma todos los nuemros racionales son cubiertos por la tabla, y por tanto son contables, se les puede asociar uno a uno un elemento de los enteros.

La tabla:


1/1  2/1  3/1  4/1  5/1  6/1  7/1  8/1  9/1  ...
1/2  2/2  3/2  4/2  5/2  6/2  7/2  8/2  9/2  ...
1/3  2/3  3/3  4/3  5/3  6/3  7/3  8/3  9/3  ...
1/4  2/4  3/4  4/4  5/4  6/4  7/4  8/4  9/4  ...
1/5  2/5  3/5  4/5  5/5  6/5  7/5  8/5  9/5  ...
1/6  2/6  3/6  4/6  5/6  6/6  7/6  8/6  9/6  ...
1/7  2/7  3/7  4/7  5/7  6/7  7/7  8/7  9/7  ...
1/8  2/8  3/8  4/8  5/8  6/8  7/8  8/8  9/8  ...
1/9  2/9  3/9  4/9  5/9  6/9  7/9  8/9  9/9  ...

Y el conteo:


0    2    5    9    14   20   27   35   44   ...
1    4    8    13   19   26   34   43   53   ...
3    7    12   18   25   33   42   52   ...  
6    11   17   24   32   41   51   ...
10   16   23   31   40   50   ...
15   22   30   39   49   ... 
21   29   38   48   ... 
28   37   47   ...
36   46   ...
45   ...
...

Resultado de la ejecucion del codigo:


(0) 1/1  = 1.000000
(1) 1/2  = 0.500000
(2) 2/1  = 2.000000
(3) 1/3  = 0.333333
(4) 2/2  = 1.000000
(5) 3/1  = 3.000000
(6) 1/4  = 0.250000
(7) 2/3  = 0.666667
(8) 3/2  = 1.500000
(9) 4/1  = 4.000000
(10) 1/5  = 0.200000
(11) 2/4  = 0.500000
(12) 3/3  = 1.000000
(13) 4/2  = 2.000000
(14) 5/1  = 5.000000
(15) 1/6  = 0.166667
(16) 2/5  = 0.400000
(17) 3/4  = 0.750000
(18) 4/3  = 1.333333
(19) 5/2  = 2.500000
(20) 6/1  = 6.000000
(21) 1/7  = 0.142857
(22) 2/6  = 0.333333
(23) 3/5  = 0.600000
(24) 4/4  = 1.000000
(25) 5/3  = 1.666667
(26) 6/2  = 3.000000
(27) 7/1  = 7.000000
(28) 1/8  = 0.125000
(29) 2/7  = 0.285714
(30) 3/6  = 0.500000
(31) 4/5  = 0.800000
(32) 5/4  = 1.250000
(33) 6/3  = 2.000000
(34) 7/2  = 3.500000
(35) 8/1  = 8.000000
(36) 1/9  = 0.111111
(37) 2/8  = 0.250000
(38) 3/7  = 0.428571
(39) 4/6  = 0.666667
(40) 5/5  = 1.000000
(41) 6/4  = 1.500000
(42) 7/3  = 2.333333
(43) 8/2  = 4.000000
(44) 9/1  = 9.000000
(45) 1/10  = 0.100000
(46) 2/9  = 0.222222
(47) 3/8  = 0.375000
(48) 4/7  = 0.571429
(49) 5/6  = 0.833333
(50) 6/5  = 1.200000
(51) 7/4  = 1.750000
(52) 8/3  = 2.666667
(53) 9/2  = 4.500000
(54) 10/1  = 10.000000
(55) 1/11  = 0.090909
(56) 2/10  = 0.200000
(57) 3/9  = 0.333333
(58) 4/8  = 0.500000
(59) 5/7  = 0.714286
(60) 6/6  = 1.000000
(61) 7/5  = 1.400000
(62) 8/4  = 2.000000
(63) 9/3  = 3.000000
(64) 10/2  = 5.000000
(65) 11/1  = 11.000000
(66) 1/12  = 0.083333
(67) 2/11  = 0.181818
(68) 3/10  = 0.300000
(69) 4/9  = 0.444444
(70) 5/8  = 0.625000
(71) 6/7  = 0.857143
(72) 7/6  = 1.166667
(73) 8/5  = 1.600000
(74) 9/4  = 2.250000
(75) 10/3  = 3.333333
(76) 11/2  = 5.500000
(77) 12/1  = 12.000000
(78) 1/13  = 0.076923
(79) 2/12  = 0.166667
(80) 3/11  = 0.272727
(81) 4/10  = 0.400000
(82) 5/9  = 0.555556
(83) 6/8  = 0.750000
(84) 7/7  = 1.000000
(85) 8/6  = 1.333333
(86) 9/5  = 1.800000
(87) 10/4  = 2.500000
(88) 11/3  = 3.666667
(89) 12/2  = 6.000000
(90) 13/1  = 13.000000
(91) 1/14  = 0.071429
(92) 2/13  = 0.153846
(93) 3/12  = 0.250000
(94) 4/11  = 0.363636
(95) 5/10  = 0.500000
(96) 6/9  = 0.666667
(97) 7/8  = 0.875000
(98) 8/7  = 1.142857
(99) 9/6  = 1.500000

El sueño de la razón produce monstruos. (Goya)

2008-11-02T08:53:00.000-08:00

GtkModelica

2008-10-07T09:39:00.000-07:00

Mmm, mezclare un poco de C++ con Anjuta, GTK con glade, FlowCanvas, y Modelica con OpenModelica, ojala salga algo bueno, y rapido!!.

Hoy encontre oro!

2008-09-22T17:58:00.000-07:00

http://en.wikipedia.org/wiki/Modelica
http://www.modelica.org/
http://www.ida.liu.se/labs/pelab/modelica/OpenModelica.html

Dynamic Programming

2008-07-14T22:20:00.000-07:00

Bien, hoy lei un poco de Dynamic Programming, es bacancito, de una pense en un decorator en Python, no es complicado, ahi esta el ejemplo:

1) Fibonacci normal:


#!/usr/bin/python
import sys

def fib(x):
    if x==0:
        return 0
    if x==1:
        return 1
    return fib(x-1)+fib(x-2)

print fib(int(sys.argv[1]))

Este codigo tiene un comportamiento como lo muestra la grafica (N vs segundos)

2) Ahora con un decorator implementando Dynamic Programming.


#!/usr/bin/python
import sys

map={}
def dynamic_programming(func):
  def decorator(*args, **kwargs):
      if (func,args) in map:
          return map[(func,args)]
      else:
          r = func(*args, **kwargs)
          map[(func,args)] = r
          return r
  return decorator


@dynamic_programming
def fib(x):
  if x==0:
      return 0
  if x==1:
      return 1
  return fib(x-1)+fib(x-2)


print fib(int(sys.argv[1]))

Como pueden ver el codigo de la funcion fib no se cambio, solo se le agrego un decorator, el cual busca primero si el argumento esta en un diccionario, para no recalcularlo. (kwargs es un dict, y no puedo usarlo como un key en otro dict, toco sacarlo del key :))

No pondre grafica de este solo esto:

$ time ./fib2.py 499
86168291600238450732788312165664788095941068326060883324529903470149056115823592713458328176574447204501

real 0m0.029s
user 0m0.024s
sys 0m0.004s

OS hello world.

2008-06-15T16:59:00.000-07:00

Je, hace rato que no escribía nada, vamos a ver si sale un hello world de un """sistema operativo""", en VirtualBox.

Primero organizamos el dispositivo de arranque, que seria una imagen de un disco. Bueno tenemos que entender bien lo de MBR, eso de cilindros, sectores y cabezas.

Sector: Agrupación lógica de 512 Bytes en el disco. La medida que muestran programas como fdisk es Sector/Track, o sea cantidad de sectores en una pista del disco.

Track: Es una "circunferencia de datos", en cada superficie magnetica. Podemos calcular la cantidad de bytes en un track con sector/track*512.

Head: Representa la cantidad de dispositivos lectores magnéticos, seria como la cantidad de superficies magnéticas disponibles en el disco, aunque la medida que entregan programas como fdisk no representa esto, puede tomarse como una especie de medida para poder ubicar un byte en el disco.

Cylinder: Estas es la cantidad total de tracks que existen en una sola superficie magnética.

Entonces podemos decir que la unidad de almacenamiento es 1 Bit, un byte lo componen 8 Bits, un sector lo componen 512 Bytes, una pista(track) la componen (sector/track) tracks o (sector/track)*512 Bytes, una superficie magnetica esta compuesta por (cylinders) tracks, o (cylinders)*(sector/track)*512 Bytes, y el disco completo esta compuesto por (heads) superficies magneticas o (heads)*(cylinders)*(sectors/track)*512 Bytes. (http://www.storagereview.com/guide2000/ref/hdd/index.html)

MBR: Master Boot Record, la misma esta ubicada en el primer sector del disco (primeros 512 Bytes).

Ahora necesitamos generar una imagen con mbr, si queremos un disco con 16 heads, 63 sectors/track y 160 cylinders:


$ dd if=/dev/zero of=disk.image bs=512 count 161280

Ahora con losetup montamos este archivo como un dispositivo de bloques:


$ sudo losetup /dev/loop0 disk.image

Creamos el MBR con:


$ sudo fdisk -u -C160 -S63 -H16 /dev/loop0

Con -C decimos los cilindros, con -S los sectores por pista, con -H las cabezas, y con -u decimos que nos muestre los tamaños en sectores y no en cilindros. Ahora creamos una particion primaria de todo el disco con el flag de boot activo, la tabla de particiones debe quedar:


Command (m for help): p

Disk /dev/loop0: 82 MB, 82575360 bytes
16 heads, 63 sectors/track, 160 cylinders, total 161280 sectors
Units = sectors of 1 * 512 = 512 bytes
Disk identifier: 0x19c7acc5

Device Boot      Start         End      Blocks   Id  System
/dev/loop0p1   *          63      161279       80608+  83  Linux

En este punto tenemos una imagen de un disco solo con MBR, en http://en.wikipedia.org/wiki/Mbr podemos ver una descripcion de las cabeceras del MBR, y podemos verificar en nuestra imagen con:


$ hd -v -n512 disk.image
00000000  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................|
00000010  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................|
00000020  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................|
00000030  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................|
00000040  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................|
00000050  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................|
00000060  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................|
00000070  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................|
00000080  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................|
00000090  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................|
000000a0  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................|
000000b0  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................|
000000c0  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................|
000000d0  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................|
000000e0  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................|
000000f0  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................|
00000100  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................|
00000110  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................|
00000120  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................|
00000130  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................|
00000140  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................|
00000150  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................|
00000160  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................|
00000170  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................|
00000180  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................|
00000190  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................|
000001a0  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................|
000001b0  00 00 00 00 00 00 00 00  c5 ac c7 19 00 00 00 01  |................|
000001c0  01 00 83 0f 3f 9f 3f 00  00 00 c1 75 02 00 00 00  |....?.?....u....|
000001d0  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................|
000001e0  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................|
000001f0  00 00 00 00 00 00 00 00  00 00 00 00 00 00 55 aa  |..............U.|
00000200

Esto corresponde al primer sector de nuestro disco, podemos ver la entrada para la particion primaria que hicimos, la cual es:


$ hd -v -s446 -n16 disk.image
000001be  00 01 01 00 83 0f 3f 9f  3f 00 00 00 c1 75 02 00  |......?.?....u..|
000001ce

Aun falta darle formato a la particion que creamos. La anterior cabecera nos dice que esta particion comienza en el sector 0x3F000000, (little-endian -> 63), y tiene una extension de 0xC1750200 (161217) sectores.

A esta particion le damos formato asi:


$ sudo losetup -o32256 /dev/loop1 disk.image
$ sudo mkfs.ext2 -b1024 /dev/loop1 80608

Sabemos que podemos comenzar a hacer uso de nuestra particion desde el sector 63, que equivale al byte 63*512=32256, y tiene una extension de 161217*512 bytes, y como manejamos bloques de 1024 en ext2, tendriamos un total de 80608.5 bloques que redondeamos a 80608.

Ahora tenemos que instalar el grub en nuestra imagen, esta es la parte mas confusa, es buena una cuanta documentacion (http://www.gnu.org/software/grub/ , http://www.gnu.org/software/grub/manual/multiboot/ )

Debemos instalar todo lo necesario en nuestro sistema de archivos en la carpeta /boot/grub.


$ mkdir mnt
$ sudo mount -text2 /dev/loop1 mnt/
$ wget http://ftp.br.debian.org/debian/pool/main/g/grub/grub_0.97-27_i386.deb
$ dpkg-deb -X grub_0.97-27_i386.deb grub_debian
$ sudo mkdir -p mnt/boot/grub
$ sudo cp grub_debian/usr/lib/grub/i386-pc/* mnt/boot/grub/
$ sudo grub

grub> device (hd0) /dev/loop0

grub> geometry (hd0) 160 16 63
drive 0x80: C/H/S = 160/16/63, The number of sectors = 161280, /dev/loop0
Partition num: 0,  Filesystem type is ext2fs, partition type 0x83

grub> root (hd0,0)
Filesystem type is ext2fs, partition type 0x83

grub> setup (hd0)
Checking if "/boot/grub/stage1" exists... yes
Checking if "/boot/grub/stage2" exists... yes
Checking if "/boot/grub/e2fs_stage1_5" exists... yes
Running "embed /boot/grub/e2fs_stage1_5 (hd0)"...  15 sectors are embedded.
succeeded
Running "install /boot/grub/stage1 (hd0) (hd0)1+15 p (hd0,0)/boot/grub/stage2 /boot/grub/menu.lst"... failed

Error 22: No such partition

grub> quit

Aqui bajamos el paquete de debian de grub, del cual copiamos todo lo necesario para que grub funcione en una arquitectura i386 que es la que pensamos virtualizar. Encontramos un eror en grub, pero si tiramos el comando con strace vemos lo siguiente(http://osdir.com/ml/boot-loaders.grub.bugs/2005-01/msg00035.html) :


$ sudo strace -o trace.log grub
...
$ cat trace.log | grep loop
open("/dev/loop0", O_RDONLY|O_LARGEFILE) = 4
open("/dev/loop0", O_RDWR|O_LARGEFILE)  = 3
open("/dev/loop01", O_RDWR|O_LARGEFILE) = -1 ENOENT (No such file or directory)

Grub cree que la particion 1 de nuestra imagen se encuentra en loop01, entonces ahi se la daremos:


$ sudo ln -s /dev/loop1 /dev/loop01

Volvemos a tirar grub:


$ sudo grub

grub> device (hd0) /dev/loop0

grub> geometry (hd0) 160 16 63
drive 0x80: C/H/S = 160/16/63, The number of sectors = 161280, /dev/loop0
Partition num: 0,  Filesystem type is ext2fs, partition type 0x83

grub> root (hd0,0)
Filesystem type is ext2fs, partition type 0x83

grub> setup (hd0)
Checking if "/boot/grub/stage1" exists... yes
Checking if "/boot/grub/stage2" exists... yes
Checking if "/boot/grub/e2fs_stage1_5" exists... yes
Running "embed /boot/grub/e2fs_stage1_5 (hd0)"...  15 sectors are embedded. succeeded
Running "install /boot/grub/stage1 (hd0) (hd0)1+15 p (hd0,0)/boot/grub/stage2 /boot/grub/menu.lst"... succeeded
Done.

grub> quit

Ya podemos verificar la instalacion del stage1 en el mbr:


$ hd -v -n512 disk.image

00000000  eb 48 90 00 00 00 00 00  00 00 00 00 00 00 00 00  |.H..............|
00000010  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................|
00000020  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................|
00000030  00 00 00 00 00 00 00 00  00 00 00 00 00 00 03 02  |................|
00000040  ff 00 00 20 01 00 00 00  00 02 fa 90 90 f6 c2 80  |... ............|
00000050  75 02 b2 80 ea 59 7c 00  00 31 c0 8e d8 8e d0 bc  |u....Y|..1......|
00000060  00 20 fb a0 40 7c 3c ff  74 02 88 c2 52 be 7f 7d  |. ..@|<.t...R..}| 
00000070  e8 34 01 f6 c2 80 74 54  b4 41 bb aa 55 cd 13 5a  |.4....tT.A..U..Z| 
00000080  52 72 49 81 fb 55 aa 75  43 a0 41 7c 84 c0 75 05  |RrI..U.uC.A|..u.| 
00000090  83 e1 01 74 37 66 8b 4c  10 be 05 7c c6 44 ff 01  |...t7f.L...|.D..| 
000000a0  66 8b 1e 44 7c c7 04 10  00 c7 44 02 01 00 66 89  |f..D|.....D...f.| 
000000b0  5c 08 c7 44 06 00 70 66  31 c0 89 44 04 66 89 44  |\..D..pf1..D.f.D| 
000000c0  0c b4 42 cd 13 72 05 bb  00 70 eb 7d b4 08 cd 13  |..B..r...p.}....| 
000000d0  73 0a f6 c2 80 0f 84 ea  00 e9 8d 00 be 05 7c c6  |s.............|.| 
000000e0  44 ff 00 66 31 c0 88 f0  40 66 89 44 04 31 d2 88  |D..f1...@f.D.1..| 
000000f0  ca c1 e2 02 88 e8 88 f4  40 89 44 08 31 c0 88 d0  |........@.D.1...| 
00000100  c0 e8 02 66 89 04 66 a1  44 7c 66 31 d2 66 f7 34  |...f..f.D|f1.f.4| 
00000110  88 54 0a 66 31 d2 66 f7  74 04 88 54 0b 89 44 0c  |.T.f1.f.t..T..D.| 
00000120  3b 44 08 7d 3c 8a 54 0d  c0 e2 06 8a 4c 0a fe c1  |;D.}<.T.....L...| 
00000130  08 d1 8a 6c 0c 5a 8a 74  0b bb 00 70 8e c3 31 db  |...l.Z.t...p..1.| 
00000140  b8 01 02 cd 13 72 2a 8c  c3 8e 06 48 7c 60 1e b9  |.....r*....H|`..| 
00000150  00 01 8e db 31 f6 31 ff  fc f3 a5 1f 61 ff 26 42  |....1.1.....a.&B| 
00000160  7c be 85 7d e8 40 00 eb  0e be 8a 7d e8 38 00 eb  ||..}.@.....}.8..| 
00000170  06 be 94 7d e8 30 00 be  99 7d e8 2a 00 eb fe 47  |...}.0...}.*...G| 
00000180  52 55 42 20 00 47 65 6f  6d 00 48 61 72 64 20 44  |RUB .Geom.Hard D| 
00000190  69 73 6b 00 52 65 61 64  00 20 45 72 72 6f 72 00  |isk.Read. Error.| 
000001a0  bb 01 00 b4 0e cd 10 ac  3c 00 75 f4 c3 00 00 00  |........<.u.....| 
000001b0  00 00 00 00 00 00 00 00  cc 69 95 5f 00 00 80 01  |.........i._....| 
000001c0  01 00 83 0f 3f 9f 3f 00  00 00 c1 75 02 00 00 00  |....?.?....u....| 
000001d0  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................| 
000001e0  00 00 00 00 00 00 00 00  00 00 00 00 00 00 00 00  |................| 
000001f0  00 00 00 00 00 00 00 00  00 00 00 00 00 00 55 aa  |..............U.| 
00000200

Comencemos con el codigo, primero necesitamos la porcion de codigo que nos salta al codigo en C que queremos ejecutar, esta primera parte viene en assembly, aqui hay uno que nos puede servir (tomado de http://javiervalcarce.es/wiki/Kernel_multiboot_m%C3%ADnimo):



/* Start-up */

#define STACK_SIZE              0x4000      /* The size of our stack (16KB) */
#define MULTIBOOT_HEADER_MAGIC  0x1BADB002      /* The magic number for the Multiboot header */

        .text
        .globl  _start

    .align  4               /* Align 32 bits boundary */
                            /* Multiboot-compliant header */
        .long   +MULTIBOOT_HEADER_MAGIC         /* magic */
        .long   0                   /* flags=0 */
        .long   -MULTIBOOT_HEADER_MAGIC         /* checksum */
_start:
        movl    $(stack + STACK_SIZE), %esp /* Initialize the stack pointer */

        push   $0                   /* Reset EFLAGS */
        popf

        call    cmain               /* C main function... */
loop:   hlt
        jmp     loop

        .comm   stack, STACK_SIZE       /* Our stack area */

lo guardamos en un archivo llamado startup.S, ahora el codigo en C:


/* kernel.c - the C part of the kernel */

/* Some screen stuff. */
#define COLUMNS         80      /* The number of columns */
#define LINES           24      /* The number of lines */
#define ATTRIBUTE       7       /* The attribute of an character */
#define VIDEO           0xB8000 /* The video memory address */

/* Variables */
static int xpos;        /* Save the X position */
static int ypos;        /* Save the Y position */
static volatile
unsigned char *video;   /* Point to the video memory */

/* Forward declarations */
void cmain (unsigned long magic, unsigned long addr);
static void cls();
static void putchar (int c);


////////////////////////////////////////////////////////////////////////////////
void cmain (unsigned long magic, unsigned long addr)
{

char msg[] = "Hello World";
char*   p;

cls();

/* printf */
p = msg;
for (p = msg; *p != 0; p++)
{
putchar(*p);
}
}

////////////////////////////////////////////////////////////////////////////////
/* Clear the screen and initialize VIDEO, XPOS and YPOS */
static void cls (void)
{
int i;

video = (unsigned char *) VIDEO;

for (i = 0; i < xpos =" 0;" ypos =" 0;" c ="=" c ="=" xpos =" 0;">= LINES)
ypos = 0;
return;
}

*(video + (xpos + ypos * COLUMNS) * 2)     = c & 0xFF;
*(video + (xpos + ypos * COLUMNS) * 2 + 1) = ATTRIBUTE;

xpos++;
if (xpos >= COLUMNS)
goto newline;
}

Guardamos esto en un kernel1.c, y compilamos :


$ gcc -c kernel1.c
$ gcc -c startup.S
$ ld -o kernel1 kernel1.o startup.o -Ttext 0x100000

Ahora terminamos de configurar el grub para que funcione con nuestro kernel1:


$ sudo cp kernel1 mnt/boot/

Creamos el archivo menu.lst con:


default         0
timeout         2
color cyan/blue white/blue

title           kernel1
root            (hd0,0)
kernel          /boot/kernel1

Y lo copiamos a boot/grub/menu.lst:


$ sudo cp menu.lst  mnt/boot/grub/

Ya deberiamos tener nuestro ambiente listo, falta convertir esta imagen en un disco que virtualbox pueda leer:


$ sudo umount mnt
$ sudo losetup -d /dev/loop1
$ sudo losetup -d /dev/loop0
$ sudo rm /dev/loop01
$ vditool DD PruebaOS.vdi disk.image

Prendemos nuestro VirtualBox, organizamos una maquina virtual cuyo disco sea el que acabamos de crear, mmmm, bueno, grub arranco, sale la lista con kernel1, intenta correrlo, pero sale un error "Error 13: Invalid or unsupported executable format" , aaaa, buen detalle:


$ uname -m
x86_64

De seguro si tienen una maquina i386 les funciona de una.
Listo, entonces compilo en chroot de i386:


$ cp kernel1.c  startup.S  /var/etch-i386-20080604/tmp/
$ sudo chroot /var/etch-i386-20080604/ /bin/bash
# cd /tmp
# gcc -c kernel1.c
# gcc -c startup.S
# ld -o kernel1 kernel1.o startup.o -Ttext 0x100000
# exit
$ sudo losetup -o32256 /dev/loop1 disk.image
$ sudo mount -text2 /dev/loop1 mnt/
$ sudo cp /var/etch-i386-20080604/tmp/kernel1 mnt/boot/kernel1
$ sudo umount mnt
$ sudo losetup -d /dev/loop1
$ rm PruebaOS.vdi
$ vditool DD PruebaOS.vdi disk.image

Otra vez abrimos virtualbox, borramos la configuracion con el disco anterior, ponemos el nuevo y corremos la maquina.

Listo eso es todo, de ahí a hacer un kernelcito decente hay un paso, jaja.

Ayuda a la comunidad

2008-06-09T17:49:00.000-07:00

Luego de 6 meses sin entradas: La pagina web de la empresa de unos amigos SIA Ingenieros, Soluciones en Ingeniería Aplicada

Botka, The Barely Standing Robot

2007-12-10T19:03:00.001-08:00

Crazy Multi-Input Touch Screen

2007-10-02T21:40:00.000-07:00

Transformando a Medellín

2007-09-10T20:35:00.000-07:00

Libertad es una palabra que el sueño humano alimenta, que no hay nadie que la explique y nadie que no la entienda.

2007-09-04T18:20:00.000-07:00

Sin embargo los ojos incerrables del che miran como si no pudieran no mirar.

2007-09-03T12:41:00.000-07:00

Che
Mario Benedetti

Lo han cubierto de afiches / de pancartas
de voces en los muros
de agravios retroactivos
de honores a destiempo

lo han transformado en pieza de consumo
en memoria trivial
en ayer sin retorno
en rabia enmbalsamada

han decidido usarlo como epilogo
como ultima thule de la inocencia vana
como anejo arquetipo de santo o satanas

y quizas han resuelto que la unica forma
de desprenderse de El
o dejarlo al garete
es vaciarlo de lumbre
convertirlo en un heroe
de marmol o de yeso
y por lo tanto inmovil
o mejor como mito
o silueta o fantasma
del pasado pisado

sin embargo los ojos incerrables del che
miran como si no pudieran no mirar
asombrados tal vez de que el mundo no entienda
que treinta anos despues siga bregando
dulce y tenaz por la dicha del hombre.

Por que software libre?

2007-08-06T18:49:00.000-07:00

If I have seen further [than others] it is by standing on the
shoulders of giants... Issac Newton

Tiempo para Tesis.

2007-05-02T21:06:00.000-07:00

Bueno, tal parece que le voy a poder sacar tiempo a la tesis el próximo semestre.
TODO:
- Comenzar anteproyecto.
- Conseguir director.
- Terminar anteproyecto.
- Revisar anteproyecto.
- Entregar anteproyecto.

Tesis TODO:
- Estudiar UML.
- Estudiar B-Method.
- Estudiar los problemas que se presentan en control.
- Model de control en UML.
- Estudiar la possibilidad y el como se haria el control distribuido.
- Modelo del framework de control en UML.
- Paso de UML a B-method.
- Verificacion.
- Implementacion de Framework de control (Sobre Pistachio y Gumstix)
- Buscar un sistema cuyo control necesite tiempo real, y controlarlo con el Framework sobre el Gumstix)

Tengo como maximo 12 meses para la tesis, toca repartir bien eso.

Tesis.

2007-04-16T19:16:00.000-07:00

Bueno, en el almuerzo se me ocurrio la idea que podria ser mi tesis (ojala algun dia decida bien que hacer). Siguiendo con la idea del anterior post, seria asi:

En gumstix se correria pistachio, en el cual se correria un Framework de Control, el cual estaria previamente diseñado , modelado con UML y verificado con B-Method, y para la configuracion del control se tendria que especificar en UML.
Se tendria algo asi:

* Una aplicacion en la cual se diseñarian los componentes del sistema: Actuadores, Sensores, y el control como tal. La salida de este seria un modelo en UML.

* Este modelo se le pasaria al Framework en Pistachio el cual lo implementaria y se comunicaria con los sensores y los actuadores.

* El Framework se diseña en UML y se verificaria en B-Method, si algo falla ya seria problema de la arquitectura del Gumstix.

* Se pondria a prueba en el control de algun sistema complejo de controlar.

Eso seria todo.

Tesis

2007-04-15T19:05:00.000-07:00

Bueno no hay GSoC para mi, entonces vuelvo con mi tesis, pensandolo bien, creo que voy a redirigir mi tesis, tengo puntos en lso cuales quiero trabajar:

* Control digital en Tiempo-Real
* Redes neuronales.
* Pistachio
* Alguna implementacion en Gumstix.

Hasta ahora no se bien que sera. Pero si se unen las partes seria correr pistachio en Gumstix, y encima correr un Framework de Tiempo Real para el control de "Algo" y ojala que con Redes Neuronales.

GSoC

2007-04-05T11:56:00.000-07:00

Some important lectures about Plan 9 and Inferno.

http://cm.bell-labs.com/sys/doc/index.html
http://www.vitanuova.com/inferno/papers/bltj.html
http://www.vitanuova.com/inferno/papers/descent.html
http://www.vitanuova.com/inferno/papers/limbo.html
http://plan9.bell-labs.com/wiki/plan9/recommended_readings/index.html

GSoC

2007-03-27T18:43:00.000-07:00

Bueno, le tire al GSoC, con lo que fue el principio de la idea de construir un sistema distribuido a partir de sistemas embebidos junto con Torres, hace como 2 años.
Al principio la idea era Inferno ( por el styx) , pero creo que limbo asusta. Seria bueno que se diera lo del GSoC, depronto saldria una implementacion de styx para pistachio. (Cumpliendo las obligaciones con el mes de marzo)

De argentina a México en bus

2007-02-02T18:22:00.000-08:00

http://www.soho.com.co/wf_InfoArticulo.aspx?idArt=5108

Corriendo wombat

2007-01-28T21:15:00.001-08:00

Bueno, despues de bajarme y utilizar un poco Fiasco, hice lo mismo con el paquete de Pistachio que estan trabajando en NICTA y UNSW, junto con Iguana y Wombat.
Lo que hice fue bajarme del svn todos los paquetes necesarios. En esta pagina muestran como bajarse los paquetes (http://www.ertos.nicta.com.au/software/baz-howto.pml) y para poder compilarlo toca hacerlo con un cross-compiler, no sirve el gcc normalcito, me baje este de nicta (http://www.ertos.nicta.com.au/downloads/i686-gcc-3.3.4-glibc-2.3.3-2006-06-02.tar.gz)
y posicionado en la carpeta donde estan todos los archivos ejecuto los siguiente:
tools/scons.py machine=pc99 toolprefix="i686-unknown-linux-gnu-" simulate
y listo, empiezan a compilarse todas las partes necesarias, pistachio, iguana, y wombat, y queda una imagen build/c.img que se corre automatica en qemu.
Una muestra:

More Lectures

2007-01-21T18:25:00.000-08:00

Operating Systems
Advanced Operating Systems
Distributed Systems

Lectures

2007-01-09T19:20:00.000-08:00

Microkernel Construction

Building Microkernel-Based Operating Systems

De como compilar Fiasco.

2006-12-30T21:01:00.000-08:00

Comence a compilar Fiasco, lo que hice fue seguir las instrucciones en http://os.inf.tu-dresden.de/fiasco/use.html.
1) Me baje los modulos grub l4 dde hello l4env dope preprocess de el cvs.
2) En el directorio l4 ejecuto lo siguiente:
make O=../build_user config
(lo dejo como esta)
make O=../build_user all
(Espero como 20 minutos y sale un error de sed, para arreglarlo toca entrar a l4/tool/gcc-wrap/Makefile y cambiar a:

$(VERBOSE)sed '\
s,@system@,$(SYSTEM),g;\
s,@ar@,$(AR),g;\
s,@nm@,$(NM),g;\
s,@ranlib@,$(RANLIB),g;\
s,@cc@,$(CC_$(DEFAULT_ARCH)),g;\
s,@cxx@,$(CXX_$(DEFAULT_ARCH)),g;\
s,@l4dir@,$(absfilename $(L4DIR)),g;\
s,@dropsdir@,$(DROPS_STDDIR),g;\
s,@compilerincdir@,$(GCCDIR_$(DEFAULT_ARCH))/include,g;\
s,@libgcc@,$(GCCLIB_$(DEFAULT_ARCH)),g;\
s,@libgcc_eh@,$(GCCLIB_EH_$(DEFAULT_ARCH)),g;\
s,@default_system@,$(DEFAULT_SYSTEM),g;\
s,@default_arch@,$(DEFAULT_ARCH),g;\
s,@default_cpu@,$(call PROCESS_SYSTEM, 's/[^-_]*_$[^-]*$.*/\1/p'),g;\
s,@default_l4api@,$(call PROCESS_SYSTEM, 's/.*-$.*$/\1/p'),g;\
' $< > $@.new

por:

$(VERBOSE)sed 's,@system@,$(SYSTEM),g;s,@ar@,$(AR),g;s,@nm@,$(NM),g;s,@ranlib@,$(RANLIB),g;s,@cc@,$(CC_$(DEFAULT_ARCH)),g;s,@cxx@,$(CXX_$(DEFAULT_ARCH)),g;s,@l4dir@,$(absfilename $(L4DIR)),g;s,@dropsdir@,$(DROPS_STDDIR),g;s,@compilerincdir@,$(GCCDIR_$(DEFAULT_ARCH))/include,g;s,@libgcc@,$(GCCLIB_$(DEFAULT_ARCH)),g;s,@libgcc_eh@,$(GCCLIB_EH_$(DEFAULT_ARCH)),g;s,@default_system@,$(DEFAULT_SYSTEM),g;s,@default_arch@,$(DEFAULT_ARCH),g;s,@default_cpu@,$(call PROCESS_SYSTEM, 's/[^-_]*_$[^-]*$.*/\1/p'),g;s,@default_l4api@,$(call PROCESS_SYSTEM, 's/.*-$.*$/\1/p'),g;' $< > $@.new

No se por que mi sed no puede con las instrucciones separadas en lineas con un \,pero se arregla el problema. Luego lo dejo otra vez compilando, y me da tiempo de salir, comprar dos libros (Calicalabozo de Andrés Caicedo, y La Despedida de Milan Kundera) , comer algo y cuando vuelvo la maquina esta en un error porque no tenia instalado el doxygen, lo instalo y vuelvo a compilar, me da tiempo de organizar los libros que tengo por leer (Calicalabozo, La despedida, Despues del Hombre de Gonzalo Arango, Gato o soledad en la lluvia de Amílcar Osorio, y El libro Negro de Orhan Pamuk), tambien juego un poco en el Wii de mi hermanito, el control es loquisimo, funciona basicamente con acelerometros. por lo que no es necesario jugar parado y moverse como loco, sino tambien sentado, la parte de infrarrojos funciona siempre y cuando no le caiga sol al sensor (pense que eso solo nos pasaba en los carritos que nos ponian a hacer en la universidad) sino se enloquece, pero casi siempre esta en lugar cerrado entonces no es problema, y me imagino que funciona con conceptos de vision esteroscopica; termino de jugar y sigue compilando, me como un postre y espero, luego sale un error porque no encuentra el comando fig2dev y me toca instalar el transfig, y comienzo de nuevo, al ratico (doy gracias que no toca compilar de nuevo lo que ya ha compilado) sale un error: no puede copiar el archivo building.tex en una carpeta, ¿¿bug del Makefile de Fiasco?? el Makefile hace ejecutar el comando, con respuesta:

cp building.tex drops-building.tex
cp: cannot stat `building.tex': No such file or directory

Lo que pasa es que el Makefile entra en el directorio de trabajo, build-user/doc/.... e intenta copiar el archivo desde esa misma carpeta y hay no esta, esta en el directorio l4/doc/building-howto, solo es copiarlo de esta posicion a build_user/doc/building-howto/ y comenzar de nuevo a compilar, y ya si termina, ¡mas largo que la primera instalada de Xfree86 fue esto!.

Ahora entro en la carpeta l4/kernel/fiasco, corro el comando:

make BUILDDIR=../../../build-kernel
cd ../../../build-kernel
make menuconfig
(lo dejo como esta)
make
(Espero mientras veo Carnivale, esta buena la serie)

Sigo los pasos de la pagina para iniciarlo, un poco distinto (pa correrlo en vmware desde una imagen de un floppy) pero en esencia lo mismo.
Y listo, tengo al buen Fiasco corriendo en una maquina virtual.

Distribuyendo mi carga

2006-12-18T16:05:00.000-08:00

Rafa, para responder a tu comentario decidi hacer 2 cosas, comenzar a bajarme el Fiasco del cvs (y no seguir jugando con esa version de Fiasco que tenia desde hace mucho), y copiar y pegar el comentario en un nedit, separarlo por espacios y abordarlo linea por linea.

1. Middleware DICE, lenguaje de IDL y generador de stub y skeleton para los servicios de userspace.

Fiasco: http://os.inf.tu-dresden.de/DICE/
Pistachio: http://www.l4ka.org/projects/idl4/
NICTA: http://ertos.nicta.com.au/software/kenge/magpie/latest/

En Fiasco van en la version 3.1.0 (Nov/2006) y en Pistachio en la 1.0.0 (2003), y el de Nicta en la 1.0 (Dic/2004) , no he trabajado con ninguna de las tres, entonces no tengo nada que decir.

2. Codigo en C++ en el kernel. Bien organizado y con un preprocesador modificado que hace el codigo aun mas limpio.

En pistachio tambien esta en C++, ahora, lo del preprocesador, no creo que exista en pistachio, y no se hasta que punto sea una gran ventaja, confio plenamente en cpp.

3. Device Driver Environment: Entorno que permite reutilizar todos los drivers de Linux en L4

Fiasco: http://os.inf.tu-dresden.de/l4env/ (no encontre pagina de DDE)
Pistachio: http://www.l4ka.org/projects/virtualization/drivers.php

He leido un poco la forma como lo implementaron en Pistachio y me convence, aunque deberia leerlo mas a fondo.

4. El Build System es demasiado poderoso y organizado. Hace checkout del CVS y veras lo que te digo

Nunca he construido Fiasco del todo, solo he usado una imagen que habia en la pagina, y fue hace un buen rato.
Pistachio lo he compilado varias veces,y he construido aplicaciones de prueba, y si es un poco complicado, pero si uno sigue las ayudas no se pierde. (de hecho ayer recompile desde cero lo del CVS y no tuve problema).

En cuanto a bajarlo del CVS, en esas ando como desde hace 2 horas y nada que baja, ojala algun dia baje.

Uno de los aspectos que mas me hacen pensar en Pistachio:
http://ertos.nicta.com.au/press/051124_L4_Qualcomm_vfinal.pdf

Plataformas soportadas:
Fiasco:

  * IA32 (Intel, AMD and compatible CPUs, i486 and above)
* ARM  (SA1110 and PXA)
* UX   (Linux usermode emulation of Fiasco/L4

Pistachio:

Alpha (21164, 21264)
AMD64 (Opteron 242, Simics)
ARM (SA1100, XScale, ARM925T)
IA32 (Pentium and higher)
IA64 (Itanium1, Itanium2, Ski)
MIPS 64bit (R4000, R5000)
PowerPC 32bit (IBM 750)
PowerPC 64bit (Power3, Power4)
((((Pronto M68K jeje))))

Al fin de cuentas, tal parece que ambos tiran para el mismo lado, pero me interesa mas pistachio, por muchas cosas, entre unas el hecho que lo sigan trabajando en UNSW y en NICTA en Australia.

Aun asi, voy a comenzar a trabajar tambien en Fiasco, para ver bien cual se acomoda mejor a lo que estoy buscando.

Pero si me tenes un puestico en Dresden me paso de una a Fiasco, y volvemos a Fiasco "El uKernel".

Volviendo a pistachio

2006-12-17T20:53:00.000-08:00

Me estoy bajando el codigo fuente de pistachio, hace un buen rato que no lo veia, y ya ni podia entender mi enredo de idrectorios, asi, que va a tocar casi desde cero. La idea seria:
1) Lograr controlar una tarjeta de red con pistachio. Primero voy a intentar con el codigo para el stack tcp/ip que viene con OSKIT, aunque igual creo que lo primero es manipular la tarjeta a bajo nivel a mano, y luego conectar esto con el oskit.
2) Hacer una aplicacion cliente/servidor entre pistachios
3) Entender bien el concepto de hilo en pistachio, la generacion de los mismos, y cuales son y para que sirven los hilos basicos de pistachio.

(Esta semana tengo que llamar a averiguar acerca del precio de las tarjetas de las CPLDs)

2006-12-04T21:05:00.000-08:00

http://groups.google.com.co/group/comp.os.linux.advocacy/browse_thread/thread/f4bf728dc301176a/a10dcb9293cd0938?lnk=st&q=&rnum=3&hl=es#a10dcb9293cd0938

Buscando un foco

2006-11-25T13:43:00.000-08:00

Bueno, no se porque busco un foco, habiendo luz ambiente, pero bueno, parece ser necesario para poder terminar haciendo algo, la cosa es la siguiente, el fin ultimo de un computador es resolver problemas, pero el meollo esta en que hay ciertos problemas que no se pueden resolver, y hay otros problemas que aunque se pueden resolver, es inutil esperar una respuesta por parte de un computador actual, entonces la cuestion radica en cambiar por completo el como se construiria el modelo matematico del sistema que quiero desarrollar, no puedo basarme en un modelo que siga a turing, pues este se basa en un modelo secuencial de calculo en el cual para resolver un problema se siguen unos pasos en donde cada paso depende de el anterior, y no es posible realizar todos los pasos en un mismo momento, de esto salen otros problemas, los problemas indistribuibles, no se bien si estos problemas se puede resolver con una maquina cuantica de calculo, pero es claro que con una maquina determinista y con un orden secuencial de calculo, no se pueden resolver ciertos problemas, la idea seria poder ya sea escoger que problemas quisiera poder resolver, o lograr demostrar cuales problemas puede resolver un sistema como el que tengo pensado.

The Kidnap Of Angry Daz

2006-10-30T17:33:00.000-08:00

The Kidnap Of Angry Daz

Que es un computador?

2006-10-19T00:05:00.000-07:00

He estado leyendo acerca de teoria de conjuntos para poder lograr a entender bien el B-Method y en general los metodos de demostracion, tambien algo de teoria de automatas. Me gustario poder comenzar a trabajar con la definicion formal de una maquina de Turing, y a partir de ella llegar al porque es necesario todas las tecnicas que se utilizan actualmente en la implementacion de los computadores.

Basicamente una maquina de turing es un automata que tiene contacto con un medio de almacenamiento. (http://www.warthman.com/ex-turing.htm)

Una buena lectura seria el "On Computable Numbers, with an Application to the Entscheidungsproblem".

Una correccion a la tarjeta de la CPLD

2006-10-15T23:23:00.000-07:00

V1.1

B-Method

2006-10-14T17:16:00.000-07:00

Existe una tesis donde hacen un trabajo con el B-Method para demostrar que pistachio esta bien diseñado ( "A Formal Model of the L4 u-kernel API Using the B Method") 2.

Acerca del B-Method.

Para instalar el BToolKit:
http://www.csse.monash.edu.au/courseware/cse4213/2006/resources/linuxInstall.html

Un resumen.
http://www.csse.monash.edu.au/courseware/cse4213/2006/resources/B-summary.pdf

Interoperabilidad de componentes.
http://www.iist.unu.edu/newrh/III/1/docs/techreports/report333/P05.pdf

Bueno, buscare mas documentacion que me sirva, la idea seria poder usar esa herramienta en el paso de volver pistachio distribuido.

Board para la CPLD y el UltraISR

2006-10-14T14:49:00.000-07:00

Bueno, por fin termine la version 0.0.01 Beta de la tarjeta para hacer pruebas con la cpld 39K200, y tambien el UltraISR para poder programarlas, toca ver cuanto me cobran por las tarjetas.

Bueno, falta revisar bien esto, y ver si le puedo hacer mejoras.

una flor que no seca el invierno

2006-10-04T22:43:00.000-07:00

una flor que no seca el invierno

UltraISR

2006-10-01T21:49:00.000-07:00

Casi que no lo encuentro, pero al fin encontre el documento para hacer un UltraISR, para poder programar las CPLD's de Cypress:
http://www2.mta.ac.il/~flaxer/edu/course/vhdl/isrdesign.pdf

Ahora puedo hacer pruebas con las cpld's, la cosa es conseguir la plata para hacer las tarjetas para las cpld's.

TODO

2006-09-27T19:18:00.000-07:00

*Sistema Embebido con el 68K.
-hacer programador de las CPLD's de Cypress.
-diseñar MMU.
-diseñar modulo de audio.
-diseñar modulo de video.
-diseñar modulo RTC.
-integrar serial, video, audio, MMU, RTC, memoria, LEDS, DIPSWITCH, LCD,
-estudiar posibilidad de un medio masivo de almacenamiento.
-Buscar un SO facil de portar, ojala Pistachio.

*Sistema embebido con el 5407
*Sistema embebido con el 5485
*Estudiar posibilidad de portar Pistachio a los Coldfire, o saltar a un ARM.
*En base a lo anterior hacer un Pistachio distribuido.
*Hacer pruebas de rendimiento.

(Demas que falta algo, si algo lo agrego luego)

La madre de los motorola

2006-09-16T15:58:00.000-07:00

Hace un tiempo me tope con un buen integrado de la familia Coldfire de Freescale, el 5407, tengo 3 de estos integrados en mi caja de electonica, y no quice dejarlos llenar de polvo, me puse a la tarea de crear un sistema embebido con el integrado.
Como no sabia como funcionaba bien, me puse a investigar y cai en la familia 68 de la cual se basa la familia coldfire, y comence a diseñar un sistema completo basado en el 68000.
El sistema consta de:
1 Procesador MC68HC000FN16
2 Memorias SRAM de 128K c/u (D431000). 256K de Ram total.
2 Memorias EEPROM de 8K c/u (AT28C64B) . 64K de memoria para el "Bios" , Bootloader, como le quieran llamar.
1 Display de Cristal Liquido. Pa ver algo
8 Leds a modo de Post Code. Pa depurar.
2 Dip Switch. Uno a modo de entrada, y otro para seleccionar el numero de interrupcion.
2 Pulsadores, pa reset y pa generar una interrupcion.
1 74LS138. pa decodificar direcciones. los 16MB que direcciona el 68000 esta partido en 8.
Glue Logic:
Un monton de compuertas.
1 2n2222 todo bamba (por no poner un IC de negadores).

Mucho cable.
Condensadores (el tipico 104)
Resistencias.

Actualmente esta problematico, creo que hay algunos cables medio sueltos y es problematico pa que funcione, pero lo hace.
Fotos:

Listo, pa que corriera algo, la cosa es sencilla, primero tengo que tener el toolchain del 68K bien compilado (demas que luego digo como, es canson el proceso), y luego hacer un archivo .S donde hay un .init, y ahi un jump a el main de algun programa en C, y listo, hay van los archivos necesarios para que corra lo que tiene actualmente.

De primero la licencia de todo el codigo a continuacion:

Program test to 68000 processor.
Copyright (C) 2005  Jorge Eduardo Cardona Gaviria

This program is free software; you can redistribute it and/or
modify it under the terms of the GNU General Public License
as published by the Free Software Foundation; either version 2
of the License, or (at your option) any later version.

This program is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
GNU General Public License for more details.

You should have received a copy of the GNU General Public License
along with this program; if not, write to the Free Software
Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA.

---------------------------------- boot.S -----------------------------------

.extern start_main
.extern int_002,int_003
.extern int_004,int_005,int_006,int_007,int_008,int_009,int_010,int_011,int_012,int_013
.extern int_014,int_015,int_016,int_017,int_018,int_019,int_020,int_021,int_022,int_023
.extern int_024,int_025,int_026,int_027,int_028,int_029,int_030,int_031,int_032,int_033
.extern int_034,int_035,int_036,int_037,int_038,int_039,int_040,int_041,int_042,int_043
.extern int_044,int_045,int_046,int_047,int_048,int_049,int_050,int_051,int_052,int_053
.extern int_054,int_055,int_056,int_057,int_058,int_059,int_060,int_061,int_062,int_063
.extern int_064,int_065,int_066,int_067,int_068,int_069,int_070,int_071,int_072,int_073
.extern int_074,int_075,int_076,int_077,int_078,int_079,int_080,int_081,int_082,int_083
.extern int_084,int_085,int_086,int_087,int_088,int_089,int_090,int_091,int_092,int_093
.extern int_094,int_095,int_096,int_097,int_098,int_099,int_100,int_101,int_102,int_103
.extern int_104,int_105,int_106,int_107,int_108,int_109,int_110,int_111,int_112,int_113
.extern int_114,int_115,int_116,int_117,int_118,int_119,int_120,int_121,int_122,int_123
.extern int_124,int_125,int_126,int_127,int_128,int_129,int_130,int_131,int_132,int_133
.extern int_134,int_135,int_136,int_137,int_138,int_139,int_140,int_141,int_142,int_143
.extern int_144,int_145,int_146,int_147,int_148,int_149,int_150,int_151,int_152,int_153
.extern int_154,int_155,int_156,int_157,int_158,int_159,int_160,int_161,int_162,int_163
.extern int_164,int_165,int_166,int_167,int_168,int_169,int_170,int_171,int_172,int_173
.extern int_174,int_175,int_176,int_177,int_178,int_179,int_180,int_181,int_182,int_183
.extern int_184,int_185,int_186,int_187,int_188,int_189,int_190,int_191,int_192,int_193
.extern int_194,int_195,int_196,int_197,int_198,int_199,int_200,int_201,int_202,int_203
.extern int_204,int_205,int_206,int_207,int_208,int_209,int_210,int_211,int_212,int_213
.extern int_214,int_215,int_216,int_217,int_218,int_219,int_220,int_221,int_222,int_223
.extern int_224,int_225,int_226,int_227,int_228,int_229,int_230,int_231,int_232,int_233
.extern int_234,int_235,int_236,int_237,int_238,int_239,int_240,int_241,int_242,int_243
.extern int_244,int_245,int_246,int_247,int_248,int_249,int_250,int_251,int_252,int_253
.extern int_254,int_255
.global _start
.global _vectors
.text

_vectors:

.long 0x320000,0x400
.long int_002,int_003
.long int_004,int_005,int_006,int_007,int_008,int_009,int_010,int_011,int_012,int_013
.long int_014,int_015,int_016,int_017,int_018,int_019,int_020,int_021,int_022,int_023
.long int_024,int_025,int_026,int_027,int_028,int_029,int_030,int_031,int_032,int_033
.long int_034,int_035,int_036,int_037,int_038,int_039,int_040,int_041,int_042,int_043
.long int_044,int_045,int_046,int_047,int_048,int_049,int_050,int_051,int_052,int_053
.long int_054,int_055,int_056,int_057,int_058,int_059,int_060,int_061,int_062,int_063
.long int_064,int_065,int_066,int_067,int_068,int_069,int_070,int_071,int_072,int_073
.long int_074,int_075,int_076,int_077,int_078,int_079,int_080,int_081,int_082,int_083
.long int_084,int_085,int_086,int_087,int_088,int_089,int_090,int_091,int_092,int_093
.long int_094,int_095,int_096,int_097,int_098,int_099,int_100,int_101,int_102,int_103
.long int_104,int_105,int_106,int_107,int_108,int_109,int_110,int_111,int_112,int_113
.long int_114,int_115,int_116,int_117,int_118,int_119,int_120,int_121,int_122,int_123
.long int_124,int_125,int_126,int_127,int_128,int_129,int_130,int_131,int_132,int_133
.long int_134,int_135,int_136,int_137,int_138,int_139,int_140,int_141,int_142,int_143
.long int_144,int_145,int_146,int_147,int_148,int_149,int_150,int_151,int_152,int_153
.long int_154,int_155,int_156,int_157,int_158,int_159,int_160,int_161,int_162,int_163
.long int_164,int_165,int_166,int_167,int_168,int_169,int_170,int_171,int_172,int_173
.long int_174,int_175,int_176,int_177,int_178,int_179,int_180,int_181,int_182,int_183
.long int_184,int_185,int_186,int_187,int_188,int_189,int_190,int_191,int_192,int_193
.long int_194,int_195,int_196,int_197,int_198,int_199,int_200,int_201,int_202,int_203
.long int_204,int_205,int_206,int_207,int_208,int_209,int_210,int_211,int_212,int_213
.long int_214,int_215,int_216,int_217,int_218,int_219,int_220,int_221,int_222,int_223
.long int_224,int_225,int_226,int_227,int_228,int_229,int_230,int_231,int_232,int_233
.long int_234,int_235,int_236,int_237,int_238,int_239,int_240,int_241,int_242,int_243
.long int_244,int_245,int_246,int_247,int_248,int_249,int_250,int_251,int_252,int_253
.long int_254,int_255

_start:

move.l #0x00100000,%a0 /* Leds */

check_ram:
move.w #0x80,(%a0)
move.l #0x00300000,%a1

bucle1:
move.l %a1,(%a1)+
cmpa.l #0x00320000,%a1
blt bucle1

bucle2:
move.l -(%a1),%d0
cmp.l %a1,%d0
bne Error
cmpa.l #0x00300000,%a1
bgt bucle2

move.w #0xff,(%a0)

move.l #0x100000,%d1 /* contador de la duracion de encendido.*/

bucle3:
sub.l #0x1,%d1
bne bucle3

move.w #0x3c,(%a0)

movea.l #0x320000,%a7
movea.l #0x310000,%a6

jsr start_main
Error:

move.w #0x00,(%a0)

forever:
jmp forever
----------------------------------------------------------------------
--------------------------- main.c ---------------------------------
#include "main.h"

void start_main(void)
{
LEDS_set_byte(0);

bienvenida();

mensaje();

}
----------------------------------------------------------------
---------------------- main.h -------------------------------
#define LEDS_ADDRESS 0x00100000
#define DIP_ADDRESS 0x00200000
#define RAM_BASE_ADDRESS 0x00300000
#define RAM_USER_BASE_ADDRESS RAM_BASE_ADDRESS+SYSTEM_INFO_LENGHT
#define RAM_LENGHT 0x00020000
#define LCD_DATA_ADDRESS 0x00400000
#define LCD_CONTROL_ADDRESS 0x00400001

#define SYSTEM_INFO_LENGHT 0x00000002

#include "func.h"
#include "LCD.h"
#include "LEDS.h"
#include "DIP.h"
#include "Int.h"

struct _SYSTEM_INFO
{
unsigned char LEDS;
};

void start_main(void);
-------------------------------------------------------------------
--------------------------- LEDS.c ----------------------------
#include "LEDS.h"
#include "main.h"

void LEDS_set_byte(unsigned char dato)
{
unsigned char *LEDS = (unsigned char*) LEDS_ADDRESS;
struct _SYSTEM_INFO * info = (struct _SYSTEM_INFO *) RAM_BASE_ADDRESS;

info->LEDS=dato;

LEDS_REFRESH(LEDS);
}

void LEDS_set_bit(unsigned char bit)
{
unsigned char *LEDS = (unsigned char*) LEDS_ADDRESS;
struct _SYSTEM_INFO * info = (struct _SYSTEM_INFO *) RAM_BASE_ADDRESS;

bit=1<

info->LEDS |= bit;

LEDS_REFRESH(LEDS);
}

void LEDS_clear_bit(unsigned char bit)
{
unsigned char *LEDS = (unsigned char*) LEDS_ADDRESS;
struct _SYSTEM_INFO * info = (struct _SYSTEM_INFO *) RAM_BASE_ADDRESS;

bit=1<
info->LEDS &= ~bit;

LEDS_REFRESH(LEDS);

}
--------------------------------------------------------------------
-------------------------- LEDS.h ------------------------------
#define LEDS_REFRESH(x) *x=info->LEDS

void LEDS_set_byte(unsigned char dato);
void LEDS_set_bit(unsigned char bit);
void LEDS_clear_bit(unsigned char bit);
--------------------------------------------------------------------
--------------------------- func.c -------------------------------
#include "main.h"

inline void prueba(unsigned char a);

/* espera(1) = 3.27315uS */
void esperar(unsigned long n)
{
do{}while(--n>0);

}

void bienvenida(void)
{
unsigned char *LEDS = (unsigned char *) LEDS_ADDRESS;
unsigned int p,q,dummy0, dummy1;
unsigned char k,j,i;

for(k=0;k<2;k++)
{
i=1;
for(j=0;j<8;j++)
{
LEDS_set_byte(i);
esperar(24000);
}
i=128;
for(j=0;j<8;j++)
{
LEDS_set_byte(i);
esperar(24000);
i=i>>1;
}
}

for(k=0;k<4;k++)
{
for(j=0;j<8;j++)
{

LEDS_set_byte(0);
LEDS_set_bit(j);
LEDS_set_bit(7-j);
esperar(48000);
}
}

prueba(j);
LEDS_set_byte(j);

}

__inline__ void prueba(register unsigned char a)
{
register unsigned long d3 asm("d3") = (unsigned long) a;

asm("move.b %d3,-(%sp) \n"
"jbsr LEDS_set_byte");
LEDS_set_byte((unsigned char) d3);

}

void mensaje(void)
{
LCD_configurar();
LCD_configurar();
LCD_configurar();
LCD_configurar();
LCD_configurar();
do
{
LCD_escribir_cadena("MC68HC000FN16",100000);
LCD_posicion(0,1);
LCD_escribir_cadena("Microprocessor",100000);
LCD_limpiar();
LCD_regresar();
}while(1);

}

---------------------------------------------------------------
------------------------- func.h ----------------------------

void esperar(unsigned long n);
void bienvenida(void);
void mensaje(void);

----------------------------------------------------------------

hay mas archivos pero luego los subo que estan malucos de pegar, por ultimo el Makefile
------------------- Makefile -----------------------------
#Compilador,Linkeo y generacion del bin para C.

all:
m68k-elf-gcc -ggdb -Os -m68000 -mshort -c func.c
m68k-elf-gcc -Os -m68000 -mshort -c -S func.c

m68k-elf-gcc -ggdb -Os -m68000 -mshort -c main.c
m68k-elf-gcc -Os -m68000 -mshort -c -S main.c

m68k-elf-gcc -ggdb -O0 -m68000 -mshort -c LCD.c
m68k-elf-gcc -O0 -m68000 -mshort -c -S LCD.c

m68k-elf-gcc -ggdb -Os -m68000 -mshort -c Int.c
m68k-elf-gcc -Os -m68000 -mshort -c -S Int.c

m68k-elf-gcc -ggdb -Os -m68000 -mshort -c LEDS.c
m68k-elf-gcc -Os -m68000 -mshort -c -S LEDS.c

m68k-elf-gcc -ggdb -Os -m68000 -mshort -c DIP.c
m68k-elf-gcc -Os -m68000 -mshort -c -S DIP.c

m68k-elf-as -m68000 -o boot.o boot.S

m68k-elf-as -m68000 -o Int_e.o Int_e.S

#m68k-elf-ld -dp -T ldscript boot.o Int_e.o main.o func.o LCD.o Int.o LEDS.o DIP.o /usr/m68k-elf/lib/gcc/m68k-elf/4.0.1/m68000/libgcc.a -o prueba.o
m68k-elf-ld -dp -T ldscript boot.o Int_e.o main.o func.o LCD.o Int.o LEDS.o DIP.o -o prueba.e

#m68k-elf-ld -dp -ggdb -T ldscript boot.o Int_e.o main.o func.o LCD.o Int.o LEDS.o DIP.o -o prueba.gdb

m68k-elf-objcopy -O binary prueba.e prueba.bin

./split prueba.bin

clean:
rm *.s *.o *.bin*
-------------------------------------------------
------------------- ldscript -------------------
ENTRY(_vectors)
SECTIONS {.text 0x0000 : {*(.text)}}
-------------------------------------------------

Basicamente es todo.

TODO:
hacer un TODO completo.

Tesis

2006-09-11T19:02:00.000-07:00

Dicen las malas lenguas que tengo que hacer tesis, y hasta el momento este va a ser mi proyecto, existe un uKernel llamado pistachio que lo han trabajado en Karlsruhe, la cosa es lograr hacer un pistachio distribuido que corra en varias maquinitas embebidas de un costo bien reducido.

De primerazo

2006-09-09T16:48:00.000-07:00

Obligado empienzo el blog, luego vere que le agrego.