My other computer is a 4000 unit Beowulf cluster

Mais um passo na minha implementação das funções de verossimilhança empírica generalizada - agora elas estão paralelizadas, e o procecimento de bootstrap pode rodar automaticamente distribuindo o processamento nos cores da maquina ou então em um cluster.
Isso é especialmente util para obter propriedades em amostras finitas usando Monte Carlo, ou então como já comentei obter intervalos e testes via bootstrap.
Mas quero ir um passo além e paralelizar mais partes do algoritmo - em especial os procedimentos de algebra linear e as funções de otimização. Nesta situação mesmo uma única estimação seria realizada de forma paralela.
Não é tão complicado, mas é um tanto trabalhoso já que preciso modificar os códigos fontes em c/c++ das biblioteca de otimização que eu uso. Já testei para um das funções e funciona bem.