Acerca de:

Este blog contiene los códigos, ejemplos y bases de datos que he usado cuando aprendía acerca de algún tema específico. En lugar de borrarlos (una vez dominado ya el tema), he decidido publicarlos :)

viernes, 12 de septiembre de 2014

Cómo ejecutar Visual Studio 6 en Windows 7 y evitar los errores con referencias y dependencias

Tengo un Visual Basic 6 Enterprise corriendo en Windows 7 de 32 bits. Todo parece funcionar bien: se ha instalado sin problemas, y se ejecuta sin inconveientes.
Pero al abrir un proyecto nuevo del tipo Exe estándar, e intentar cargar librerías y ocx de terceros (no son de Microsoft, o no están firmadas por Microsoft) me salta este error "Error al intentar tener acceso al registro del sistema":


(No lo parece, pero esto es un Windows 7, sólo que con el tema de Windows Clásico. Los colorcitos sólo sirven para consumir memoria, otra razón es que así puedo aluscinar que tengo un Windows Server)


Este fue un error muy feo que me tuvo ocupada durante horas (además de hacerme instalar y desinstalar las ocx varias veces), y al final la solución era muy simple, ejecutar el Visual Studio 6 como Administrador:


Ya con todos los permisos disponibles, el Windows nos deja importar componentes, referencias y dependencias. Y como se muestra en la imagen, el Visual Studio 6 puede convivir perfectamente con el Visual Studio 2013 (yo tengo la versión Express).
Un último consejo: si se tiene un proyecto de Visual Basic 6 con dependencias de terceros, es mejor ejecutar primero el IDE como administrador, y luego abrir nuestro proyecto. Se evitarán así muchos errores.

domingo, 24 de agosto de 2014

Algoritmos para calcular la factorización QR de una matriz

La factorización QR se utiliza para descomponer una matriz (que puede representar un sistema de ecuaciones a resolver) en dos matrices, con las cuales es más fácil de operar y hallar soluciones al sistema de ecuaciones.
Más información en la WIkipedia.

Los siguientes códigos están basados en los algoritmos mostrados en el libro "Numerical Linear Algebra" de Trefethen y Bau.

Algoritmo de Gram-Schmidth Modificado (devuelve la factorización QR reducida):



El código en Matlab es:

function [Q, R] = gs(A)
[m, n] = size(A);
R = zeros(n, n);
V = A;
Q=zeros(m, n);

for i =1:n
    R(i,i)= norm(V(:,i));
    Q(:,i)= V(:,i)/R(i,i);
   
    for j=i+1:n
       R(i,j)= (Q(:,i)')*V(:,j);
       V(:,j)=V(:,j) - R(i,j)*Q(:,i);
    end
end

"A" es una matriz de rango completo, la función devuelve dos matrices: Q (la matriz ortonormal) y R (la matriz triangular).

Problemas con este algoritmo: Sólo da la factorización QR reducida, y sólo funciona para matrices de rango completo (es decir: ninguna columna es múltiplo de, y/o resultado de, sumar otras columnas de la mariz entre sí).

Modifiqué este algoritmo para que funcionara también con matrices de "rango deficiente" (lo opuesto a "rango completo"), y que devolviera la factorización completa; esto último se hace completando las columnas de la matriz Q, las cuales son ortogonales entre sí. En palabras del curso de Álgebra Lineal: "completando la base ortonormal del espacio vectorial generado por las columnas de Q". Dicho de otra forma: aplicando el Teorema de la Base Ortonormal Incompleta.

El código en Matlab es:

function [Q, R] = gsC(A)
clc;
[m, n] = size(A);
R = zeros(m, n);
V = A;
Q=zeros(m, n);

k=1; % inicializo eta variable, sólo la uso para matrices de rango deficiente

for i =1:n
    R(i,i)= norm(V(:,i));  
   
    if abs(R(i,i)) > 1.0e-10
        Q(:,i)= V(:,i)/R(i,i);
    else  % salió cero! La matriz A es de rango deficiente
        r= zeros(1, m); % este va a ser el vector que reemplazará a la fila problemática de A
              
        Q(:,i)= r;
        Q(k, i)=1;
        k=k+1;
        Q2 = Q; % Q2 tendrá los valores de Q, en Q irá la nueva base ortonormal

        % ortogonalizo otra vez aplicando Gram-Schmidt a la matriz Q para
        % hallar un vector con el qué reemplazar a la fila problemática y
        % seguir teniendo una base ortonormal
        for j = 1 : i-1
            r(j) = (Q2(:,j)')*Q2(:,i);
            Q(:,i)=Q(:,i) - r(j)*Q2(:,j);
            Q(:,i)=Q(:,i)/norm(Q(:,i));
        end
    end
   
    for j=i+1:n
       R(i,j)= (Q(:,i)')*V(:,j);
       V(:,j)=V(:,j) - R(i,j)*Q(:,i);
    end
end

% hasta aquí Q y R son la descomposición reducida
% acá hallo la completa, Gram Schmidt otra vez
if m > n   
    I = eye(m, m - n);
    Q=[Q I];
    V=Q;
    r= zeros(1, m);
   
    k=0;
    for i=n+1:m
        for j = 1 : n + k
            r(j) = (V(:,j)')*V(:,i);
            Q(:,i)=Q(:,i) - r(j)*V(:,j);
            Q(:,i)=Q(:,i)/norm(Q(:,i));
        end
        k=k+1;
    end
end

"A" es una matriz cualquiera, la función devuelve dos matrices: Q (la matriz ortonormal) y R (la matriz triangular). 

Otra forma de calcular la factorización QR (y la que utiliza internamente Matlab) es usando la factorización de Householder:


El código en Matlab es:

function [Q, R] = houseHolder(A)
clc;
[m n] = size(A);

R= A;
Q = eye(m,m);

n1=min(m,n);

for i= 1:n1
    x = R(i:m, i);
    s=0;
   
    if x(1)==0
        s = 1;
    else
        s= sign(x(1));
    end
   
    e = zeros(m - i + 1, 1);   
    e(1)=1;   
    u = s * norm(x) * e + x;
    u = u / norm(u);
    R(i:m, i:n) =  R(i:m, i:n) - 2*(u*u')*R(i:m, i:n);
    Q(:,i:m) = Q(:,i:m) - Q(:,i:m)*(2*u*u');   
end


Igual que en los casos anteriores, "A" es una matriz cualquiera, la función devuelve dos matrices: Q (la matriz ortonormal) y R (la matriz triangular).

Los tres algoritmos han sido probados en Matlab 2008a, así que deben correr en cualquier versión posterior, y no tan posterior :)

miércoles, 9 de julio de 2014

Crystal Reports Runtime: No se ha podido cargar el informe

El error completo es:
"CrystalDecisions.Shared.CrystalReportsException: No se ha podido cargar el informe. Un documento procesado por el motor JRC no se puede abrir en la pila C++".
 En concreto, salta cuando se utilizan las librerías de CrystalReports para Net.


Es un error muy feo para el motivo que lo produce: el archivo rpt que debe cargarse simplemente no está en la ruta que le indicamos. Exactamente en el método Load del objeto ReportDocument:

Private WithEvents oRpt As ReportDocument
oRpt = New ReportDocument
oRpt.Load("../informes/CrystalReportAutores.rpt")


Código sacado de la Güeb de Joaquín :D
Agradezco a este enlace por darme la solución :D

martes, 3 de junio de 2014

Npgsql lanza una System.OutOfMemoryException

Jugaba un día con el C#, las librerías del Npgsql y con una base de datos con dos tablas vinculadas entre sí por una llave foránea. Luego de insertar unos cuantos registros en ambas tablas, intenté borrar los registros de una de ellas y entonces saltó esta excepción:


Saltó una "System.OutOfMemoryException". Llegué a pensar que mi programita tendría fugas de memoria en algún lado, pero al intentar borrar los mismos registros desde el pgAdmin III saltaba el verdadero error: violación de las reglas de la llave foránea.

Me pareció curioso así que, luego de experimentar un poco, descubrí que el Npgsql lanza la System.OutOfMemoryException cada vez que se viola una relación de llaves foráneas (en mi caso: intentar borrar los registros cuyas llaves primarias estaban vinculadas hacia la otra tabla por medio de la llave foránea).
En general, esta excepción suele saltar si se violan las reglas establecidas para las secuencias (datos de tipo serial), llaves primarias y llaves foráneas.

Pero me quedaba una pregunta: ¿no debería saltar una InvalidOperationException?

Mi teoría es la siguiente: una llave foránea es en realidad un puntero hacia los datos de otra tabla. Al querer realizar una operación no permitida (por ejemplo intentar borrar los datos a donde apunta el puntero pero sin borrar el puntero) se tiene un puntero que apunta (valga la redundancia) hacia un área de memoria que ahora es inválida. Para colmo, el sistema operativo (si es código no administrado no hay colector de basura) no puede reasignar esa memoria para su uso posterior ni liberar el puntero. En resumen: resulta una fuga de memoria, la cual se traduce (al pasar ya al mundo del código administrado del .Net) como una OutOfMemoryException.

Para saber más del tema de los punteros y sus problemas, recomiendo este pdf.

viernes, 9 de mayo de 2014

¿Cuál es la diferencia entre “Null” y “Nothing” en VB6?

Esto lo pregunté en StackOverflow hace un tiempo. Acá pongo la traducción:

Tengo un recordset como el siguiente:


Dim rs as Recordset
Set rs as New Recordset

'... a lot of coding ...

if Err.Number <> 0 Then ' oops, algo salió mal
    If rs.State <> adStateClosed Then rs.Close
    Set rs = Nothing
end if

' quiero evaluar si rs es Nothing o Null

if rs is Nothing then 
' esto no tira errores y funciona bien :D
end if

if rs is Null then
' esto tira un error de "tipos no compatibles"
end if

if rs = Null then
' esto tira un error de "tipos no compatibles"
end if

if isNull(rs) then
' nunca entra aquí, isNull(rs) evalua a Falso
end if 

La respuesta que me dio Mark Bertenshaw es:
"Null" es un subtipo específico de Variant. No existe fuera de Variant, y se creó para permitir modelar un dato de tipo nulo en una base de datos.
"Nothing" es un valor de tipo Object. Identifica un puntero nulo, es decir: no hay objeto.
Al evaluar un recordset con el operador "is" tira un error porque sólo se puede usar con datos de tipo Object.
IsNull(rs) retorna falso porque evalúa variables del tipo Variant.
IsNull(rs) equivale a VarType(rs) = VBNull

sábado, 22 de marzo de 2014

Dibujando matrices con Latex

El código básico para dibujar una matriz lo conseguí de este enlace. Pero yo necesitaba dibujar una matriz como la siguiente:


El código en Latex para dibujarla es:

\[ \left(
\begin{array}{cccc}
b\textsubscript{11} & b\textsubscript{12} & ... & b\textsubscript{1n} \\
b\textsubscript{11} & b\textsubscript{12} & ... & b\textsubscript{1n} \\
\vdots & \vdots & \vdots & \vdots \\
b\textsubscript{n1} & b\textsubscript{n2} & ... & b\textsubscript{nn}
\end{array}
\right)\]



Luego tenía que dibujar una matriz como combinación de columnas multiplicadas por un escalar, como la siguiente:

El código es:

\[b= \left(
\begin{array}{c|c|c|c|c}
&&&& \\
&&&& \\
u\textsubscript{1}\mbox{. \={v}} & u\textsubscript{2}\mbox{. \={v}} & u\textsubscript{3}\mbox{. \={v}} &...& u\textsubscript{n}\mbox{. \={v}}\\
&&&& \\
&&&&
\end{array}
\right)\]


jueves, 20 de febrero de 2014

Continuando con el Conjunto de Wirth...

Actualización de Electrónica a Martillazos, hallé algunas cosas interesantes con respecto al algoritmo de PuntoPeek para calcular el Conjunto de Wirth.
El enlace aquí.
Un mirror aquí.

lunes, 13 de enero de 2014

Resetear campo autoincrementable (serial) en Postgre 9.3

Tengo una base de datos llamada "Prueba" con una tabla "phonebook", la cual tiene un campo autoincrementable (tipo serial) llamado "ID":


Cuando se crea un dato serial usando el pgAdmin, Postgre automáticamente nos crea una secuencia (en este caso llamada "phonebook_ID_seq") a la cual el campo serial está vinculado mediante la siguiente instrucción:

ALTER TABLE phonebook ALTER COLUMN "ID" SET DEFAULT nextval('"phonebook_ID_seq"'::regclass);

Para Postgre, un dato serial es un dato Integer con la función "nextval".

Esta es mi situación:
Luego de añadir y borrar varios registros a la tabla "phonebook", y finalmente dejarla vacía, quiero resetear el valor del campo "ID" a 1.

Al igual que en MySql, Postgre tiene una función llamada "setval" que me permite ponerle un valor arbitrario a la secuencia desde el cual continuará con los incrementos.
En el pgAdmin seleccionamos el campo "ID", abrimos la ventana de queries mediante el botón "Execute Arbitrary Queries" (el que tiene una lupa y el texto "SQL"), e ingresamos el siguiente query:

select setval('"phonebook_ID_seq"',1);

Hay que notar que lo que se altera es la secuencia, no el campo "ID". Hay que tener cuidado de poner el nombre de la secuencia entre comillas dobles dentro de las comillas simples. Como Postgre es Case Sensitive, saltarán errores si omitimos las comillas dobles y si el nombre de la secuencia contiene mayúsculas.

En la ventana de queries quedará así:


Si todo sale bien, nos aparecerá el valor al que hemos reseteado la secuencia. No tiene que ser 1, puede ser cualquier valor que queramos.

miércoles, 8 de enero de 2014

Postgre 9.3: Errores 42501 y 42703:

Tengo la siguiente base de datos Prueba instalada en un servidor local de Postgre versión 9.3, la cual tiene una sola tabla "phonebook" cuya estructura copié de algún sitio de Internet (al instalar Postgre se tiene la opción de instalar también el pgAdmin):


El Owner de la base de datos y de la tabla es el superusuario "postgre". La base de datos Prueba tiene como esquema por defecto (Default Schema) igual a público (public), también sus privilegios son públicos.

También he creado un usuario "user01" con las siguientes propiedades:


Al conectar una aplicación en C# mediante el usuario "user01", usando Npgsql (el cual se puede instalar junto con el Postgre usando la aplicación Stack Builder que también viene con el instalador del Postgre) y luego de ver los datos de la tabla phonebook todo iba bien; pero cuando quería modificar la tabla phonebook con un update, insert ó delete me arrojaba el siguiente error:


"Error 42501: permiso denegado a la secuencia phonebook_ID_seq"

Este error salta si el usuario con el que intento modificar la tabla no es su dueño (owner). No interesa si la base de datos tiene como dueño al superusuario "postgre", la tabla debe tener como dueño al usuario que intenta modificarla.

Para cambiar al dueño de la tabla phonebook de "postgre" a "user01" basta este script:

ALTER TABLE phonebook
  OWNER TO user01;

Y entonces al usar la siguiente sentencia insert desde C#:

string query = @"insert into phonebook (phone, firstname, lastname, address) values ('1234', 'Wert', 'Asdlo', 'Fake Street 123')";

Los datos se guardaban correctamente, pero al querer actualizar o borrar: 

string query1 = @"Update phonebook set phone ='0000', firstname ='Carry', lastname ='Doe', address = 'No signed' where ID = 1";
string query2 = "delete from phonebook where ID = 1";

Me saltaba el siguiente error:


"Error 42703: no existe la columna id"

Lo extraño es que la tabla phonebook sí tiene una columna llamada "ID". Buscando por internet, hallé aquí la solución: se debe poner el nombre del campo que contiene mayúsculas entre comillas:

string query1 = "Update phonebook set phone ='0000', firstname ='Carry', "
" lastname ='Doe', address = 'No signed' where \"ID\" = 1";
string query2 = "delete from phonebook where \"ID\" = 1";

El caracter \ le indica al C# que la comilla que le sigue no indica el final de la cadena. El motor del Postgre convierte todos los nombres de las columnas en una sentencia sql a minúsculas. Por ahí leí que, internamente, Postgre es Case Sensitive.