翔〜　望♭

progress tracking

2010-03-24T15:41:00.000-07:00

created a script to automatically generate progress bars

yay.

spring break

2010-03-22T05:32:00.000-07:00

ambitious ( = i won't be surprised if only 50% gets done) spring break goal

1. Brahms

- mvt 1

- mvt 2

- mvt 3 (low priority)

- mvt 4

2. Elasticity

- controlled bending energy, derivation and implementation

- bending experiments

- re-implement tet mesh + minimal topology code

- re-implement elasticity code on tets

- re-implement dynamics on tets

- re-implement geodesic interpolation on tets

- re-implement parameterization idea on triangle meshes

3. Read The Pragmatic Programmer

4. Project Ashita

- project vbw

- project 4st

5. get cramster job @_@

- and start $_$ing

6. table tennis

7. IST lunchbunch

- abstract

- slides

what does this function do?

bool lolcopter(double x) {
　　return x != x;
}

image surgery

2010-03-21T19:35:00.001-07:00

original
seam carve
seam carve
seam carve
gimp
gimp

kansei

Statistically Significant

2010-03-09T18:19:00.000-08:00

(06:17:20 PM) yjiatwo: uuuuuuuuuuuuuuuuuuuuuuuuyjjjjjjjjjjjjjjjjjjjjjjjjjjjjjjjjjjjjhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh
(06:17:29 PM) yjiatwo: i was cleaning keyboard
(06:17:32 PM) yjiatwo: dropped some cheese on it
(06:17:35 PM) me: lol
(06:17:43 PM) me: uhujuhjjjjjuhujh
(06:17:53 PM) me: i was simulating cleaning that batch of keys
(06:18:07 PM) me: i find that your exaggeration is statistically significant
(06:18:24 PM) yjiatwo: your simulation seems to be wrong

Engrish

2010-03-06T18:59:00.000-08:00

(06:20:11 PM) yjia2: if my techcomm teacher reads the first paragraph he'll be like EF EFEF YOUR ENGLISH SUCKS
(06:20:18 PM) me: uh...
(06:20:27 PM) me: i dont need your techcomm teacher to tell me that...

i've been speaking Engrish (and not English) my entire life, hard to change that fact, unfortunately..

Anyway, doing election fraud project... why were the LDP missing so many 4s in 1996?
4にたくないから?

Also trying out machine learning ideas for fraud detection, seems to be working out great..

=)

Simplified/Traditional Chinese -- cussing or not?

2010-02-28T00:41:00.000-08:00

Too often you see Engrish.com posting something to the effect of
"干货" translating into "fucking goods"

That's ok -- it is funny at the level where people are lazy/cheap enough to just throw text into online translators.. and the translator recognizes that 干->幹->fuck, while it is really 干->乾->dry

so they don't mean to cuss at you when they use 干

However, when you see it in traditional, no excuse, they are cussing at you.

I didn't order fucking noodles from them.

Buzz

2010-02-10T17:41:00.001-08:00

From now on you can use "twitter" inside the gmail interface
Not only are you allowed to delete posts, you can edit them..

for instance you can do something like this..
originally talking about taking a midterm exam

... to something else...

have fun =)

IT ojisan

2010-02-09T23:22:00.001-08:00

Latest comic from a Japanese programming magazine, mentioned in a previous post --
with *fansubs* lol,

enjoy.. perhaps someone can digg this and rankmaniac 2010 rankmaniac 2010 it up

Full listing of the original comic series in Japanese here

@codezine.jp: contact me if you wish this comic be taken off..

Viewer overhaul

2010-02-08T11:50:00.000-08:00

rewriting all my mesh viewer code ...
implemented arcball in C++.. going to re-implement topology construction as the next step.

mal yeliah is also a Rankmaniac 2010!
patrick is still a Rankmaniac 2010!
I'm not a Rankmaniac 2010, though.

By the way, i have a twitter account, which is very fun (and easy) to follow.. you should look at it

By the way #2, project euler is getting more popular in Taiwan .. she just found its 101st participant!

Taiwan 101!

CS2 Robot Arm

2010-02-07T01:19:00.001-08:00

This is the second assignment i wrote for cs2..
robot arm's been in the cs2 curriculum, but in java.
With this year's aim to overhaul cs2 and to switch to python, everything has to be rewritten.

The basic framework wasn't much of a problem, since I have rewritten the same assignment for cs171 last term already. The problem is, how to make it more understandable and -- if you will -- more hand-holding for the students?

I ended up adding several components that will check their answers, so they can have more milestones..

.. but that sounds funny.. if i check their answers, that means the code we provide must compute the answer! Yes it does, but since I can write obfuscated code without much effort, that wasn't much of a problem. Actually, that's not the reason why I dont worry -- the answer is computed implicitly using some opengl tricks that i doubt anyone will bother decoding..

anyway that took forever to write, but its a nice assignment with key concepts

Continuation methods, pseudoinverses, finite differences, and rankmaniac 2010

rankmaniac 2010

2010-02-05T14:21:00.000-08:00

Search for "rankmaniac 2010" on google/bing..
you'll see some heated war among caltech cs/ee students..

of course i'm not involved .. aside from the fact that i publicly support Patrick's Rankmaniac 2010 conquest.

I am NOT a Rankmaniac 2010.

Running on the clock

2010-02-05T11:02:00.000-08:00

How do you compute the last three digits of 2007^2008^2009^2010.... by hand?

even python cannot do it without a little bit help from you ..
>>> pow(2007, pow(2008, pow(2009,2010)), 1000)
zzz

A few tools we will need...
Modular arithmetic (同餘)
Euler's Totient Theorem (歐拉定理)
a^phi(n) = 1 mod n
for gcd(a, n) = 1
The Chinese Remainder Theorem (中國餘數定理)
x = a_1 mod n_1
x = a_2 mod n_2
and gcd(n_1, n_2) = 1,
then there's a solution x = y mod (n_1 * n_2) satisfying the two above constraints.

to start, when we say we want to compute the last 3 digits of something, we are computing mod 1000. Why? because when we find the remainder from dividing something by 1000, we are filtering out eveything but the last three digits :
x = 182412somthingsomething1823012312 * 1000 + 231 has a remainder of 231 when divided by 1000, for example.

so the problem becomes 2007^2008^2009^2010 mod 1000

2007^2008^2009^2010 mod 1000
= 7^2008^2009^2010 mod 1000

This is by simple rules of modular arithmetic, consider

(2*1000 + 7)^a mod 1000 and look at its binomial expansion:
(2*1000)^a + a * (2*1000)^(a-1)*7 + ... + a * (2*1000)*7^(a-1) + 7^a
all term but the last have a factor of 1000, and are hence modded out.

next step..
gcd(7, 1000) = 1, so we can apply euler's theorem.
phi(1000) = 400 (very easy to remember.. phi(10) = 4, phi(100) = 40, ... so on)

so
2007^2008^2009^2010 mod 1000
= 7^2008^2009^2010 mod 1000
= 7^(2008^2009^2010 mod 400) mod 1000
= 7^(8^2009^2010 mod 400) mod 1000

why is that?
from euler 7^phi(1000) = 7^400 = 1 mod 1000
if the exponent is A = 400 * m + (A mod 400)
then
7^A mod 1000
= 7^(400 * m + (A mod 400)) mod 1000
= (7^400)^m * 7^(A mod 400) mod 1000
= 1^m * 7^(A mod 400) mod 1000
= 7^(A mod 400) mod 1000

So now let's first compute

A = 8^2009^2010 mod 400

(btw, once you have reduced your problem to this size, python can solve it for you already:
for those who want to do the computation by hand, do not run this python script

>>> pow(7, pow(8, pow(2009,2010), 400), 1000)
answer here
)

We invoke Chinese remainder theorem here to compute the exponent A

first, our base is 8, which is not relatively prime to 400. Factor 400 = 16 * 25
and we know
A = 8^2009^2010 = 0 mod 16 (since there are wayyyy too many factors of 2)
and
A = 8^2009^2010 = y mod 25.

once we find y, we can use the Chinese remainder theorem to get A mod 400.

So first we need to compute
8^2009^2010 mod 25

gcd(8, 25) = 1, compute phi(25) = 20, (and phi(20) = 8)

8^2009^2010 mod 25
= 8^(9^2010 mod 20)
= 8^(9^(2010 mod 8) mod 20) mod 25
= 8^(9^2 mod 20) mod 25
= 8^(81 mod 20) mod 25
= 8^1 mod 25
= 8

so

A = 8^2009^2010 = 8 mod 25
and
A = 8^2009^2010 = 0 mod 16

or in other words, A is a multiple of 16, and A is 8 more than another multiple of 25

A = 16*m = 25*n + 8

m = (25n + 8)/16 = n + (9n + 8)/16

what form must n take so that m is guaranteed to be an integer?
n is an integer, so for m to be integral, so must (9n + 8)/16

we can find n to be
n = 16t + 8,
since
9*(16t + 8) + 8 = 16*(9t) + 72+8 = 16*9t + 16*5 is divisible by 16.

A = 25*n + 8 = 25*(16t + 8) + 8 = 400t + 208
so
A = 8^2009^2010 = 208 mod 400

then all that's left is to compute

7^208 mod 1000

... not much you can do here but to brute force

= (7^2)^104
= 49^104 mod 1000
= ((50 - 1)^2)^52 mod 1000
= (2500 - 100 + 1)^52 mod 1000
= 401^52 mod 1000
= ((400 + 1)^2)^26 mod 1000
= (160000 + 800 + 1)^26 mod 1000
= 801^26 mod 1000

(side note: to do 801 * 801 mod 1000,
just do regular multiplication, but only look at last three digits
801
801
-------
801
00
8
-------
601

)

= 601^13 mod 1000
= 601 * (601^2)^6 mod 1000
= 601 * (201^2)^3 mod 1000
= 601 * 401^3 mod 1000
= 601 * 401 * 401^2 mod 1000
= 601 * 401 * 801 mod 1000
= 801 mod 1000

last three digit of 2007^2008^2009^2010 is 801

(we didn't really have to go through all that work, since there's some nice property/pattern hidden in there that we could have used. Nonetheless that was a great exercise using number theoretical tools...)

日本語から... English

2010-02-01T03:29:00.000-08:00

Yea.. got tired of using Japanese.. 其實又有點想用中文...
Anyway, found this interesting IT related manga series online.. but one problem is that they didn't index their manga in an orderly fashion like XKCD.....

So by design, they always have a pointer to the previous manga, but never a forward pointer, because they have no clue a priori what the posting number is gonna be...
The best one can do is to read backwards using their "previous" button.. but who wants to read backwards? So hence the journey of grep-ing for all the manga..

First, they do have a page that contains all the article number to the manga series, so i won't have to read through all the articles to search for them..
All that remains then is to open up those article html and grep for images... but again failed design.. the image names do not follow a pattern..

oh well, do the best i can and probably do some manual filtering..

on the terminal, run a single(?) command:

base="http://codezine.jp"; for i in `seq 1 4`; do for page in `wget -O - -q $base"/article/corner/210?p=$i&" | egrep -o -i /article/detail/[0-9]+ | uniq -`; do echo $base`wget -O - -q "$base$page" | egrep -o -i "/static.*[0-9]+\.jpg"` >> itlist.txt; done ; done

that will throw all the candidate manga image links into itlist.txt

(now do some manual filtering -- which turns out to be fairly minimal)

with a clean list of image links, I can get a javascript-formatted list by running a python one-liner

print [line.strip() for line in file("itlist.txt")][::-1]

recall that [::-1] is the python idiom for getting a reversed copy of a list.

Review:
L[a:b:c] is getting a sublist of L from a (inclusive) to b (exclusive) every other c elements. Using the default value for a and b will give you the full range, and using c = -1 reads the list backwards.

so yea, there i have it, a list of the manga images and i can now read forwards!

罰金

2009-07-14T21:55:00.000-07:00

先月　ケイタイの通話制限を越えたから、何百円も請求された

composition いっぱいつかった

cubes = takeWhile (<10^13) . dropWhile (<10^9) $ map (^3) [1..]
l = filter ((5==).length) . groupBy (\(x,_) (y,_)->x==y) .
　　　sort . map (\x->((sort . show) x, x)) $ cubes
ans = head . sort . map snd . concat $ l

PETSc + TAO が強い

カン

2009-07-12T01:00:00.000-07:00

近所のスーパーには自動リサイクル機があるのを気付いた
高校時代の友達と一緒に行って、コイン返却口がないのも気付いた
気になったから、プラスチック缶を探すために周りをくまなく探した
最初はどこにも見られなかったが
　やっと大きいゴミ箱で一本を探した
ゴミを見ただけでそんなに楽しくなれること、知らなかった
やってみて、レシートが出て、それで店員にお金を要求するというわけだ

project euler 197, 204　を解けた
197は無限リストを使った

u :: [Double]
u = -1 : (map f u)

f :: Double -> Double
f x = 1E-9 * m
　　　where b = 30.403243784 - x^2
　　　　　　　m = fromIntegral $ floor (2 ** b)

/_\

2009-06-30T07:23:00.001-07:00

夏休みに家でゆったり過ごしてるのに
頭痛の一日だった

作文

2009-06-09T15:04:00.001-07:00

済まない、君を誤解した

永遠な仲間　ワン・ジャックソン　に捧ぐ

　人生には不如意なことがあるが、それこそが人生だ。一生懸命頑張ったのに、なお上手くできないこともあるし、他人の成功を羨んで自分に不満するし、どうにかして「なぜ神様は私に不幸を与えるだろうか」と思ってしまう。そうすると意気消沈したり、人生への夢を諦めたりしてしまう。ただし、自分のことであるからこそ、責めをほとんど他人に負わせるわけにはいかないだろう。自分だけではなく、異なる問題でも、みんな同じ課題に直面している。挫折に遭ったいま、頽廃するか、それとも逆境に対して何かをするか自分次第だ。
　
　挫折は必ずしも悪いことだ。挫折がなかったら、まるで自分ができた人で改善する必要のあるところもないように、向上心もなくしてしまう。「失敗から学ぶ」ということわざを多分知っているが、よほど不如意なことを失敗のように扱わないと適用できないと思ってしまう。挫折はまるで転げ落ちるようだ。落ちるともちろん立ち上がらなければならない。しかし、それだけすれば、本当に骨折り損のくたびれもうけだ。いたずらに傷つき、学ぶ機会も無駄にしてしまう。転げ落ちたら、なぜ落ちた原因を探さないだろうか。できるだけどこが間違っていたのをのかを調べる。もちろん、自分が問題があるわけがないけれど、常に学べるものがあるだろう。例えば、自分の何かに不満しているとき、他人の意見を求めば、それがそんなに悪くないと教えてくれるかもしれないし、自分が見過ごす長所も教えてくれるんじゃないか。そうすると、気持ちがよくなるし、自信も帰ってくるし、自分を再び認識する。
　自分の問題に取り付かれている際、他の人も対称的に何かに直面しているのを忘れてしまう。他人の長所を羨むとき、他人も同時に自分の何かを羨んでいるだろう。自分が持っている長所に慣れたので、あまり気付かないけれど、他の人に取って達成できない特徴なのかもしれない。人はそういう風に他人の長所を大げさに見て、自分のを無視してしまい、逆に自分の短所を無視して、他人のを大げさに見てしまう。自分自身を向上させるためにはそもそもよいことなのにもかかわらず、自分の存在価値がないと思ってしまうことになる。とにかく、挫折に対して肯定的な態度を取らなければならない。また、自分だけがいわゆる「不幸」を与えられるという考えを持ってはいけない。
　何をしてもしなくても時間が過ぎるんだ。ある人は時間を解毒剤として使う。確かにそれぞれの事情や心の傷などは流れ時間とともに治してくるけれど、ある事情は逆に悪化するし、手がつけられないほどになってしまう。水野敬也の「夢をかなえるゾウ」と言う本によると、誰にでも一日二十四時間を与えられているが、其れの使い方によって人生に影響する。この二十四時間は深い自責の念に苦しみ、他人を責め、神様に文句を言うためなのか？もちろんそうではない。こういう二十四時間はまさに挫折から生まれた課題を発見するために、大地を踏みしめてやり直すために、生活の速いリズムに追いつくために与えられたんだ。そうしないであくまで悲しんだら、利益のみならずもっと大きいストレスに苦しむようになってしまう。
　仲間のない方の人生の旅路より仲間のある方の旅路の方が楽だ。なぜなら、人生を悲観している際、もし誰かが「自分も色々な原因で困ってる」と教えてくれれば、自分は独りではないのに気が付く。友達とは互いに励ましあうしたり、サポートしたりして、人生は楽になる。もし何だかんだで悩んでいたら、絶対に思いを伝えず心に抑え込んではいけない。憂慮を抑え込んだら、時間が長ければ長いほど、事情が悪くなり、悲観から逃げにくくなってしまう。こういうときは気持ちを友達に伝えたら、状況がよくなってストレスも軽減される。仲間と共に歩けば、旅路はなとかスムーズになる。
　
　人生の不如意に対しては、工夫できる点が色々あるけれど、限界を越えている場合もある。例えば、仲間の質と言うことこそがすべて自分の手につけられるわけがないだろうか。確か自分が友達を選択できるが、いい友達がいるのかいないのか、ほとんど運だろう。大体結構いい友達をできるけれど、足りるかどうか分からない。どんな個性がよいか、状況によって違っている。あるとき、ただの抱擁をしてあげて笑わせたらいいけれど、またあるときは人の心の奥底に行って悩みを解き放させてあげないといけない。これが友人として一番難しい課題である。人の必要に応じて自分の快適帯から出なければいけないかもしれない。
　ある人は朗らかで陽気なので、いつも他人を笑わせられる。殆どの人にとって、こんな友達がいれば十分だ。でも、別の場合でこんな友達は「数秒間の笑い」しか与えられない。なぜなら、不安感がそんなに強いから、笑いの力は一瞬の効果しかなくて、すぐ無くしてしまう。こんな人にとって、不真面目な冗談をやってあげても役を立てない。好意からの行為だと分かるが、状況が全然改善していないので、どうしようもなく思って「アホなことしか言わないバカな友達だけいるのか」と思ってしまう。この場合だったら、是非勇気を集めて、真剣に会話をしてあげなければならない。
　多くの人が必要となるとき快く他人を助けるんだ。誰かが助けに来たら、是非手伝うけれど、もし聞かなければ、何の行動も起こさない。人を思いやるのに、積極的に思いやるほどではない、というわけだ。皮肉にも、多くの人は率先して助けに来ない。なぜなら、自分の問題を他人に教えることは大きな勇気が必要だ。社会にはこの原因で多くの人が鬱病を患うだろう。もしかして、皆が人々に思いやりをもっと持って接したり、自発に人々を見舞いに行ったりすれば、鬱病を患っている人の人数は何とか減るだろう。懸念の持つ人はもっと早く発見し、自分が救けにきたくない人にも役を立つし、この取り組みにおいて「地球上にはまだいい人がおる」と思わせて、何とか慰安になる。この点で人の友達として一番難しい使命を担われさせる。なぜ難しいと聞くなら、人は大体自分勝手なので自分の事柄を優先にするべきだ。他人を慰安するのを先取りするのは自分の事柄を延期しないと無理だ。

　今時の社会人は殆どの責任を「関係者」に帰してしまう。この課題をしなければならないとか、こちは工夫をしなければならないとか、とにかく自分と全然関係ないと考えるふしがある。自分が社会の一部分としていいサポーターにならなければならないという責任を全然無視してしないんだ。だからこそ今時の社会はまだいろいろな悲劇が起こしている。こういう問題は見た目に個人の事柄なので自分で頑張って解決することであるが、むしろ友達の方から工夫しなければ問題は消滅しない。皆、身勝手さを落とすように、世界に一層の注意を傾けるように頑張ろう。

久しぶり？

2009-06-08T07:30:00.000-07:00

もう物凄く長い時間　アップデートしてなかった…
私を思ってる？

最近しんどいな

発表二つあって、作文も書かなきゃ
CUDA期末プロジェクトもあった…

雪うさぎ〜

雪アルマジロ〜

メッシュと粒子のコリジョン検出を問題としてやってみた

見た目が簡単だけど、メッシュの面の数が多ければ多いほど、プログラムが遅くなっちゃう

よほど巧妙に実装しないと　ゴミになっちゃう。

とにかく、面のソーティングを必要とする。

今日、全然寝てなくて、それを実装した。
で、今日の研究の発表もぼやぼやで準備した…　大丈夫かな
3D剛体変形　やっとできた

　

次は…　もっと怖い数学だ…

サボる日

2009-05-29T23:45:00.001-07:00

今週の水曜日、本番の「サボる日」だった
本当におもしろかった

写真がまだ手に入れてない＝お＝

2009-05-22T04:09:00.000-07:00

ジョギング

2009-05-21T10:52:00.000-07:00

クラス登録

2009-05-20T20:44:00.000-07:00

CUDAでスペクトログラムを作った。
CPUでやったDSPのより速かった。

おもしろい音声を分析した…　6~10khz フィルターにかけた

何と、顔を入ってた…　音楽のプロヂューサー　余裕だな

来学期のクラスもそろそろ選ばなければあかん…

なんに取ればええ…

EE 51
EC 121a
E 10
L 130
PE 84
CS 90
(CS 156)
CS 171 TA
CS 11 TA

炎の泉

2009-05-16T11:14:00.000-07:00

Particle system + billboard + fire texture =

スイカ？

2009-05-16T00:41:00.000-07:00

夏に食べる果物といえば…　西瓜！
昨日、西瓜真っ二つに切って半分をそのままで食べた。
今日、創造性アップ
　ミルクを加えて、西瓜ミルク飲み物

結構うまかった！

Ditch Day いよいよくるぞ
今朝は　模擬Ditch Dayがあって、「ドラえもん」テーマのStackに参加した

121/130

2009-05-12T22:30:00.000-07:00

中間試験の添削がやっと終わった。
何人が締切りを過ぎて試験出してないのを除いて
　121/130 全部添削した…　
後は…　4週間の宿題と　期末試験を残っている。

Cudaで3Dフラクタルを実装した
まずQuaternionで4Dフラクタルを作って（2Dのと同じように）
一つ3Dの切断面を表現する…

フロヨー　とちょっと似るから
宿題のウェブサイトを基で、パロディーウェブサイト　を作った

87/130

2009-05-12T01:30:00.000-07:00

また　得点

87分の中間試験の高さ　それだけかい？？

cs11 完成！

2009-05-11T16:24:00.000-07:00

Haskell 三角形ゲーム　宿題　をやっと解けて完成した！
最後の宿題なので　この授業が正式に終わった　＾＾

複雑はずだが、Haskellで　僅かな行数　びっくり

形態定義

type Move = (Int, Int, Int)
type Board = [Int]
type Solution = (Board, [(Move, Board)])

合法な動きかどうかの判断

valid_move :: Board -> Move -> Bool
valid_move b (f, p, t) = (elem f b) && (elem p b) && not (elem t b

動く

make_move :: Board -> Move -> Board
make_move b (f, p, t) = fh ++ [t] ++ sh -- add the 'target' peg
　　where
　　　　fh = [x | x <- b, x < t, x /= f && x /= p]
　　　　sh = [x | x <- b, x > t, x /= f && x /= p]

全部可能な動きを取る

make_moves :: Board -> [(Move, Board)]
make_moves b = let valid_mv = (filter $ valid_move b) moves
　　　　　　　in
　　　　　　　　　do mv <- valid_mv
　　　　　　　　　return (mv, (make_move b mv))

完成したかどうかの判断

is_solution :: Board -> Bool
is_solution b = length b == 1

主なロジック　Board_i のソリューション　は反復で
Board_i から Board_{i+1} の動き　次に　Board_{i+1}のソリューション

update :: Board -> [[(Move, Board)]]
update b | is_solution b = [[]]
　　　　　| otherwise =
　　　　　　　let step = make_moves b
　　　　　　　in
　　　　　　　　　do x <- step
　　　　　　　　　　(map (\y -> x : y) $ update (snd x))

図で説明すれば
　もし Board_{i+1}のソリューション　もう知ってたり
　　　　Board_{i+1}^1 ------ Path 1 ----> Goal
　　　　...
　　　　Board_{i+1}^k ------ Path k ----> Goal
Board_i　から　Move 1 ... Move k を取れば　
　　　　　　Board_{i+1}^1 ... Board_{i+1}^k
になることも知ってたら
　Board_i のソリューション　は
　　Move 1 + Path 1 :
　　　　Board_i -- Move 1 --> Board_{i+1}^1 ---Path 1---> Goal
　　…
　　Move k + Path k :
　　　　Board_i -- Move k --> Board_{i+1}^k ---Path k---> Goal

とにかく　cs11 　完成！

でも　ec 11採点　まだ終わってない
61/130　だ
今日　85/130　が目標だ

チェンバー

2009-05-11T02:11:00.000-07:00

コンサートがまずまずで終わった

Smetanaの一曲よー出演したと思ったけど
　Reichaで　前に出たことのないミスが浮いてきた気がした？
　　まぁぁ　聴衆が分かんないかもしれない　小っちゃいなミスだから
　
　満足した＾＾

クインテット〜〜〜　→ジョンくん　→ジュンちゃん　→翔　→エリザベスさん　→ローラさん



ジョンくん　<3　　〜　ジュンちゃん　　<3

母と昼食をした　　「ハ宝氷」というものも食べちゃった

エラい疲れたから　寝てしまって　一日で別になんもしてなかった

チェンバーの練習もうないから　その時間をよく利用しよう…
遅らせたワーク　追って追って追って　　頑張りまーーーす
今週　いっぱいできたようにしよう！

緊張？

2009-05-10T02:05:00.000-07:00

やっと　45/130
まだ遠いけど
前作ったヤツ↓
得点
　　長い宿題をするとき、進捗の動きを見るため作ったんだ
　　　（達成感と自信を高める）

言うことで、今は　九点　！　まだまだまだ…………

明日、コンサートだ　すぐ練習しにいくつもりやで

すると　cs11 haskell lab5 やり続きたい
まぁ…　ある問題に　にっちもさっちもいかなかったが……

ホナ

採点採点採点

2009-05-09T17:34:00.000-07:00

130学生の中間試験を採点しなきゃT_T
130!!!何年も掛かるでしょ

進捗： 20/130
まだまだ　終点から遠いな　T_T

日曜日は母の日　全世界のハハ　おめでとうございます　エラいお疲れ様です
その日　コンサートがあって

2曲を弾く予定：
1. Smetana, String Quartet No.1 　楽章1 楽章3

2. Reicha, Quintet for Oboe and String Quartet in F, Op. 107
楽章1 Allegro non tanto
楽章2 Andante siciliano
楽章3 Minuet - Trio - Minuet
楽章4 Finale. Vivace

来年はオフ-キャンパス

2009-05-07T15:30:00.000-07:00

来年　240 Michiganに住むことになる

この感じだ



キャンパスから近いし、そばに公園もあるし
結構よい選択だと

最近付いてないな　
ナムさんがアイリッシュウィスキー買ってくれた＾ろ＾

PC復活！ピカ復活！！

2009-05-06T04:25:00.000-07:00

前　nvidia driver 何やかんやで死んじゃって
　宿題とか研究とか　やりたくてもやれなかったから
　　困ってた。　二三日の労力で　やっと直した　ホホホホ〜

でも　直すと　なんと　やる気がなくした　XD
言うことで　僕はやはり普通の学生じゃん

最近　子供のころずっとやってたゲーム
　を再発見しちゃった！　ピカチューバレーボールだった

まさか　ピカに毒されたんだ？

(じ)ぷー gpu!

2009-05-05T21:46:00.000-07:00

ラブでGPU宿題をやった　　CUDAで N-body simulation の実装だった

こういう感じ

結構かっこよかった

不幸続き

2009-05-05T05:31:00.000-07:00

なんてこったい！！！

gpuクラスは新しい領域に進んだ　……　　CUDA!!

でも　CUDAを使うため　nvidia driver を壊しちゃって
　　　　　　　ーー　openGL　さえ使えない　（なんのこと？？？？）
　　　　　　　　　　　＝　研究できない
　　　　　　　　　　　＝　なにもできない！！！
　　　ああああああああああああああああああああああ
ああああああああああああああああ
　　あああああああああああああああああああ

すぐ直したい…

もう　今日、卓球総当たり戦　パート２　だった
　　なんやかんやで　しゃがみ込みサーブに負けちゃった
　　　　ーーー　最初は０ー２で　落後…
　　　　ーーー　で　２ー２
　　　　ーーー　５ゲーム　１０ー４でリードしてたが、　自分もサーブミスだし、レシーブもミスだし
　　　　　　　　ジュース　二回で　負けちゃった
　　　　　　（●｀ε´●）ぶーーーーーーーーー

　　でも　左手結構進歩したと思う

後は、経済学の中間試験を取った。授業を取ってないが、
教授は「ソリューションを書いてください」と頼むと
　…　取らなきゃ　ということ

幸いだが、答えは全部正解だった…　

コンピューター直せぇぇぇぇぇ　　(｀□´) コラッ！

未来の…いや　今時の強者

2009-05-03T04:09:00.000-07:00

張怡寧 vs 石川佳純　＠準々決勝

残念ながら　結果的に　1-4で張に負けた　T_T
でも　１ゲームを女王の手から奪うこと　驚かせた
第一ゲーム　4-11で取られたが
第二ゲーム　11-7でいいプレイをやってくれた
　　　なんと　サーブで女王にミスを押し付けた
　　　　　↑↑↑　多分　高いトス　サーブを利用して　回転の変化幅が大きい　返しにくい
第三ゲーム　ターニングポイントだと思う
　　　　　一度　１対７で落後　その時から
　　　　　彼女の眼から　なんか　恐れが見えた
　　　　　　最後　4-11で取られた
第四ゲーム　また　4-11で
　　　　　これは　大きいな衝撃だった　取られた３ゲーム　全部 4-11で
　　　　　「限界か」と思わせたかも知れない
最後に　7-11で　負けた
　　　　…　勉強になっただろう
　　　　　　いい表現だと思う
　　　　　ー＞　レシーブミス　…　３つだけ（？）
　　　　　　　とにかく少なかった
　　　　　ー＞　いいサーブ　相手が混乱された
　　負けた理由は…　多分　張選手の激しさにおびえて　なんと自信がなくなって
　　　　　　　　　　　　　　　通常な自分なりのプレイを出せなかった
　　　　　　　　……　惜しかった

やりたくなる

2009-05-02T18:05:00.000-07:00

女子シングル
日本選手の佳純ちゃん　準々決勝へ！！

昨日の事だったが　　今日　放送を観た
感動！　次は　張怡寧と対戦！絶対刺激
張選手は「卓球女王」と呼ばれ　世界ランキング１だが
なんと　意外が発生…　かも？？

グーグルでこの写真を探した…　　…　超サイヤ人のよう！

世界卓球！

2009-05-02T05:55:00.000-07:00

寝坊しちゃった…
起きた時、もう2時…　　経済学の助教やってるけど、recitationという、学生を助ける時間…　それが2時だった。
まぁ…　教室に飛んで（間に合った？）　先週の宿題のし方法を教えてあげた。
今週＝中間試験　　皆　頑張れ！

…世界卓球の選手たちと一緒に！
ittf.comで
幾つのゲームを実況で観た

女子シングル　張怡寧　vs tóth krisztina
実力がちょっと遠かったけれど　いいゲームだと思った
張怡寧、めっちゃ強かった
　　　集中力　すごかった
　　　一貫なヒット

でも tóthとは　とても面白い選手である
なぜなら　珍しい男性っぽい卓球やってた　
----> 素晴らしいロブ
ボール感、いいプレーを見せてくれた　

男子シングル…
……台湾選手は全部敗れた…ちょっと早かったが…　みんな強すぎ　しょうがない？
大好きな　水谷隼選手も…　中国のChen Qiに　4対0で　KOされた…　残念
日本の松平健太　惜しかった！！！　３ー４で中国のMa Linに負けた　T_T
しかし　　吉田海偉　準々決勝に進出！　次のゲーム　世界一のWang Haoと勝負！
絶対に見る！

自分も卓球をやった　結構楽しかった
その前　バドミントン　Nam様とチェンさんとやった

今日、新しいProject Euler問題が出て　　やってみて　すぐ解けた！

この問題と似るから　困難無しで

最初は紙筆で解けたけど、Haskellでも実装した

１．素数の生成

primes :: [Integer]
primes = sieve [2..] -- manually exclude 1 as a prime
　　　where
　　　　　sieve [] = []
　　　　　sieve (x:xs) = x : sieve [y | y <- xs, y `mod` x /= 0]

で…　phi(N)/N の　最小値を探す

divv x y = fromInteger x / fromInteger y

loop :: [Integer] -> Integer -> Integer -> Int -> Double -> Int
loop (x:xs) n phi l c
　　| divv count (prod - 1) < c = l
　　| otherwise = loop xs prod count (l + 1) c
　　where prod = n * x
　　　　　　count = phi * (x - 1)

p243 :: Double -> Integer
p243 m = s * denom
　　where c = loop primes 1 1 0 m
　　　　　l = take c primes
　　　　　numer = product (map (\x -> x - 1) l)
　　　　　denom = product l

　　　　　looper k | r < m = k
　　　　　　　　　| otherwise = looper (k + 1)
　　　　　　　　　where r = divv (k * numer) (k * denom - 1)
　　　　　s = looper 2

getAnswer = p243 (15499 /94744)

これでもうすぐ半分の問題を解答したように

でも…残してるのは…　大体難しい…

今日から、生活は楽になる…　　と思う！
でも　予報によると　雨降りそう　T-T
まぁ　関係ない？　オタクである僕は部屋から出ないから

もうちょっと

2009-04-29T10:55:00.000-07:00

研究アドバイザー　ドイツから帰った
進行状況を調べることになる
そのため　寮の部屋に来るそうです

------------------------------
結果的に、Deformation数学導出の討論に打ち込んで、
自分の部屋にいかなかった。

Laplace Beltramiの討論だった。

次の内方法で、同じ結論になるかい？
１．　Discrete Exterior Calculus
-- △ = d★d
２．　direct algebraic derivation

....次の研究だ

CS11 Haskell の宿題をした。

あまり難しくないが、謎のような言語…　一部分は下記

main = do args <- getArgs
　　　　　　if length args == 0
　　　　　　　　　　then interact id
　　　　　　　　　　else mapM_ (\f -> (readFile f) >>= putStr) args

夜は「夢をかなえるゾウ」二章を読んだ

今回二つの課題は
１．腹八分にすませる
ー　特にこんな時間（~午前4時）

２．人が欲しがっているものを先取りする
ー　難しいなぁ〜でも、やってみる！

ー＞　宿題　＜ー

長い日

2009-04-28T03:38:00.000-07:00

まだ全然寝てない　もう36時間

今日　午餐会があって　大喰い　ラッキー！

後は　卓球クラスだ
教練が世界卓球大会のため　横浜に行った
なので　こっちは　総当たり戦をやった

その前　学校のラケット（シェーク）で　ウォーミングアップした
ペん（持ち方）だけど、意外にやり易かった！
重心が自分のラケットのより後ろ　かも…
なので、　ラケットの重心を調整した

やっぱりやり易くなった　＝DD

研究を遅らせて、代わりにHaskellでProject Eulerをやった。
今日はこの問題に挑戦！

基本的に、文字二つ｛O,A｝で構成されてAが三回連続で出てない文字列　幾つある　
　　という問題を手前に解決しなきゃ

で　Lが出ない場合と一つ出る場合　を別々で考える。

という訳で、Haskellで実装：

tribonacci :: [Integer]
tribonacci = 1 : 2 : 4 :
　　　　zipWith (+) (zipWith (+) tribonacci (tail tribonacci))
　　　　　　　　　　　(tail (tail tribonacci))

p191 :: Int -> Integer
p191 n = noL + oneL
　　　　where noL = last (take (n + 1) tribonacci)
　　　　　　　　l = take n tribonacci
　　　　　　　　oneL = sum (zipWith (*) l (reverse l))

wii project　ガンバレ！

いつもの

2009-04-27T06:09:00.000-07:00

またオールした

昨日（土曜日）　珍しいパーティーがあって
ちょっと　飲み過ぎちゃって　二日酔いだった
２１才になったから　初めのアルコール　（合法的に！）

アイリッシュコーヒー飲みたいな〜

Colors何回も聴いた　今までも聴いています

研究やってる！As-rigid-as-possible Deformation

使用者が幾つの頂点を固定して、頂点を移動して
最初のメッシュと一番似るようにするため　固定してない頂点の配置を計算する　
という研究だ。

先学期２Dでやってた。今、３Dに進む！
ネックになってたが…　なので　毎日オールする！

フォォォォォォ！

新しい出発

2009-04-25T15:01:00.000-07:00

今日から　なんと　日本語でブログを書くことになる
日本語がたどたどしいけど　もっと使えば上達する…　だろう？

まぁ　学業で忙しい！　いくつの授業を落としたが、まだ締切りに遅れた研究（プログラム実装）をしなきゃ
来週過ぎると　楽になる？　そう望んでるけど

今学期　haskellを習うことになった
めっちゃおもしろいプログラミング言語だと思う
今日は　初めて（haskellで)Project Eulerの一つの問題を解けた
http://projecteuler.net/index.php?section=problems&id=74
あまり難しくないけど（最近の問題と比べると）　haskellの勉強になった

GPU授業を取ってるが、取ったことがあるから　楽勝！
何といっても、興味があっておもしろい

最近の実装宿題

日本語授業のため　「夢をかなえるゾウ」という本を読みかけた
'Little Tokyo'にある本屋さんで思い切って買ったが　T_Tたかぁぁ

…というわけ　第一ポストです

研究をしなきゃT_T

翔〜 望♭