やまかわのログ http://antas.jp/blog/yamakawa/ ja Copyright 2010 Thu, 14 Oct 2010 16:43:21 +0900 http://www.sixapart.com/movabletype/?v=3.33-ja http://blogs.law.harvard.edu/tech/rss Foundation of bezier curve (3) 実は、これについても簡単に求める方法がなく、以下の式の根を求めることになります。 nearest.gif ここでB(t)は1回目で紹介したベジェ曲線の多項式です。上記の式は5次方程式になり、解の公式が存在しないため数値解析等を利用して解くことになります。 解き方としては以下のものが代表的ですが、他にも色々あるようです。
  • Newton法
  • Bezier Clipping
  • Bisectionning
参考： Accurate and Efficient Algorithm for the Closest Point on a Parametric Curve他 この論文ではBezier Clipping同様にベジェ曲線のConvex hull性を利用しています。]]> http://antas.jp/blog/yamakawa/2010/10/foundation_of_bezier_curve_3.html http://antas.jp/blog/yamakawa/2010/10/foundation_of_bezier_curve_3.html math Thu, 14 Oct 2010 16:43:21 +0900 Foundation of bezier curve (2) この式に出てきたベクトルを下記のように選んで新しいベジェ曲線を２つ作ります。すると、曲線の形を維持したままtの位置で分割された2つのベジェ曲線を得ることができます。 splitted.gif 参考： The Essentials of CAGD http://www.farinhansford.com/books/essentials-cagd/]]> http://antas.jp/blog/yamakawa/2010/10/foundation_of_bezier_curve_2.html http://antas.jp/blog/yamakawa/2010/10/foundation_of_bezier_curve_2.html math Thu, 14 Oct 2010 16:22:24 +0900 手稲高校インターンシップ 2010 (http://scratch.mit.edu/) を利用することにしました。 カリキュラムとしては物体の動かし方を学んでもらってから、最終的にpongっぽいものをつくるというかんじでやってみました。 ことしのみなさんはLOGOやC/C++、なでしこ（！）の経験者という強者ぞろいだったため、事前に用意したネタが尽きてしまうというトラブルもありましたが、おおむね喜んで頂けたようです。 当日作ったスケッチを以下に貼ります。実行にはブラウザへのSun Java プラグインインストールが必要です。 Scratch Project (英語が間違ってるのに気付きましたが直してもサムネイルが更新されません！ぎゃー）]]> http://antas.jp/blog/yamakawa/2010/08/_2010.html http://antas.jp/blog/yamakawa/2010/08/_2010.html Tue, 24 Aug 2010 12:21:52 +0900 Foundation of bezier curve 三次ベジェ曲線は次の式によって与えられます。 ここでtは長さを表すパラメータで0≦t≦1です。P0（始点）が0、P4（終点）が1になります。P1,P2は制御点です。 ベジェ曲線はふつうの直線の描画とことなり、媒介変数tを動かしてそのtの時のx,y座標を求める、というやり方でしか描画ができなかったりします。 さて、このパスを適当なtの位置（長さ）を決めて、パスの形をそのままにそこで切断するとした時、どのようにしたらよいでしょうか。 これには、de Casteljau's algorithmというのを使います。de Casteljauはベジェ曲線を考えた人ですが、これを使うと新しいコントロールポイントの位置が多項式を逐次的に計算していくだけで簡単に決まります。(続く) 参考： The Essentials of CAGD http://www.farinhansford.com/books/essentials-cagd/]]> http://antas.jp/blog/yamakawa/2010/03/foundation_of_bezier_curve.html http://antas.jp/blog/yamakawa/2010/03/foundation_of_bezier_curve.html math Thu, 25 Mar 2010 12:43:00 +0900 NSImageに文字を描画するとさかさまになる？ cocoa-devで同じようなハマりをしている人が見つかりました。 focusが当たっているビューの座標系が上下反転させた状態で、 NSImageにlockFocusでfocusを当てて画像上に文字を書くと、この症状が起きるようです。 Cocoaの仕様的にどうにもならないようでしたので、 NSImageにsetFlipped:YESを設定した状態で反転したビットマップを作り、 さらにそのあとにNSAffineTransformで上下反転マトクリクスを作成して 現在のグラフィクスコンテキストに結合し、 上下反転コピーを行ったビットマップをもう一個作ることで なんとなく対処できました。 一応予定の見栄えにはなったのですが、 2回ビットマップができるので大変無駄臭いような気がします…。]]> http://antas.jp/blog/yamakawa/2010/01/nsimage.html http://antas.jp/blog/yamakawa/2010/01/nsimage.html OSX Wed, 06 Jan 2010 18:34:02 +0900 scipyの統計関数が爆速なことについて JNI-BLASバインディングを使わなくてもよさそうです。そのうち比較してみますが。 Javaの場合演算子オーバーライドがないですので微妙にまどろっこしい印象はいなめませんがしかたないですね。]]> http://antas.jp/blog/yamakawa/2009/11/scipy.html http://antas.jp/blog/yamakawa/2009/11/scipy.html Wed, 25 Nov 2009 17:22:42 +0900 PythonでHTMLのスクレイピング：scrapemark http://arshaw.com/scrapemark/docs/examples/というのを使うと大変簡単にできました。 使い方はまあ簡単で
>>> scrape(""" <div id='content'> {{ before_table }} <table /> </div> """, url="http://example.com") >>> {'before_table': 'Look at these data points'}
みたいなかんじです。 HTMLの繰り返しになっているところの中身をリスト化してとってくることもできます。これは
>>> scrape(""" <body> {* <div class='section' id='{{ [section_ids] }}' /> *} </body> """, html) >>> {'section_ids': ['nav', 'content', 'footer']}
というかんじで大変使い勝手が良いですねー。開発元のexampleを見ると、もっと構造化したデータも取ってこれるようです。 RubyでNokogiriとかするより楽な気が致しますねえ。 あと、スーパー便利な事に結果にはユニコード化された文字列が入ってきますので扱いも楽ちんです。]]> http://antas.jp/blog/yamakawa/2009/10/pythonhtmlscrapemark.html http://antas.jp/blog/yamakawa/2009/10/pythonhtmlscrapemark.html Wed, 21 Oct 2009 17:20:04 +0900 SQLAutocode http://turbogears.org/2.1/docs/main/Utilities/sqlautocode.html 使い方は
easy_install sqlautocode sqlautocode -d -o model.py sqlite:///moviedemo.db
とかするだけです。 べりー簡単ですねえ。 追記：relationはどうにもならないので手で書いて上げる必要があります。]]> http://antas.jp/blog/yamakawa/2009/10/sqlautocode.html http://antas.jp/blog/yamakawa/2009/10/sqlautocode.html Wed, 21 Oct 2009 17:10:07 +0900 Ruby + GDB http://github.com/michaelklishin/gdb-macros-for-ruby などから取得してきたスクリプトを~/.gdbinitに追加しておけば、あとはデッドロックしているプロセスにgdb -p pidでアタッチして、GDB内でrb_backtraceなりrb_raiseなりを実行するとバックトレースが取れます。Googleで検索すると、他にもいくつかスクリプトがあるようです。 Pythonでも同じようにGDB用のスクリプトが用意されているようです。 http://wiki.python.org/moin/DebuggingWithGdb 不意にRubyがC Extention内などでSEGVってしまった場合などでも同様にインスぺクションが可能ですので、なかなか便利なように思います。 追記：コアを取るようにしていた場合はもちろんコアダンプの解析にも使えますね。]]> http://antas.jp/blog/yamakawa/2009/10/ruby_gdb.html http://antas.jp/blog/yamakawa/2009/10/ruby_gdb.html Tue, 06 Oct 2009 13:41:14 +0900 インターンシップ受け入れ ball.png http://www.openprocessing.org/visuals/?visualID=3128 時間が足りなかったため、 ↑の打ち出し機は実際には固定になってしまいました。 来年があったらハードウェアにしますかねー。 AVRとか。 （こっそり最近ICE買ったので。。） ]]> http://antas.jp/blog/yamakawa/2009/07/post_1.html http://antas.jp/blog/yamakawa/2009/07/post_1.html Fri, 31 Jul 2009 19:57:58 +0900 Velocity 2009(2) Velocity 2009の続きです。 Building OpenDNS Stats OpenDNSというDNSを外部提供してる会社の人のログ捌きの話です。ログが一日140億件くらいあるそうです…。生ログは阿呆ということでとりあえずグラフ、というのと、ログを一個のテーブルに入れるのは問題外なのでネットワーク毎に別テーブルにしてる、オートインクリメントはロックがかかるのでタコだからSHA1でPKを作ってる、iowaitの80%はファイルのopen/closeなんでmysqld_safe に ulimit -n 600000をセットしとく、ログ集計処理は専用のサーバを3台用意してバッチ処理してる、集計クエリはLIMITつけないで発行してmemcachedに丸ごと入れておき、ページングは手元でやるなどなど。 Death of a Web Server: Crisis in Caching あんまり中身がなかったので省略。珍しくWindows ServerベースのWebサーバの話だったのでなんか特有のスケーリングてくにっくとかあるのかと思ったらそういう話も特になく…。 Fixing Twitter: Improving the Performance and Scalability of the World's Most Popular Micro-blogging Site えーまず2008年に752%（げえ）ユニークビジターが増えたという話から始まります。ていうかグラフの角度が怖いです。このくらいになるとちゃんと計算機科学を使わんと無理、とか消火活動するより予測を行って事前にCapacity Planningしておきましょう、とか。Webサーバのログを拾って503エラーとか500エラーの件数はグラフ作ってるそうです。これはいいですねー。サーバの設定変更はレビューを必須にしてて、履歴をレビューツールに記録してて、Subversionはchangelogに"reviewed"っていうのが入ってないとはじくようにしてる、等々。 プログラム側の方の話ではなんでも計測できるようにしておこう、とかメトリクスは推測じゃなくて決断のために使いましょう、あとキュー使おうなど。色々書いてあって忘れてきました。ActiveRecordのキャッシュライブラリCache Moneyというのが割と気になりました。 キャッシュの話。全体設定の変更とかでキャッシュが一気に整合性壊れるのとかどうすんべ、memecachedもたまにSEGVるなど。]]> http://antas.jp/blog/yamakawa/2009/07/velocity_20092.html http://antas.jp/blog/yamakawa/2009/07/velocity_20092.html Wed, 15 Jul 2009 11:50:24 +0900 Velocity 2009 Asian network内ので配信されているBhangra番組のクオリティの高さにびっくりしたりしていましたが、ぼちぼちVelocity2009の資料が公開されてきているので、上から順番に眺めてみた感想をつらつら書いていこうと思います。 10+ Deploys Per Day: Dev and Ops Cooperation at Flickr The Art of Capacity PlanningでおなじみFlickrのマネージャJohn Allspawさんらによる、運用チームと開発チームの関係を如何に円滑にマネージして行くかという話です。兎角運用は「オレのサーバは悪くない！おまえのヘボコードのせいだ！」となり一方開発は「オレのコードは悪くない！おまえのヘボサーバのせいだ！」となりがちですが、そこらへんの責任の投げ合いをどうやって回避していくかというような大変心に響く内容でした。開発の方も運用の方も一読をお勧めいたします。 2 Years Later, Loving and Hating the Cloud 一年半くらいオンラインレタッチつき写真カタログサービスをAmazonのCloud関係で運用してるPicnikの方のAmazonCloud感想です。正直あんまり印象が残らず…。「 S3に要らんデータを残しとくと無駄にお金がかかるけど、あんまりうまい削除の方法が思いつかんので放置中」という話くらいでしょうか。後はサーバリソースの使用状況に大きいスパイクがあるような運用ケースでは結構使い出があるというグラフなど。 Ajax Performance JavaScript: The Good Partsという著書（この本は知らなかったので後で調べます。。）のあるDouglas CrockfordさんのAjaxでどうやって性能を出していくかという話です。ブラウザのヘッポコ加減やDOM操作の遅さに嫌気がさしている方は結構いらっしゃるかと思いますが、DOMにノードを追加する前にプロパティの設定は終わらせておくなど、「あるある」ネタが割と書いてあって良かったと思います。 Beyond Gzipping HTTPストリームをGzipというかdeflate対応にしておくと転送量が減るのはすでにおなじみですが、なんかgzipで転送されないユーザーエージェントが結構あるなー(15%くらい）と思ったら、パーソナルファイアウォールとかのプロクシ機能を持つセキュリティソフトが、HTTP通信の内容をgzipされると中身を読めないので勝手にHTTPヘッダの「Accept-Encoding:gzip, deflate」を削除しくさってた、というかなり「まじでー」という話でした。 続きは後日！]]> http://antas.jp/blog/yamakawa/2009/07/velocity_2009.html http://antas.jp/blog/yamakawa/2009/07/velocity_2009.html Fri, 03 Jul 2009 14:07:31 +0900 Rackデーモンの起動スクリプトを書いてみる 割合最近マルチスレッドかつそこそこ高速に、簡単なHTTPリクエストをDBを引いて捌くサーバプログラムを書く機会がありました。

そういうのにMerbをとかを使うのもヘビーウェイトだなあと思い、Ramazeとかもいんじゃないかなあと少し思いましたが折角なのでRack+Sequelでいってみることにしました。 RackについてはMerbでおおよそ知っていましたのでさっくり作ることができたのですが、デーモンの起動スクリプトを作るところで少し困ってしまいました。

Rackは出来合いのrackupというコマンドがあり、これでサーバを起動します。-dオプションを付けるとデーモンとしてバックグラウンド起動するのですが、たいへん不便なことに停止するコマンドがありません。なのでとめる時はpidをもとにSIGINTをおくってやらないといけないです。

また、rackupには実行ユーザーを指定する方法がありません。したがって、rackupコマンドを実行する際にsu -cを利用して、実際の実行ユーザーでデーモンを起動する必要があります。

これらにより以下のようなスクリプトを書きました。出来あがったプログラムはRubyとは思えない速さでウルトラ高速に動作して大変よかったですので機会を見つけてまた使ってみようと思います。

]]> http://antas.jp/blog/yamakawa/2009/05/rack.html http://antas.jp/blog/yamakawa/2009/05/rack.html Wed, 13 May 2009 18:46:51 +0900 mod_xsendfile+mod_proxyの小ネタ mod_xsendfileが利用できる場合はmod_xsendfileを使って、 そうでない場合は自前で送るようにしたいなあ、というときの小ネタです。

そういうことをするには、mod_xsendfileが有効であるかどうかをなんとか 判別しないといけないわけですが、そのためにmod_headerを使って、 mod_xsendfileが有効な時だけ付加的なヘッダをねじこんでやるとうまいことが判別ができました。

]]> http://antas.jp/blog/yamakawa/2009/04/mod_xsendfilemod_proxy.html http://antas.jp/blog/yamakawa/2009/04/mod_xsendfilemod_proxy.html 小ネタ Wed, 15 Apr 2009 13:24:48 +0900 RubyでUDPブロードキャストを送受信してみる その２ 同一マシン内外の複数プロセスにブロードキャストメッセージを送る必要があり、 さて不特定多数に送るならUDPブロードキャストかなあと思ったのですが、 同じポートにbindできないといけないからダメだなーとと思ったところ、 実はUDPだとsetsockoptでSO_REUSEADDRをセットすると 複数プロセスで同一ポートを待ち受けるようになり、 マルチキャストやブロードキャストもそれぞれのbindしているプロセスに届けられるようです。

ということを知ったので以前のRubyでUDPブロードキャストを送受信してみるを改修してみました。ついでにアドレスなどをSocketクラスの定数を使うように置き換えています。

]]> http://antas.jp/blog/yamakawa/2009/03/rubyudp_1.html http://antas.jp/blog/yamakawa/2009/03/rubyudp_1.html Wed, 25 Mar 2009 15:23:44 +0900