Testing ChemDoodle Web Components

For more information check The website of ChemDoodle Web Components

A better way to do this with static HTML is probably loading the data into a hidden frame.

Lactococcus chaodoufuii

This strand of Lactococcus is surely funny. It's in Taiwanese with an accent.
Via Yakult central institute:

Exploring novel LAB strains from fermented brines used in producing stinky tofu (「臭豆腐」 調製用発酵液中の新規乳酸菌の探索)
Shiou Huei Chao*・富井 康明・笹本 正枝・工藤 裕子・藤本 淳治・渡辺 幸一・Ying-Chieh Tsai* (*台湾国立陽明大学)
　台湾の伝統的食品「臭豆腐」調製用の発酵液中の乳酸菌の多様性の解析により、新規乳酸菌を発見し、Lactobacillus capillatusおよびLactococcus chaodoufuiiと命名し、新種提案を行った。

Now it seems that phage might be back on stage?

To my surprise, almost no progress on the FDA side in 10+x years. Hard to believe that paradigm differences are so hugh.

Jackpot

甄試委員會：「這位同學請大概說明如果進入本所碩士班，你會想做些什麼？」
回答：「在大學三年間我在博士班學長的指導下，我參與了一些研究計畫，其中我們有發現了一個基因在很多物種之間的相似度幾乎百分之百，我想那一定是非常重要的基因，接下來我想繼續完成這樣的研究。」
.....
.....
It turns out that I wasn't kidding them. I lost the chance to continue to be involved in this project and for some reason we failed to publish what we found before others did. The gene was really close to the holy grail. The protein, a ribosomal protein, was found to play an important role recruiting mRNA into ribosome. I was so close to a gold mine.

本日最佳金句。

“We don’t teach alchemy along with chemistry.” - Dr. Francisco J. Ayala

Via new york times

Brain Science is Fun

Apparently the embed player is not working for this talk. Check out here at TED talk page for the video.
Dr. Jill Taylor described her observation on her own stroke as a neuroanatomist! Way better than analysis on functional MRI data by other projects!

Traveling Salesman Problem

Stop Comic Sans!

Ref: ban comic sans.

Using head tracking for molecular modeling

Understanding DGETRF/SGETRF interface

I was working on the inverse matrix the other day, it can be done through DGETRI of LAPACK. The DGETRI requires DGETRF to do LU decomposition first in order to calculate the inverse matrix. Then I found DGETRF interesting to start with since it has been awhile for me to touch linear algebra stuff. I wrote a "hello world" code for it. I was stuck on the issue to construct permutation matrix from IPIV, basically IPIV is not like vecP in Mathematica. This is how I understand it, I am not confident about the P matrix, this requires some LAPACK friends' help. On Mac OS X, just do g95 -framework Accelerate hello.f90 -o hellolapack and run the program by ./hellolapack , sorry I am too lazy to install a software lapack. :p

program hello
   integer INFO
   integer IPIV(3)
   real*8 A(3,3),B(3,3),C(3,3),D(3,3)
   real*8 P(3,3),L(3,3),U(3,3)
   real*8 myVEC(3)
   real*8 WORK(64*3)
   A(1,1)=1d0
   A(1,2)=2d0
   A(1,3)=2d0
   A(2,1)=2d0
   A(2,2)=7d0
   A(2,3)=7d0
   A(3,1)=2d0
   A(3,2)=7d0
   A(3,3)=9d0
   P=0d0
   P(1,1)=1d0
   P(2,2)=1d0
   P(3,3)=1d0
   ! dgetrf(M,N,A,LDA,IPIV,INFO)
   ! input
   ! M          The number of rows of the matrix A.  M >= 0. Integer
   ! N          The number of columns of the matrix A.  N >= 0. Integer
   ! A          The DP array
   ! output
   ! A          The DP array
   ! LDA        The leading dimension of the array A.  LDA >= max(1,M).
   ! IPIV       INTEGER array, dimension (min(M,N))
   !            The pivot indices; for 1 <= i <= min(M,N), row i of the
   !            matrix was interchanged with row IPIV(i).
   print *,"A"
   write (*,'(F8.3,F8.3,F8.3)') A(1,1:3)
   write (*,'(F8.3,F8.3,F8.3)') A(2,1:3)
   write (*,'(F8.3,F8.3,F8.3)') A(3,1:3)
   B=A
   call dgetrf(3,3,A,3,IPIV,INFO)
   if (INFO.ne.0) stop
   !IPIV(1)= 3 implies that the first and third rows were interchanged
   !when factoring the first column; IPIV(2)= 3 implies that the second
   !and third rows were interchanged when factoring the second column.
   !In this case, IPIV(3) must be 3 because there are only three rows.
   if (IPIV(1).eq.3) then
      myVEC=P(1,1:3)
      P(1,1:3)=P(3,1:3)
      P(3,1:3)=myVEC
   elseif(IPIV(1).eq.2) then
      myVEC=P(1,1:3)
      P(1,1:3)=P(2,1:3)
      P(2,1:3)=myVEC
   endif
   if (IPIV(2).eq.3) then
      myVEC=P(2,1:3)
      P(2,1:3)=P(3,1:3)
      P(3,1:3)=myVEC
   elseif(IPIV(2).eq.1) then
      myVEC=P(2,1:3)
      P(2,1:3)=P(1,1:3)
      P(1,1:3)=myVEC
   endif
   print *,"P"
   write (*,'(F8.3,F8.3,F8.3)') P(1,1:3)
   write (*,'(F8.3,F8.3,F8.3)') P(2,1:3)
   write (*,'(F8.3,F8.3,F8.3)') P(3,1:3)
   print *,"L"
   L(1,1)=1d0
   L(1,2:3)=0d0
   L(2,1)=A(2,1)
   L(2,2)=1d0
   L(2,3)=0d0
   L(3,1:2)=A(3,1:2)
   L(3,3)=1d0
   write (*,'(F8.3,F8.3,F8.3)') L(1,1:3)
   write (*,'(F8.3,F8.3,F8.3)') L(2,1:3)
   write (*,'(F8.3,F8.3,F8.3)') L(3,1:3)
   print *,"U"
   U(1,1:3)=A(1,1:3)
   U(2,1)=0d0
   U(2,2:3)=A(2,2:3)
   U(3,1:2)=0d0
   U(3,3)=A(3,3)
   write (*,'(F8.3,F8.3,F8.3)') U(1,1:3)
   write (*,'(F8.3,F8.3,F8.3)') U(2,1:3)
   write (*,'(F8.3,F8.3,F8.3)') U(3,1:3)
   call dgemm("N","N",3,3,3,1d0,L,3,U,3,0d0,C,3)
   print *,"L.U"
   write (*,'(F8.3,F8.3,F8.3)') C(1,1:3)
   write (*,'(F8.3,F8.3,F8.3)') C(2,1:3)
   write (*,'(F8.3,F8.3,F8.3)') C(3,1:3)
   call dgemm("N","N",3,3,3,1d0,P,3,L,3,0d0,C,3)
   call dgemm("N","N",3,3,3,1d0,C,3,U,3,0d0,D,3)
   print *,"P.L.U"
   write (*,'(F8.3,F8.3,F8.3)') D(1,1:3)
   write (*,'(F8.3,F8.3,F8.3)') D(2,1:3)
   write (*,'(F8.3,F8.3,F8.3)') D(3,1:3)

   call dgetri(3,A,3,IPIV,WORK,64*3,INFO)
   if (INFO.ne.0) stop
   call dgemm("N","N",3,3,3,1d0,B,3,A,3,0d0,C,3)
   print *,"A.A^-1"
   write (*,'(F8.3,F8.3,F8.3)') C(1,1:3)
   write (*,'(F8.3,F8.3,F8.3)') C(2,1:3)
   write (*,'(F8.3,F8.3,F8.3)') C(3,1:3)
end program hello

世界最大火力發電廠

Kornberg (1918-2007) 去世

Requiescat in pace.

Bacterial flagellum

- via random flickr search

Entropy and Time

Just joking

One of my random thought in 1998 after a debate on SciFi issues:

整體 △Ｓ 會增加定義了時間的方向。
^^^^^^^^^^^^^^^^      ^^^^^^^^^^
      ｐ          →      ｑ
所以非 p 無法得到非 q
    熵流失       時間倒流

整體熵快速的流失的狀況無法推得時間在倒流。

我的物化是不及格了三次的。 (The posting date is an estimation.)