Don't panic !

USB to RS485/RS422 with protection circuit and Isolation housing

2012-02-15T09:40:00.000-08:00

http://www.kisansystem.kr/index/shop/item.php?it_id=1263307062
http://www.eleparts.co.kr/front/productdetail.php?productcode=035011001000000074

합리주의 혹은 중용의 도

2012-02-15T08:14:00.000-08:00

중용과 합리주의 이 둘의 언급 순서를 어떻게 해야 할지 먼저 생각하게 된다. 이 둘은 울림은 개인적으로 비슷하게 느껴지지만, 아직은 중용 보다는 합리주의라는 말이 가슴에 와 닿는 말이다. 아무튼 같다고는 할 수 없지만 이 둘은 공통된 부분을 공감하며 살아간다는 일들이 왠만히 부지런하지 않으면 쉬운 일이 아닌 것 같다. 당연히 지켜져야 할 일들이 짜증나게 사소한 것들에서 부터 추상적인 담론까지 챙기지 못한다면 못된 사람들의 꽃노름에 휘둘이기에 쉽다. 너무나도 혼잡하다. 개인적인 내적 탐구도 힘든데 외부적으로 바른 길을 찾는 길도 힘들다.

AVR development environment with eclipse in debian

2012-02-13T23:44:00.000-08:00

* install "eclipse"
# apt-get install eclipse-cdt

* install "AVR-plugin"
Execute "eclipse" -> Help -> Install New Software... -> Add -> type "AVR Plugin" in Name tab and insert "http://avr-eclipse.sourceforge.net/updatesite/" in Location tab -> install "AVR Eclipse Plugin"

*install "C99 LC Parser" in eclipse

* install gcc-avr, task-c-devel, gcc-doc, gcc-4.2 avr-libc, binutils-avr, and avrdude, avrdude-doc libftdi1
# apt-get install gcc-avr task-c-devel, gcc-doc, gcc-4.2, avr-libc, binutils-avr avrdude avrdude-doc libftdi1

미국 일하러 가기 준비

2012-02-11T07:16:00.000-08:00

* 보통 연봉 : 5만달러 (현재 시점 6천만원)

* H1비자는 세금 25%

* 9시 출근 5시 퇴근

* Silver Spring에 위치한 Lado라는 학교와 University of Maryland 의 어학 연수프로그램, 전자는 회화위주의 사설 교육기간 Beltsville에서 차로 15분 정도거리 후자는 차로 13분 정도의 거리이다.

* 가정용 공구 : 몽키 스패, 펜치, 프라이어, 니퍼, 롱노즈프라이어, 망치, 드라이버 +,- x 대소너, 육각렌치

* 전기 장판 : 110, 220 겸용 http://itempage3.auction.co.kr/DetailView.aspx?ItemNo=A594372381&frm3=V2 한일 전기요 대형 180x 135

* 주거지
rollins park
http://www.rollinspark.net/floorPlans.php
congressional tower
http://www.congressionaltowers.net/
mira mont
http://www.southernmanagement.com/miramont/
https://www.southernmanagement.com/miramont/apply/

* 이민가방
http://www.brightonmall.co.kr/mall/index.php

* 인터넷 전화기 준비 (LG U+ 070, 기본 요금: 2200원)

[번역] Chrome 브라우져의 cache 영역을 ram으로

2012-02-11T07:11:00.000-08:00

다음 문서는 Chrome의 cache영역을 ram disk(/dev/shm)밑으로 설정하고 전원이 꺼질때 유실되는 데이터를 하드 디스크로 복사하고 부팅시 복사된 하드 디스크의 데이터를 다시 ramdisk 밑으로 복사하는 과정을 거치는 노력을 하였다.

출처 : http://www.anders.bennehag.com/blog/2010/putting-chromes-cache-in-the-ram/
Chrome 웹 브라우져는 속도를 향상하는 방향으로 변화하고 있다. 그러나 지속적인 디스크 읽기와 쓰기로 내 Laptop 컴퓨터의 배터리와 하드드라이브를 소모하게된다.

다행이도 리눅스는 Ram으로 Mount 폴더를 가지고 있다. 이건 읽고 쓰는 작업이 하드 디스크에 영향을 미치지 않는다는 것을 의미한다. (실제로 내 Laptop의 디스크보다 2배정도 빠르다) 단점은 전원을 끄게되면 데이터가 사라지게 된다. 이 폴더는 통상 /dev/shm으로 불린다. 나는 이 폴더를 Chrome 브라우져의 cache로 사용하고자 한다. 그러나 나는 컴퓨터의 전원을 끌때 환경 설정, 확장 프로그램, 그리고 저장된 데이터의 잃어 버릴 준비가 되어 있지 않다. 그래서 나는 Ram으로 mount된 /dev/shm에 있는 폴더를 디스크로 복사하는 스크립트를 작성하였다. 그리고 이 스크립트는 Crome이 종료 될 때마다 실행된다.

이것은 Chrome의 cache를 Ram으로 설정하는 것이다.

# mkdir /dev/shm/.personal_synced # The folder in ram
# mkdir ~/.personal_synced # The mirrored folder on disk

매 시간마다 실행되는 간단한 스크립트를 작성하였다. (나는 좀더 쉬운 방식으로 작성되어 질 수 있다고 생각한다. 하지만 나는 아직 bash에 초보자 이다.)

#!/bin/bash
# Sync the ramdisk to harddrive

foldername=”.personal_synced”
ramdisk=”/dev/shm/”$foldername
harddisk_copy=”/home/anders/”
backup=”/home/anders/”$foldername”.old”

# Remove the old backup
echo “removing $backup”
rm -fr “$backup”
# Make a copy of the old one
echo “Copies the latest backup $harddisk_copy$foldername”
mv -f “$harddisk_copy$foldername” “$backup”
# Copy from ram
echo “Copies $ramdisk to harddrive, $harddisk_copy”
cp -R “$ramdisk” “$harddisk_copy”

나는 이 sync 스크립트가 Chrome이 종료 될때 실행되기를 원한다. 먼저 나는 나는 전원이 종료될때 실행되도록 시도하였다. 다음 스크립트는 Chrome이 종료될때 Ram 디스크를 동기화하는 것이고 난 그져 이 스크립트를 가지고 /usr/bin/google-chrome에 있는 symbolic link를 대체하였다.

#!/bin/bash
# Start Chrome with cache in the ram

/opt/google/chrome/google-chrome -user-data-dir=”/dev/shm/.personal_synced/google-chrome” “$@” && /home/anders/Scripts/syncRamdisk

소실된 모든 것들은 컴퓨터가 부팅될때 디스크 폴더의 것들을 복사하는 스크립트가 보상한다.

# The synced folder
cp -R /home/anders/.personal_synced /dev/shm/
chown -R anders /dev/shm/.personal_synced

# The chrome cache. I don’t want to sync this so it’s placed outside .personal_synced
mkdir /dev/shm/google-chrome
chown anders /dev/shm/google-chrome

그리고 단순히 시작 프로그램에 이 스크립트를 추가하였다. System->Preferences에 해당하는 곳에 하면 된다.

Chrome의 cache를 우리의 자동 synced 폴더로 이동하자.

.config/google-chrome-folder를 synced 폴더로 이동하고 symbolic link로 대체함으로써 작업은 완벽해 졌다. (나는 또한 cache를 ramdisk로 옮겼으나 synced 폴더를 그렇게 하지 않았다. 나는 cache를 유지하는데 필요한 공간과 소비 시간에 아무런 변화점을 보지 못했다.)

# mv ~/.cache/google-chrome /dev/shm/google-chrome
# ln -s /dev/shm/google-chrome ~/.cache/google-chrome

# mv ~/.config/google-chrome /dev/shm/.personal_synced/
# ln -s /dev/shm/.personal_synced/google-chrome ~/.config/google-chrome

주의 사항들:

/dev/shm/는 모두에 의해서 읽기가 가능하다. 이 폴더에 어떤 것을 넣는 것은 매우 안전한 작업은 아니다.

추가작업: 나는 스크립트를 약간 수정하였다. root 소유의 폴더에서 문제점이 발생하였다.

/dev/shm 무엇에 쓰는 물건인가?

2012-02-08T07:10:00.000-08:00

POSIX 기반의 공유 메모리를 위한 파일 시스템으로 응용프로그램들이 데이터 공유를 위해서 잠시 사용하기는 용도로 규정되어 있는 것으로 보인다. 문제는 이 파일 시스템이 메모리를 저장 공간으로 사용하기 때문에 저장및 읽기 속도는 빠르지만 hibernation상태로 진입하게 될때 처럼 전원이 차된되는 경우 데이터를 유실할 가능성이 있다. 그래서 swap영역을 사용하도록 규정되어 있다고 언급된다. 그레서 shm의 데이터는 hibernation상태로 진입할때 하드 디스크영역에 해당하는 swap영역에 데이터를 저장한다고 여겨진다.

RAM upgrade 후 (2g -> 4g) tmpfs의 용량 바꾸기

2012-02-07T20:06:00.000-08:00

tmpfs로 mount된 /dev/shm의 용량은 ram의 용량이 2g였을 때의 설정을 그대로 가지고 있다.

# df -h /dev/shm

Filesystem Size Used Avail Use% Mounted on
tmpfs 760M 1.1M 759M 1% /run/shm

음. ram이 2g 였었지만 비디오 카드가 사용하던 용량이 300M정도 였었음을 감안하면 얼추 맞는것 같다. 90M가 모자르지만 맞다고 하자.

그럼 용량을 바꾸어보자.

# vi /etc/fstab

어라 /dev/shm의 mount 정보가 없다.

좀더 찾아 보니 초기화 스크립트로 옮겨 갔다.

# ls /etc/rcS.d

README S09checkfs.sh S17nfs-common

S01mountkernfs.sh S10mountall.sh S18mountnfs.sh

S02udev S11mountall-bootclean.sh S19mountnfs-bootclean.sh

S03mountdevsubfs.sh S12mountoverflowtmp S20console-screen.sh

S04bootlogd S13ifupdown-clean S21kbd

S05keyboard-setup S13pppd-dns S22console-setup

S06hdparm S13procps S23alsa-utils

S06hostname.sh S13udev-mtab S23bootmisc.sh

S06hwclockfirst.sh S13urandom S23fuse

S07checkroot.sh S13x11-common S23lm-sensors

S08hwclock.sh S14ifupdown S24stop-bootlogd-single

S08module-init-tools S15networking

S08mtab.sh S16rpcbind

여기서 S03mountdevsubfs.sh에 해당 내용이 들어 있다.

# Mount a tmpfs on /run/shm

if [ ! -d /run/shm ]

then

mkdir --mode=755 /run/shm

[ -x /sbin/restorecon ] && /sbin/restorecon /run/shm

if [ yes = "$RAMSHM" ] || read_fstab_entry /run/shm tmpfs; then

domount "$MNTMODE" tmpfs shmfs /run/shm tmpfs "-onosuid,nodev$SHM_OPT"

else

chmod "$SHM_MODE" /run/shm

# Migrate early, so /dev/shm is available from the start

if [ "$MNTMODE" = mount_noupdate ] || [ "$MNTMODE" = mount ]; then

run_migrate /dev/shm /run/shm ../run/shm

음.. 잘 모르겠다 용량 설정하는 내용이 없다. 아마도 SHM_OPT라는 변수에 있을 것 같다.

여기(/lib/init/tmpfs.sh)에 다있네

# Set tmpfs vars

#에

# Set defaults for /etc/default/tmpfs, in case any options are

# commented out which are needed for booting. So that an empty or

# outdated file missing newer options works correctly, set the default

# values here.

TMPFS_SIZE=20%

TMPFS_MODE=755

RUN_SIZE=10%

RUN_MODE=755

LOCK_SIZE=5242880 # 5MiB

LOCK_MODE=1777

SHM_SIZE=

SHM_MODE=1777

TMP_SIZE=

TMP_MODE=1777

# Source conffile

if [ -f /etc/default/tmpfs ]; then

. /etc/default/tmpfs

RUN_OPT=

[ "${RUN_SIZE:=$TMPFS_SIZE}" ] && RUN_OPT=",size=$RUN_SIZE"

[ "${RUN_MODE:=$TMPFS_MODE}" ] && RUN_OPT="$RUN_OPT,mode=$RUN_MODE"

LOCK_OPT=

[ "${LOCK_SIZE:=$TMPFS_SIZE}" ] && LOCK_OPT=",size=$LOCK_SIZE"

[ "${LOCK_MODE:=$TMPFS_MODE}" ] && LOCK_OPT="$LOCK_OPT,mode=$LOCK_MODE"

SHM_OPT=

[ "${SHM_SIZE:=$TMPFS_SIZE}" ] && SHM_OPT=",size=$SHM_SIZE"

[ "${SHM_MODE:=$TMPFS_MODE}" ] && SHM_OPT="$SHM_OPT,mode=$SHM_MODE"

TMP_OPT=

[ "${TMP_SIZE:=$TMPFS_SIZE}" ] && TMP_OPT=",size=$TMP_SIZE"

[ "${TMP_MODE:=$TMPFS_MODE}" ] && TMP_OPT="$TMP_OPT,mode=$TMP_MODE"

아 잘 모르겠다. 시간을 두고 살펴봐야지.

RAMDISK, ramfs, 그리고 tmpfs사이에 생기는 혼동!

2012-02-07T00:38:00.000-08:00

이후의 글은 다음 링크를 참조하였다.
http://www.linuxfordevices.com/c/a/Linux-For-Devices-Articles/Introducing-initramfs-a-new-model-for-initial-RAM-disks/
ramdisk와 ramfs의 차이

initrd같은 ramdisk는 ram 기반 block device이다. 이건 disk같이 mount되고 포멧된 고정된 크기의 메모리 공간을 의미한다. 이것은 ramdisk의 내용들은 일정한 형식을 갖는 형태로 되어 있어야 하고 특별한 도구(mke2fs, losetup)로 준비되어야한다. 그리고 모든 block device들과 같이 실행시 데이터를 해석하기 위해서 filesystem 드라이버가 필요하다. 이건 또한 낭비되는 공간(ramdisk가 꽉차지 않으면, 나머지메모리는 어떠한 것으로 사용될 수 없다.)이나 공간 제한(ramdisk가 꽉 채워졌지만 나머지 메모리 공간이 여유가 있더라도, 포멧없이 ramdisk의 확장은 불가능하다.)인위적인 크기 제한을 의미한다.

그러나 ramdisk들은 실제로 caching 때문에 많은 공간을 낭비한다. Linux는 block device를 쓰거나 읽기 위한 모든 파일과 디렉토리 entry를 cache로 처리하도록 고안되어져 있다. 그레서 Linux는 ramdisk에 쓰거나 읽게된 데이터를 "page cache"나 "dentry cache"로 복사한다. ramdisk의 아랫쪽 부분은 block device처럼 다루어 진다.

몇년 후 Linus Torvalds는 잘정리된 아이디어를 가졌다. 리눅스의 cache시스템이 filesystem처럼 mount되어 진다면 어떨까?하는... cache에 파일을 유지하고 시스템이 재시동하거나 지워질때까지 절대 지우지 안으면 어떨까? Linus는 "ramfs"이라 불리는 cache 시스템을 둘러싸는 작은 껍데기를 작성하였다. 다른 커널 개발자들은 "tmpfs"이라 불리는 개선된 버전을 만들었다. (이것은 데이터를 swap 공간에 저장할 수 있고 모든 가용 메모리를 다 소비하기 전에 공간이 다 채워지는 mount 지점의 크기를 제한한다. tmpfs의 예는 Initramfs이다.

이러한 ram 기반 filesystem들은 자동적으로 그들이 가지고 있는 데이터의 크기에 맞게 작아지거나 커진다. ramfs에 파일을 추가하거나 있는 파일의 용량을 키우는 작업은 자동적으로 더 많은 메모리를 할당한다. 그리고 파일을 지우거나 줄이는 작업은 그 메모리를 되돌려준다. 현시점에서 block device가 없기 때문에 block device와 cache사이에 어떠한 복사도 없다. cache에서의 복사는 데이터의 단순한 복사이다. 모든 훌룡한 부분들은 새로운 코드들이 아니다. 그러나 기존의 Linux caching 코드들을 위한 새로운 응용프로그램은 매우 단순하고 극도록 잘 검증된 기반구조에 기초한다.

root file system으로 initramfs을 사용하는 시스템은 filesystem으로 해석될 block device가 없기 때문에 하나의 filesystem 드라이버만으로 구현된 커널이면 안된다. 그져 메모리에 파일이 존재하게된다.

이 둘의 차이는 가상적인 disk와 filesystem의 차이이다. 서로 다른 차원에서 형성된 개념으로 비교가 불가능하다. 예를 들어 이야기 하자면 시스템의 RAM의 한부분을 RAMDISK로 사용하고 파일 시스템을 ramfs으로 포멧한다라고 하면 이해가 되려나?

이후의 내용은 다음 링크의 내용을 참조 하였음

http://www.thegeekstuff.com/2008/11/overview-of-ramfs-and-tmpfs-on-linux/

ramfs과 tmpfs의 차이

기본적으로 tmpfs과 ramfs는 몇가지 세세한 사항을 제외하고 거의 동일하다.

ramdisk는 동적으로 공간을 할당한다. 그레서 ramfs에 데이터를 쓰는 작업은 시스템에서 사용가능한 메모리의 크기를 넘어가지 않도록 해야한다. 이야기 해보자면 여러분이 시스템에 2G의 RAM을 가지고 있고 1G의 ramfs을 생성하고 /tmp/ram의 지점에 mount되어 있다고 하자. /tmp/ram의 전체 크기가 1G를 상회할 때도 여전히 데이터를 저장할 수 있다. 시스템은 1G를 넘는 데이터 저장 행위를 멈추지 못한다. 그러나 시스템의 RAM의 전체 크기 2G를 넘는 데이터 저장을 하게 되면 시스템은 멈추고 RAM에는 더이상 데이터를 저장할 공간이 없어지게 된다.

tmpfs은 동적으로 공간을 할당하지 않는다. tmpfs을 mount할 때 설정한 용량을 넘어서서 데이터 저장을 허용하지 않는다. 그레서 데이터를 쓰는 작업이 용량이 넘지 않도록 신경쓰지 않아도 된다. 하지만, 용량을 넘어서는 작업을 하게되면 “No space left on device”라는 메시지를 내보낸다.

tmpfs은 swap을 사용한다. (swap은 파일 시스템의 swap 디스크를 일컷는다.)

ramfs은 swap을 사용하지 않는다.

ramfs과 tmpfs의 단점

ramfs과 tmpfs의 단점은 시스템이 재시동하거나 멈추게 될 때 다 지워지게 된다. 그레서 주기적으로 disk에 저장하는 작업을 해주어야 한다. 시스템이 재시동할 때 ramfs과 tmpfs의 데이터를 disk에 저장하는 작업을 할 수 있지만 시스템이 멈출때는 이러한 추가적인 도움을 받을 수 없다.

Table: Comparison of ramfs and tmpfs

Experimentation Tmpfs Ramfs
허용량이상의 데이터 쓰기 에러 메시지 쓰기 작업 허용
고정된 총 용량 고정된 총크기 변동적인 총 크기
Swap 디스크 사용 사용 무사용
휘발성 저장 예 예

여러분이 데이터를 쓰는 작업을 좀 더 빠르게 하자고 한다면, 시스템 주의사항을 제공하는 tmpfs이 좀더 낳다.

RAM 업그레이드 후 Debian (wheezy) 시스템의 자원 사용 현황

2012-02-06T22:38:00.000-08:00

먼저 gnome-system-monitor에서 보여준 결과는 다음과 같다.

실행한 프로그램은 rythmbox, crome(active 30 tabs), virtualbox, 관리자 터미널, libreoffice등이다. 평소에 사용하는 조합이다. 가장 큰 가변 요인은 virtualbox의 램 사용률이다 설정은 1G로 설정해 놓았지만, 실제 그렇지는 않는 것 같다. 유동적으로 나눠 쓰는 모양인지 모르겠다. 아무튼 살펴본 봐는 위와 같다.

swap 파일 시스템은 전혀 사용하지 않는다. CPU는 어쩔 수 없지만 메모리는 1G정도 여유가 있는 것 같다. 그레서 crome의 cache 저장용으로 ramdisk나 tmpfs을 만들어서 사용해야 겠다.

우선 ramdisk와 tmpfs의 차이점에 대해서 살펴 봐야 겠다. 둘다 RAM을 디스크로 사용하는 기법인 것 같은데 구체적인 구현이 다른 것 같다. 그리고 둘을 지원하는지 살펴 보아야 겠다. 이것 당장 해 보았다.

tmpfs을 지원하는지 살펴보기
# grep tmpfs /proc/filesystems
nodev tmpfs
nodev devtmpfs

RAMDISK를 지원하는지 살펴보기

# dmesg | grep RAMDISK

[ 0.000000] RAMDISK: 36cc6000 - 3765b000

음. tmpfs이나 RAMDISK는 지원하는것 같다. 아래 출력은 이미 사용중이라는 생각이 들긴하지만 자세히 무엇을 의미하는지 모르겠다. 용량으로 살펴 볼땐 그래픽 모듈이 사용하는 것 같긴하다.

다음은 RAMDISK 만드는 법이다.
http://shibuvarkala.blogspot.com/2009/09/howto-create-ramdisk-in-ubuntu-debian.html
http://www.techtach.org/using-and-mounting-a-ram-disk-in-ubuntu-the-easy-way
http://securfox.wordpress.com/2009/03/14/how-to-create-a-ram-disk/

Ubuntu나 Fedora 배포판의 Firefox나 Crome에서 RAM을 cache로 설정하는 법
http://www.linuxreaders.com/2011/01/11/firefox-chrome-cache-on-ram-drive-fedora-ubuntu/

인터넷 전화기(LG U+ 070) 준비

2012-02-06T18:45:00.000-08:00

새것으로 구매하려고 하였으나 매달 행사로 나오는 물건들도 따지고 보면 많은 손해를 보도록 되어 있다.
삼년약정에 공유기 대여비 추가 부담.
그러니 행사가 아닌 물건은 도대체 누가 사는 것인지 궁금하다.

중고로 물건을 샀다. 평균가격이 3만원이었다만 여기저기 전화하고 그러는게 귀찬아 그냥 4만원짜리 샀다.
인내심을 시험하는 ARS 전화를 거쳐 중고기계를 사용한다는 사정 설명과 귀찮은 개인정보 동의 같은 것들.

그래서 결론은 월 2200원에 설치비 11000원이 든단다.
이놈의 설치비 라는 것은 전화한통으로 해결되는 mac주소 불러주는 것을 이야기하는 것이 었다.

음. 근데 동작을 안한다. MAC주소 받아 적는데 착오가 있었거나 뭔가 전산 착오라고 유추한다.
어쨌든 월요일에 설치기사가 온단다. 기다려 보자.

* mylg070로 기본 SSID가 설정되어 있는것 해제하여 다른 공유기를 사용하기
http://clien.career.co.kr/cs2/bbs/board.php?bo_table=lecture&wr_id=41561

Suspend to Ram에 생긴 그간의 문제 해결

2012-02-05T07:42:00.000-08:00

그 동안 debian wheezy의 hibernation은 잘 동작하였다.
그러니 suspend에 필요한 기본적인 동작들은 완벽하다.
하지만 Suspend to Ram에서 복구하는 기능은 마우스 커서가 살아있는 상태에서
멈춰서 버린다. 그래서 Ctrl+Alt+F1으로 관리자 터미널을 띄워서 재부팅하였다.
그리고 /var/log/pm-suspend.log의 기록들을 확인해 보았다.
다음과 같은 실패 항목들이 나타난다.

/etc/pm/sleep.d/10_unattended-upgrades-hibernate suspend suspend: success.
Running hook /usr/lib/pm-utils/sleep.d/55NetworkManager suspend suspend:
Having NetworkManager put all interaces to sleep...Failed.

/usr/lib/pm-utils/sleep.d/55NetworkManager suspend suspend: success.
Running hook /usr/lib/pm-utils/sleep.d/60_wpa_supplicant suspend suspend:
Failed to connect to wpa_supplicant - wpa_ctrl_open: No such file or directory

자세한 내용들을 찾아 보려고 googling을 해보았지만, 시원한 해결책은 없었다.
그져 네트워크 문제라고 유추해 볼 뿐이다.

네트워크의 가장 기초적인 설정을 살펴 보았다. 다음과 같이 되어있었다.

# This file describes the network interfaces available on your system
# and how to activate them. For more information, see interfaces(5).

# The loopback network interface
auto lo
iface lo inet loopback

# The primary network interface
allow-hotplug eth0
#NetworkManager#iface eth0 inet dhcp

#auto wlan0
#iface wlan0 inet dhcp
# wpa-ssid networkname
# wpa-psk password

아. 이 모든 것들은 필요에 따라 조금씩 수정해 보왔던 내 탓이다.

그래도 잘 동작했었다. 가끔씩 resume이 안되었을 때가 있었지만.

그리고 suspend to Ram은 잘 사용하지도 않는다.

그래서 다시 되돌려 놓았다. 문제되는 NetworkManager와 wpa-ssid, wpa-psk항목들을 제외하고, 그러니 되돌아 왔다.

# This file describes the network interfaces available on your system

# and how to activate them. For more information, see interfaces(5).

# The loopback network interface

auto lo

iface lo inet loopback

# The primary network interface

allow-hotplug eth0

#NetworkManager

iface eth0 inet dhcp

auto wlan0

#iface wlan0 inet dhcp

# wpa-ssid networkname

# wpa-psk password

이로써 Suspend to Ram도 잘 동작하였다.

미래에 대처하는 나의 자세

2012-02-02T09:27:00.000-08:00

나는 불안정하고
때로는 위로 받고 싶다.

가끔은 이런게 모자라 정체되곤한다.
그래서 바보가 되곤한다.

그래도 계속 가야한다.
그러니 가끔 멍하게 지내야한다.

Eigenvector(고유벡터) 와 Eigenvalue(고유값) 구하기

2012-01-24T17:52:00.000-08:00

책속에서 존재하는 eigenvector와 eigenvalue구하기의 난제는 다항식을 푸는 문제로 귀결된다.
하지만, 컴퓨터로 구현하고자 한다면 다항식을 푸는 문제 자체도 대단히 어려운 일이다.
해가 있는지 (허근인지 실근인지) 무엇보다도 소인수 분해와 같은 암호 해독과 같이 어려운 문제에 봉착하게 된다.
그래서 다항식 해법으로 접근하지 않고 eigenvetor와 eigenvalue의 property(속성)을 이용하여 반복적인 연산을 통하여 실제 eigenvector와 eigenvalue에 수렴할 수 있는 근접한 값을 찾아 내는 방법을 사용한다. 여기에는 여러가지 방법들이 존재한다. power iteration method, shifted inverse iteration method, Rayleigh quotient method, simultaneous iteration method, QR method
이와 같은 반복적인 연산을 이용한 방식은 eigenvalue와 eigenvector가 존재한다는 가정하에 다음과 같은 속성을 이용한다.

존재한다고 가정한 eigenvalue들은 {q_{i}}이고 eigenvalue들은 {lambda_{i}}이다.
여기서, | lambda_{1} | > | lambda_{2} | > ... > | lambda_{n} |의 조건을 갖는다.

그리고 추정하고자하는 eigenvector는 실제 eigenvector의 span으로 표현할 수 있다. 왜냐하면 eigenvector는 해당 공간의 basis가 되기 때문이다. 그래서 다음과 같이 eigenvector의 linear combination으로 표현할 수 있다.
v_{0} = c_{1} * q_{1} + ... + c_{n} * q_{n}.

이와같은 추정치가 eigenvector라고 한다면 실제의 eigenvector의 정의에 해당하는 조건식에도 부합되어야한다.

A * x = lambda *x 에서 부터
A * v_{0} = c_{1} * lambda_{1} * q_{1} + ... + c_{n} * lambda_{n} * q_{n}

그리고 A^k * x = lambda^k * x 의 확장된 속성으로 부터
A^k * v_{0} = c_{1} * lambda_{1}^k * q_{1} + ... + c_{n} * lambda_{n}^k * q_{n}
= lambda_{1}^k * [ c_{1} * q_{1} + c_{2} * (lambda_{2} / lambda_{1})^k * q_{2} + ... + c_{n} * (lambda_{n} / lambda_{1})^k * q_{n} ]

으로 전제할 수 있다. 여기서 주의 깊게 살펴보아야할 부분은 (lambda_{i} / lambda_{1})^k의 부분이다. 분모에 해당하는 lambda_{1}는 가장 큰 절대값을 갖는 eigenvalue이기 때문에 k값이 커질 수록 0의 값으로 수렴하게 된다. 그럼 나머지 항들은 사라지고 lambda_{1}에 해당하는 eigenvector의 값으로 v_{0}의 값은 수렴하게 된다. 기본적으로 이러한 개념하에 위의 반복적인 연산이 이루어진다.

그리고 이러한 기본 개념하에 power iteration방식과 shifted inverse iteration방식이 설명된다. 기타 다른 방식들도 이 방식에서 확장해 나아가는 형식을 취하고 있다.

power iteration
이 방식은 eigen value중 가장 큰 값에 해당하는 eigen vector만 구할수 있다.

* psudo code of power iteration
v = random vector with its norm "1"
do {
w = A * v;
if (w == v) break;
v = w / |w|;
} while (1);

inverse iteration
이 방식은 eigen value중 가장 작은 값에 해당하는 eigen vector만 구할 수 있다. power iteration 방식은 A의 eigen value가 lambda_{i}라면 A^{-1}의 eigen value이 1/lambda_{i}가 되는 특성을 이용하여 A^{-1}에 power iteration을 적용한다.

* psudo code of inverse iteration

v = random vector with its norm "1"
do {
w = A^{-1} * v;
if (w == v) break;
v = w / |w|;
} while (1);

shifted inverse iteration

A의 특정 eigen value에 근접하는 값을 임의로 선택한 값이 u라고 하면 (A - u * I)의 eigen value은 lambda_{i} - u가 된다. 그리고 (A - u * I)에 inverse iteration을 적용하면 u값에 근접하는 eigen value에 해당하는 eigen vector를 구할 수 있다.

* psudo code of shifted inverse iteration (finding eigen vector corresponding to eigen value closed to u)

v = random vector with its norm "1"
do {
w = (A - u * I)^{-1} * v;
if (w == v) break;
v = w / |w|;
} while (1);

Rayleigh quotient iteration

이 방식은 shifted inverse iteration에서 u값과 v값을 동시에 근사화 시켜서 연산 시간을 줄여주는 방식이다.

* psudo code of Rayleigh quotient iteration (finding eigen vector corresponding to eigen value closed to u)

v = random vector with its norm "1"
do {
w = (A - u * I)^{-1} * v;
if (w == v) break;
v = w / |w|;
u = v^{T} * A * v;
} while (1);

Simultaneous iteration

이 방식은 이전의 방식들이 한번에 하나의 eigen vector만을 구할 수 있지만 한번에 모두 구할 수 있다.

* psudo code of simultaneous iteration

V = [v_{i} : random vector with its norm "1" and i = 1..n];

// QR decomposite

Q * R = QR_decomposite(V);

do {

W = A * Q;

Q * R = QR_decomposite(W);

if (Q is not changed) break;

} while (1);

the columns of Q will converge towards a basis of eigenvectors of A

The QR method
이 방식은 QR 분해를 반복적으로 수행하는 것이다. 이론이고 뭐고 orthogonal column vector들은 항상 Q 행렬에 모이게 되고 Upper triangular 행렬은 R에 모이는 QR 분해를 Q 행렬에 반복적으로 recursive하게 적용하면 뭔가 유일한 orthogonal columns vector가 나올 법하게 느껴진다.

* psudo code of QR method

B = A;

do {

Q * R = QR_decomposite(B);

B = R * Q;
if (Q is not changed) break;

} while (1);

무모한 클릭질

2012-01-24T03:41:00.000-08:00

모니터에 불을 댕기고
수십번 클릭질을 해보아도
맘맞는 사람과 이야기하는 것만큼
허기가 풀리지는 않는구나

Must buy !

2012-01-15T06:12:00.000-08:00

잠비나이 - 나부락

이슈타르 - All of You ( 여자 찰스디킨스 )

쫄지마 ! (Don't panic !)

2011-12-30T08:18:00.001-08:00

Debian(wheezy)의 Gnome3 기본 테마 adwaita 변경하기

2011-12-12T17:50:00.000-08:00

# apt-get dist-upgrade후 gnome3의 기본 테마인 adwaita에 변화가 생겼다. taskbar 의 메뉴의 글자색이 taskbar와 같아져 버려 확인이 불가능하고 풍선 도움말의 글자색도 풍선 도움 바탕색과 같아서 마찬가지로 확인 불가능하다. 이래저래 기본 테마를 변경해야하는 상황이 되어 버렸다. 우선 gnome-look.org에 가서 적당한 gtk3 theme를 다운로드하고 ~/.themes밑에 복사한다. 그리고 테마를 변경하기 위해서 gnome-tweak-tool을 설치하자

# apt-get install gnome-tweak-tool

그럼 taskbar의 프로그램 -> 시스템도구 -> 기본설정 밑에 Tweak Advanced Settings라는 항목이 생긴다. 이를 실행하고 gtk+ 테마와 current theme를 다운로드한 테마로 변경하자.

내가 선택한 테마는 clearwaita-master라는 테마다.

그리고 Gnome Shell titlebar의 테마를 변경하기 위해서 Mutter라는 것을 건들여야한다. 나의 경우는 Gnome Classic으로 사용중이기 때문에 변경 하지 않아도 된다. 하지만 Gnome Classic 버전의 Metacity의 변경을 해주어야한다.

이와 같이 Mutter와 Metacity를 변경하기 위해선 gconf-editor를 사용해야한다. 설치되지 않았다면 설치하자.

# apt-get install gconf-editor

다음은 Mutter의 테마와 Metacity의 테마를 변경하는 상황이다.

해놓고 확인된 사항이지만 gnome-tweak-tool로 설정해 놓으면 위의 두 사항은 자동으로 설정되는 것 같다. 막상 변경하려고 gconf-editor를 실행 시켜 보았더니 이미 변경되어 있었다.

Gnome3 의 바탕화면에서의 아이콘과 마우스 오른쪽 버튼 살리기

2011-12-12T02:40:00.000-08:00

# dconf-editor

바탕화면에 아이콘이 보고 싶다면 다음 사항을 체크하고

org -> gnome -> background 밑에 show-desktop-icons [V]

쓰레기통이 보고 싶다면 다음 사항도 체크하자.

org -> nautilus -> desktop밑에 trash-icon-visible [V]

[펌] 한미 FTA; Raising the Flag on Korea

2011-12-06T23:36:00.001-08:00

출처 : http://clien.career.co.kr/cs2/bbs/board.php?bo_table=image&wr_id=3647134

Xorg의 touchpad설정하기 (클릭과 스크롤 설정하기)

2011-12-01T07:37:00.001-08:00

http://www.reddeadresolve.com/2011/01/configure-synaptic-touchpad-enable.html

Debian Wheezy에서 사라졌었던 fglrx-driver와 fglrx-control 다시 깔기

2011-11-30T03:30:00.001-08:00

$ sudo apt-get remove --purge xserver-xorg-video-radeon
$ sudo apt-get install fglrx-driver fglrx-control fglrx-glx fglrx-atieventsd fglrx-modules-dkms

그리고

$ sudo aticonfig -initial 로 xorg.conf파일을 /etc/X11밑에 생성 시킨다.

만약 laptop의 경우 dual헤드를 사용하여 모니터를 두 개 사용하고자 한다면

다음과 같은 방식으로 xorg.conf파일을 생성하자.

$ sudo aticonfig --initial=dual-head --screen-layout=right -v

새로운 3CCD 영상 분할 프리즘 구조

2011-11-25T01:58:00.001-08:00

출처 : http://www.nitto-optical.co.jp/english/products/basic_prism/index.html

참조 : http://ntuzhchen.blogspot.com/2011/01/toshiba-3ccd.html

[번역] Multispectral line-scan imaging system for simultaneous fluorescence and reflectance measurements of apples: multitask apple inspection system

2011-11-24T18:42:00.001-08:00

Multispectral line-scan imaging system for simultaneous
fluorescence and reflectance measurements of apples: multitask
apple inspection system

Moon Sung Kim, Kangjin Lee, Kuanglin Chao, Alan M. Lefcourt,
Won Jun, Diane E. Chan

이 방법론적 논문에서는 고속으로 움직이는 물체로부터 다파장 반사및 형광 영상을 동시에 획득할 수 있는 최근 개발된 라인-스캔 방식의 영상 시스템이 소개 한다.

시스템은 다파장 영역의 유동적인 선택으로 사과의 품질과 안정성에 대한 다중 검사 작업을 할 수 있도록 해준다.

상업적인 사과 분류 기계에 장착된 시스템은 초당 3 ~ 4개의 선별 속도에서 사과의 영상 해상도를 결정하도록 개발되었다.

분류 시스템에 적재된 사과들은 약 15cm의 간격을 갖도록 하였다.

이런 인자들을 가지고 결과 영상들은 사과당 40개의 라인스캔 영상들을 보여준다.

수평 방향 해상도는 픽셀당 2mm이고, 수직 방향으로 각 라인 스캔 영상당 1002픽셀이고,

시스템의 공간 해상도는 binning의 선택에 의존하는데 픽셀당 0.2mm 정도이다.

종합적인 공간 해상도는 사과의 표면 오염이나 결함의 온라인 영상 기반 검사에 적합하고 이전의 우리의 연구에 비해서도 구분된다.

Keywords
Online inspection, Multispectral imaging, Imaging system, Fluorescence, Reflectance, Food safety and quality, Apples

경향신문 1면 FTA 찬성자 명단

2011-11-24T01:07:00.001-08:00

[번역] LINE‐SCAN HYPERSPECTRAL IMAGING PLATFORM FOR AGRO‐FOOD SAFETY AND QUALITY EVALUATION: SYSTEM ENHANCEMENT AND CHARACTERIZATION

2011-11-23T23:06:00.001-08:00

LINE‐SCAN HYPERSPECTRAL IMAGING PLATFORM FOR
AGRO‐FOOD SAFETY AND QUALITY EVALUATION:
SYSTEM ENHANCEMENT AND CHARACTERIZATION

M. S. Kim, K. Chao, D. E. Chan, W. Jun, A. M. Lefcourt, S. R. Delwiche, S. Kang, K. Lee

Keywords. Fluorescence, Hyperspectral, Line?scan imaging, Reflectance, Spectral calibration.

라인 스캔방식의 초분광 영상화 기술들은 고속 온라인 응용분야에 대한 단 몇 가지의 파장대에 기반한 고속 다파장 기법 개발에 대한 연구 도구로 자주 제안되어 진다.

지속적인 기술적 발전과 광전자 영상화 소자(optoelectronic imaging sensors)와 파장 영상용 파장 분해 장치(spectral imaging spectrographs)에 대한 좀 더 낳은 접근성, 유용성이 증대 됨에 따라서

초분광 영상 장치의 적용 분야는 식품및 농업 산업에서 품질및 안정성 검사 수요에 걸쳐 확대되어 왔다.

우리는 적은 빛에도 감응도를 갖는 EMCCD를 사용한 초분광 라인 스캔 영상 장치에 기반한 식품 검사 영상 시스템을 계속해서 개발해 왔다.

이 방법론적 기사에서 ARS 초분광 영상 시스템의 가장 최근 세대(가시 및 근적외 파장대 영역 400 - 1000nm에서 반사와 형광 측정이 가능한)의 파장 및 공간적 성능이 평가 되었다.

결과는 시스템의 파장 분해능은 우리의 일반적인 동작 모드 (6 픽셀의 파장 binning)에서 6nm의 FWHM과 FWHM 구간에서 4.4nm의 간격을 보인다.

우리는 A/D 변환의 동적인 범위를 더 활용하기 위해서 파장공간에서의 가중 필터를 사용한 시스템 출력 응답성을 개선하였다.

이러한 시스템 출력 조정으로 표준 방사 반사체의 영역에 대한 근적외선 영역 대에서의 잡음으로 추정되는 지수를 측정한 결과 50%정도의 감소 되었다.

[번역] DETECTION OF ORGANIC RESIDUES ON POULTRY PROCESSING EQUIPMENT SURFACES BY LED‐INDUCED FLUORESCENCE IMAGING

2011-11-23T20:53:00.001-08:00

DETECTION OF ORGANIC RESIDUES ON POULTRY PROCESSING
EQUIPMENT SURFACES BY LED‐INDUCED
FLUORESCENCE IMAGING

J. Qin, K. Chao, M. S. Kim, S. Kang, B.‐K. Cho, W. Jun

가금류 처리 공장에서 장치 표면에 있는 Organic residues(유기물 찌꺼기)는 교차오염을 유발 시킬 수 있고 소비자에게 안전하지 않은 음식에 대한 위험을 증가 시킨다.

이 연구는 가금류 처리 장비 표면에 있는 유기물 찌꺼기 신속한 검사를 위한 LED로 생성된 형광 영상 기술의 가능성을 조사하는 목적을 가지고 있다.

400nm의 파장대의 출력을 갖는 고출력 파란색 LED가 라인 스캔 방식의 초분광 영상 시스템을 위한 활성화 반응의 광원으로 사용되어 졌다.

지방, 피, (결장, 십이지장, 장)배설물을 포함한 공통적인 닭의 찌꺼기 표본들은 스테인레스 철판에 준비하였다.

형광 방출 영상들은 500에서 700nm의 파장대에서 각 찌꺼기의 종류별 20개씩 총 120개의 표본의 영상이다.

실험에 사용된 표본의 LED광원 형광 특성들이 판정 되었다.

주성분 분석(PCA)은 형광 스펙트럼 데이터를 분석에 사용되었다.

두 개의 SIMCA 모델들이 유기물 찌꺼기와 스테인레스 표본을 구분하는데 개발되었다.

유기물 찌꺼기와 스테인레스 강의 2 부류와 스테인레스강, 지방, 피, 배설물등의 4부류에 대한 SIMCA 모델의 분류 정확도는 각각 100%와 97.5%였다.

666nm로 선택된 파장대를 이용한 최적 단일 밴드 접근은 닭의 피 검사에 있어서 false negative errors를 생성할 수 있었다.

503과 666nm의 파장대를 이용한 두 파장의 비율(F503/F666) 영상은 스테인레스 강 표면에 다양한 닭 찌꺼지를 검출하는데 가능성을 보인다.

이 파장대 조합은 가금류 처리 장치 표면 검사를 위해서 LED 기반의 휴대용 형광 영상화 장치 개발에 적용 될 수 있다.

keywords. Food safety, Sanitation monitoring, Hyperspectral imaging, Fluorescence, LED, Poultry.

한글 번역 참조

2011-11-21T03:31:00.001-08:00

HowTO : 도움글
Residues : 찌꺼기
Optoelectronic Imaging Sensors: 광전자 영상화 소자
Spectral Imaging Spectrographs: 파장 영상용 파장 분해 장치
Push Broom : 밀대
KnowHow: 비법
Threshold Value: 문턱값

ASUS 1215T 여러가지 문제점들

2011-11-19T23:57:00.001-08:00

전에는 touchpad 동작에 문제가 많아 분해하는 과정에서 유격이 발생하고 몇몇 힌지가 부러져 버렸다. 문제 해결은 touchpad의 금박에 접촉되는 두개의 패드만 제거 하면 되었다.
그리고 그 이후 데비안에서 xorg.conf를 변경하고 테스트 하던중 마우스가 혼자서 춤을 추길래 재부팅을 하였다. 화면이 나오질 않았다. 부팅과정은 외장 모니터를 연결해본 결과 정상적으로 동작하는 것을 확인 할 수 있었다. 우선 변경하였던 Xorg.conf를 원상 복구하고 외장 모니터를 연결해서 원상 복구 시켰다. 하지만, 노트북 lcd는 복구 불가능 한 것 같아 보였다. ASUS A/S에 갔다 첫날엔 부품 없다고 다음날 오라는 통보를 받았다. 다음날 갔다. 1시간 이상 걸린것 같다. 엔지니어 분도 미안했는지 LCD 연결 케이블 까지 같이 갈아 드리려다 보니 시간이 지체되었다고 하시며 교체한 케이블도 보여주었다. 우선은 잘 동작하는 것 같아 만족스럽게 가지고 왔다.

몇일 사용하던중 노트북 Lid "closed"가 동작하지 않는 것을 발견하였다. 뚜껑이 닫혀도 화면은 여전히 밝게 나오고 있었다. 우선은 ACPI의 동작을 의심해 보았다. 코딩으로 LID의 동작을 모니터링 해 보았다.

#include

void main(void)
{
char state[80];

FILE *fd;

while (1) {
fd = fopen("/proc/acpi/button/lid/LID/state", "r");
fgets(state, 80, fd);
printf("%s", state);
fclose(fd);
sleep(1);
}
}

위와 같이 코딩하고 다음과 같이 컴파일하고 실행 하였다.

# gcc test_lid.c

# ./a.out

그리고 뚜껑을 닫아 보았다. 여전히 "state: open"이라는 메시지만 나온다.

하드웨어 적인 문제로 의심을 확정하게 되었다. 우선 뚜껑 주변에 스위치를 찾아 보았다.

어디에도 뚜껑이 닫히는 것을 감지할만한 스위치는 없었다. 검색해 보니 LCD하단에 자석이 있어야 동작을 한단다.

출처: http://www.insidemylaptop.com/images/Asus-Eee-PC-1201n-screen/replace-damaged-laptop-screen-11.jpg

위 그림의 노란색 화살표가 바로 뚜껑이 동작하게 해주는 자석이다. 뜯어 보니 자석인 사진의 위치보다 한참 위에 있었다. A/S기사의 실수 였다. 핀셋으로 위치를 내렸지만 볼트 머리때문에 사진의 위치까지는 내려가지 않았다. 하지만, 동작은 하였다. 그냥 덥었다. 하도 뜯었다 조립했다를 반복했더니 유격이 있다. 에휴,. 깨끗하게 쓰고 싶었는데 걸래가 되었다.

이질적인 우분투 11.10의 가장먼저 설치할 것들

2011-11-13T08:46:00.001-08:00

# sudo apt-get install gnome-session-fallback gnome-tweak-tool

로그인시 그놈 클래식으로 선택하고 gnome-tweak-tool에서 눈 안아픈 테마를 설치하자

추천하는 테마는 adwait-white 인가하는 것이다.

막가파 불도저의 무의식이란 ?

2011-11-04T06:30:00.000-07:00

무의식이란 단어의 의미는 의식이 없다는 뜻이다.

이 말은 본질적으로 잘못된 말이다.

굳이 이야기 하자면 자각과 무자각의 대립상태를 말하는 것이다.

이러한 무의식의 개념은 프로이드가 많이 연구하였다.

하지만 여기서 이야기하고자 하는 내용은

에리히 프롬의 무의식에 대한 재해석을 기반으로 한다.

나의 경험이다

커다란 일을 시작하다 보면 초반에 여러가지의 문제에 봉착하기도 하고

이 일의 당위성에 위배되는 많은 일들을 겪게 된다.

그래서 아주 짜증나고 골치아픈 일이지만

그러한 일들을 피해갈 수 있는 길들을 고려하면서 일을 진행한다.

그렇게 시간이 지나고 일이 끝날때 쯤 생각을 정리해 보면

위와 같은 일들은 까맣게 잃어 버리고 있었던 나를 발견한다.

싫은 일은 보지 않고 넘어간다는 인간의 능력,

즉 인간은 인식하려고하면 인식할 수 있는 일을

진지한 얼굴로 부정하는 경우도 있다는 사실을 잊었던 것이다.
(H.S. 설리번은 이현상을 매우 훌룡히 표현하여 '선택적 부주의'라는 술어를 만들었다.)

나의 경우는 일을 깔끔하고 추진력있게 진행해야한다는 압박 때문에

이러한 선택적 부주의가 작동한 것같다.

지금에와서 자각해 버렸지만, 지나온 일들은 우려하는 목소리를 전혀 듣지 않았다.

그져 미안한 마음뿐이다. 정말 악의적이지는 않았다.

우리의 가카도 이와 같지 않았나 하는 추정을 해본다.

처음에 반대의 목소리를 귀기울여 듣겠다고 명심했지만,

꼭 임기 중에 해치워야 하는데 하는 압박감이

'선택적 부주의'라는 매카니즘을 작동 시켰다고.

그래서 '유체이탈 화법'은 의도적인 것이 아닌것이다.

그러한 면에서 나도 조심해야 하겠다고 생각한다.

반대의 목소리를 선택적으로 부주의하게 하지 않도록....

하지만, 우리나라에서 성공했다는 그 많은 사람들은

알게 모르게 위와 같은 '선택적 부주의' 매카니즘을 효율적으로 사용하였다.

사실인지 모르겠지만 공식적으로 이야기 한다.

반대의 안된다는 목소리는 신경쓰지 않았다고.

뭐 하고 싶은 일이고 당위성이 있었다면 문제가 없는 건가 ?

왜 국가는 왜 주식시장에 민감한가 ?

2011-11-04T03:32:00.000-07:00

실업율보다 더
물가 보다 더
양극화 지수보다 더
주택비용보다 더
빈민층의 비율 보다 더

신경쓰는 것 같다.

반대로 생각해야 하는게 맞는거 아닌가.

Enable MSN, yahoo, etc in Empathy (gnome default messenger)

2011-11-03T10:43:00.001-07:00

# apt-get install telepathy-haze

Install Mobile Atlas Creator on Debian

2011-11-02T21:24:00.000-07:00

Mobile Atlas Creator is offline atlases for GPS handhelds and cell phone applications.

Installation
# su
# apt-get install mobile-atlas-creator

Execution
# su
# mobile-atlas-creator

EEEPC(1215t) hang on with Caplock led blinking in Debian Wheezy.

2011-10-31T08:09:00.000-07:00

This symptom is because of kernel panic crashing with wireless module.

# lspci
...

02:00.0 Network controller: Broadcom Corporation BCM4313 802.11b/g/n Wireless LAN Controller (rev 01)

...

As you can see, my wireless card is named as "BCM4313 802.11b/g/n Wireless LAN".

Googling provides that it is necessary to upgrade linux-kernel as 3.0.x.

자전거란?

2011-10-28T09:06:00.000-07:00

나의 경우 자전거를 사용하게된 계기는

내가 생활하는 영역안에서 조금이나마 자동차의 사용을 줄이려는 목적에 의해서 였다.

다른 사람들도 이와 크게 다르지 않을 것이라 생각한다.

등교 길에서 출퇴근 길에서 자전거 도로가 구비 되고 환경이 갖추어진다면 좋다.

그런 의미에서 "MB의 자전거 길이란 ?"

주말에 잠시 자전거를 타기 위해서 만든

마치 놀이 공원에 가서 자전거를 잠시 타고온 듯한 효과를 노려

놀이 공원이 아주 친환경적인 듯한 착각을 불러 잃으키게 만드는

자전거의 용도를 오직 레져로서만 바라본

극히 끼워 맞춘듯한 그 많은 공사에 대한 변명거리의 용도로 밖에 볼 수 없다.

이와 같은 예는 얼마든지 찾아 볼 수 있다.

핵연료가 단지 이산화탄소만 배출하지 않는 다는 이유로

폐기물 처리에 들어가는 환경적인 피해나 관리 비용에 대한 언급은 없이

친환경 에너지라고 선전하는 것이라던가.

수소에너지에 대한 선전을 위해서 항상 예를 들던 맹물 자동차에 대한 허풍이라던가.

물에서 수소를 분리해내기 위해선

얼마나 많은 양의 석유로 발전된 전기로 분해내야 한다는 것은 전혀 고려하지 않는다.

이런식이다. 항상.

After removing AMD/ATI Catalyst radeon driver, experimental trials for fixing problems드라이버의 문제해결을 위한 혹시나 하는 재설정

2011-10-26T01:51:00.000-07:00

# apt-get update

# apt-get install --reinstall libgl1-mesa-glx libgl1-mesa-dri <--- libGL.so + mesa3d DRI drivers

# apt-get install --reinstall xserver-xorg-core <--- libglx.so of xserver

# apt-get build-dep libdrm <--- Userspace interface to kernel DRM services (development files, radeon-specific, etc.)

# apt-get build-dep mesa <--- A free implementation of the OpenGL API (DRI modules, GLX runtime, GLX development files, etc.)

# apt-get build-dep xserver-xorg-video-ati <--- X.Org X server -- AMD/ATI display driver

# apt-get build-dep xorg-server <--- Xorg X server (core server, development files, etc.)

Reference

http://www.x.org/wiki/radeonBuildHowTo#Removing_AMD.2BAC8-ATI_catalyst.2BAC8-fglrx_driver_.28closed_source.29

ATI radeon 비디오 카드를 위한 Kernel ModeSetting (KMS) 비활성화 하기

2011-10-26T01:20:00.000-07:00

부팅시 "nomodeset"를 설정해서 시스템을 가동하고

/etc/modprobe.d/radeon-kms.conf를 편집한다.

"options radeon modeset=0"

Debian Wheezy(7.0) 에서 suspend 동작 살펴보기(debugging 하기)

2011-10-26T01:02:00.000-07:00

[출처] http://wiki.debian.org/Suspend

Wheezy는 아직도 testing 상태에 있다. 그레서 그 구성은 언제든지 갑자기 변경될 수 있다.

가장 주목할 만한 변화는 HAL 패키지가 단계적으로 배제된다는 것이다. 여러분이 아직도 HAL 패키지가 설치된 상태로 있다면, 제거해야 합니다. 그렇지 않으면 suspend상태 동안에 pm-utils과 충돌하게 될 겁니다.

suspend / resume 이 여러분의 시스템에서 잘 동작한다면, 여러분은 다행이도 이 페이지의 어떤 것도 읽을 것이 없습니다. 그렇지 않다면 맨 먼저 debug 해야 할 것은suspend와 resume을 제어하는 pm-utils의 debug를 활성화 하는 것이다.

pm-utils의 debug를 활성화 하기
suspend와 resume process의 동작 기록은 /var/log/pm-suspend.log 파일에 있다. 이것은 초기설정으로 절제된 동작 정보를 가지고 있다. 더 많은 정보는 /usr/lib/pm-utils/pm-functions 파일의 앞부분에 "export PM_DEBUG=true"를 추가해서 활성화된 debug를 통해서 활성화 된다.

resume시 오류가 있는 video화면 고치기
컴퓨터가 resume 후에 생기는 가장 공통적인 문제는 오류가 있는 video나 검은 화면, 백라이트가 없는 LCD화면이 나타나는 것이다. 먼저 점검해 봐야 할 점은 시스템이 여전히 동작하는지 확인해 보는 것이다. 이건 "Capslock"버튼을 눌러서 간단히 확인해 볼 수 있다. Capslock LED가 적절히 변화는지를 점검하면 된다. 시스템이 여전히 동작한다면, 대부분의 경우 여러분의 비디오 카드를 위한 "video quirk"를 추가하는 것이 필요하다.

Debian은 현재 Intel, nVidia와 ATI 비디오 카드에 대해서 초기적으로 활성화된 kernel mode setting(KMS)을 가지고 있다. 그러나 pm-utils의 video quirk는 KMS를 이미 지원하고 있다. 그래서 대부분의 경우 여러분은 KMS를 먼저 비활성화하는 시도를 해야한다. 여러분의 특정 비디오 카드에 대한 세세한 내용은 온라인을 검색해서 쉽게 찾아 볼 수 있다.

KMS를 비활성화한 후에도 resume후 video에 여전히 오류가 있다면, 여러분은 어떤 video quirk를 사용함으로써 시스템이 suspend하도록 시도할 수 있다. 모든 사용 가능한 quirk의 자세한 설명은 pm-suspend 프로그램의 manpage를 읽어라. 그리고 명령행에서 그들의 조합을 시도하라. 여러분이 성공적으로 시스템에서 잘 동작하는 quirk의 조합을 찾았다면, 영구적으로 동작하도록 /usr/lib/pm-utils/video-quirks에 추가 할 수 있다. 그리고 동시에 많은 사람들에게 도움이 될 수 있도록 여러분의 변경사항에 대한 patch로 pm-utils 패키지의 bug 파일에 도움을 달라.

Debian의 예전 버전에서 upgrade된 시스템에서 생기는 공통된 문제점은 suspend동안 시스템을 멈추어 버리게 하는 quirk-s3-bios의 활성화 되는 것이다. 여러분의 시스템이 suspend동안 멈추어 버린다면, debug를 활성화 한후에 pm-suspend.log 파일을 주의 깊게 점검하자. 그리고 quirk-s3-bios를 사용하지 않는지를 확인하도록 하자.

10.26. 투표 안내문 없이 투표장 찾기

2011-10-25T21:35:00.000-07:00

마이크로 소프트 사의 인터넷 익스플로에서

http://www.nec.go.kr/

위의 주소로 이동한다.

다음과 같은 화면에서 빨간 네모의 "내 투표소 찾기"를 선택한다.

내 거주 지역의 주소와 이름 주민번호를 입력하면 해당 투표소의 주소와 지도를볼수 있다.

Orthogonal Partial Least Squares Calibration

2011-10-18T09:42:00.000-07:00

* Expression of matrix

(MATRIX:[n x m]) ; n is the number of rows and m is the number of columns.

* Expression of column vector

(MATRIX)(j,*) ; j-th column vector of MATRIX

* Expression of row vector

(MATRIX)(*,i) ; i-th row vector of MATRIX

* Operation for matrix

(MATRIX)' ; transpose of Matrix
(MATRIX)^(-1) ; inverse matrix

* Algorithm

Inputs

n ; the number of samples
m ; the number of dimension of independent value
(X:[n x m]) ; independent values
(y:[n x 1]) ; dependent values
interation ; the number of iteration for calibration

Step 1.

Mean centering
(mx:[1 x m]) ; Mean of (X:[n x m])(*,i) for i = 0 ... n-1
(X_{0}:[n x m])(*,i) = (X:[n x m])(*,i) - (mx:[1 x m]) for i = 0 ... n-1

Step 2.

Loading Weight
(w_{a}:[m x 1]) = (X_{a-1}:[n x m])' * (y_{a-1}:[n x 1]) /
sqrt((y_{a-1}:[n x 1])' * (X_{a-1}:[n x m]) * (X_{a-1}:[n x m])' * (y_{a-1}:[n x 1]))

Step 3.

Scoring
(t_{a}:[n x 1]) = (X_{a-1}:[n x m]) * (w_{a}:[m x 1])

Step 4.
Spectral Loading
(p_{a}:[m x 1]) = (X_{a-1}:[n x m])' * (t_{a}:[n x 1]) / ((t_{a}:[n x 1])' * (t_{a}:[n x 1]))

Step 5.

Chemical Loading
(q_{a}:[1 x 1]) = (y_{a-1}:[n x 1])' * (t_{a}:[n x 1]) / ((t_{a}:[n x 1])' * (t_{a}:[n x 1]))

Step 6.

Updating
(X_{a}:[n x m]) = (X_{a-1}:[n x m]) - (t_{a}:[n x 1]) * (p_{a}:[m x 1])'
(y_{a}:[n x 1]) = (y_{a-1}:[n x 1]) - (t_{a}:[n x 1]) * (q_{a}:[1 x 1])
a = a + 1

Step 7.

Until 'a' is equal to "iteration" repeat Step 2, 3, 4, 5, and 6.

Step 8.

(b:[m x 1]) = (W:[m x iteration]) * ((P:[m x iteration])' * (W:[m x iteration]))^(-1) * (Q:[iteration x 1])
(b0:[1 x 1]) = (my:[1 x 1]) - (mx:[1 x m]) * (b:[m x 1])

where

(W:[m x iteration])(a,*) = (w_{a}:[m x 1])
(P:[m x iteration])(a, *) = (p_{a}:[m x 1])
(Q:[iteration x 1])(a, *) = (q_{a}:[1 x 1])
(my:[1 x 1]) ; Mean of (y:[m x 1])(0,i) for i = 0 ... n-1

* Prediction

(y^:[1 x 1]) = (b0:[1 x 1]) + (x:[1 x m]) * (b:[m x 1])

Debian Wheezy(7.0)의 repository에서 사라진 fglrx-driver, fglrx-control ???

2011-10-18T04:30:00.000-07:00

여느날과 다름없이 업데이트 알림에 대해 의심없이 찬성을 했다가 다음 부팅시 Xorg 작동이 멈추어서 터미널로 로그인하게 되었다. 우선 /etc/X11/xorg.conf를 지우고 부팅하면 작동은 되었다만 여간 답답한게 아니다. 동영상도 느리고 화면 반응도 답답하다. 찾아 본 봐 AMD proprierty 드라이버가 충돌이 잃어난다는 보고에 패키지가 내려갔다고 한다. 그 동안은 xorg-server-radeon으로 버텨야 해야한다. http://web.archiveorange.com/archive/v/rzXEMJaA6QnUnWIezYjn 에서 알 수 있듯이 fglrx-driver 11-10이 나온후 packaging과 migration 을 거친후 11월 중순에나 만나 볼 수 있겠다고 한다.

MingW + GIT + Emacs

2011-10-06T19:24:00.000-07:00

Previous post mentioned about the installation of MingW.
And I add further setting to use Emacs;

To use git tool with MingW in windows XP I installed windows-ported version.

http://code.google.com/p/msysgit/

Then there are two the same shell commanders; Git Bash and Msys Bash contained in MingW package.
I would like to use Git under Msys Bash. However, in Msys Bash I can't access Git binary.
So, I need the below procedure.

1. Create .profile file in your home directory

After installation of Msys, set up your .profile something like this.
Create a .profile file in your home directory

Put the content needed into it. An example might be:

#!/bin/sh.exe
#
# Get the aliases and functions
#
if [ -f ${HOME}/.bashrc ]
then
. ${HOME}/.bashrc
fi

export CVS_RSH="ssh"
export EDITOR="/usr/bin/vim"
export HISTSIZE="500"
export PATH="/c/emacs-23.3/bin:/c/Program\ Files/Git/bin:${HOME}/bin:${PATH}"

2. Create a .bashrc file in your home directory :

Put the content needed into it. An example might be:

# .bashrc

# Source global definitions
if [ -f /etc/bashrc ]
then
. /etc/bashrc
fi

# User specific aliases and functions
#

Close and restart the MSYS Shell.

Install MingW and SDL.

2011-09-25T22:36:00.000-07:00

* Goto http://sourceforge.net/projects/mingw/files/Automated%20MinGW%20Installer/mingw-get/

* Download mingw-get-inst-xxxxxx.exe

* Execute MingW Shell

* Execute "mingw-get update" and "mingw-get upgrade"

* Goto http://www.libsdl.org/download-1.2.php

* Download SDL-X.X.X-win32.zip under Win32 of Runtime Libraries and SDL-devel-X.X.X-mingw32.tar.gz under Win32 of Development Libraries.

* Extract SDL-devel-X.X.X-win32.zip.

* Copy libSDLmain.a, libSDL.la, and libSDL.dll.a into into $(MingW Directory)/lib/*.a.

* Copy include/SDL/* into $(MingW Directory)/include/SDL.

* Copy bin/SDL.dll and bin/sdl-config into $(MingW Directory)/bin.

Using delaunay trianglulation, thinning of a mass object.

2011-09-25T03:08:00.000-07:00

Thinning algorithm based on morphology operation is not adopted for a mass object on image. It needs too many operations. To down-size a processing time one of the other methods is based on delaunay triangulation.

For a triangle having a edge in outside of a mass object, it is out of consideration.

Blue-colored triangles are out of consideration. And then it needs analysis the relationship of connectivity between triangles. A few definitions are enlisted. the End Triangle is defined as two edges of triangle meet outside of mass object, the Edge Triangle is defined as one edge of triangle meets outside of mass object, and the Branch Triangle is defined as none edge of triangle meets outside of mass object.

Blue-colored triangle indicates a End Triangle, Green-colored triangle indicates a Edge Triangle, and red-colored triangle indicates a Branch Triangle. According to connectivity among three types of triangles, the center of innercircle of End and Branch Triangle forms thinning lines.

Fixing delaunay triangulation

2011-09-21T09:15:00.000-07:00

During delaunay triangulation unexpacted line crossing is performed.
Given points are like this.

And Its result is

Its code is

void delaunay_list_create_2d_network(delaunay_list_t *list, point_list_t *points, bitmap_t *roi)
{
dlink_t *x, *y, *z;
point_t *o, *a, *b, *c;
point_t *v1, *v2, *v3, *org;
delaunay_net_t *neta, *netb, *netc, *net;
//dlist_t *triangles;
triangle_list_t *triangles;
triangle_t *tri;
//dlist_t *lines;
line_list_t *lines;
line_t *line, *l;
real_t da, db, dc, s;
real_t val, vmin;
real_t radius, area;

assert(list);
assert(points);

// allocate delaunay network list
//list = delaunay_list_new();
for (x = points->tail->next; x != points->head; x = x->next) {
net = delaunay_net_new((point_t *)x->object);
delaunay_list_insert(net, list);
}

// starting point
neta = delaunay_list_get_net(0, list);
a = delaunay_net_get_coordinate(neta);
// find minimum distance between points
for (netb = NULL, x = list->tail->next->next; x != list->head; x = x->next) {
net = (delaunay_net_t *)x->object;
o = delaunay_net_get_coordinate(net);
val = distance_between_points(a, o);
if (val < REAL_EPSILON) continue;
if (netb == NULL || val < vmin) { netb = net; vmin = val; }
}

b = delaunay_net_get_coordinate(netb);
// find minimum radius of circumscribed circle for triangle
for (netc = NULL, x = list->tail->next->next; x != list->head; x = x->next) {
net = (delaunay_net_t *)x->object;
if (net == netb) continue;
o = delaunay_net_get_coordinate(net);
// Heron's formula to find the area formed by three points
da = distance_between_points(a, b);
db = distance_between_points(b, o);
dc = distance_between_points(o, a);
s = (da + db + dc) / 2;
area = sqrt(s * (s - da) * (s - db) * (s - dc));
if (area < REAL_EPSILON) continue;
// calculate radius of circumscribed circle for triangle
radius = (da * db * dc) / (4 * area);
if (netc == NULL || radius < vmin) { netc = net; vmin = radius; }
}
c = delaunay_net_get_coordinate(netc);

v1 = point_new();
v2 = point_new();
v3 = point_new();
org = point_new();

// re-ordering
point_subtract(v1, b, a);
point_subtract(v2, c, a);
point_xproduct(v3, v1, v2);
if (v3->z > 0) { // counterclockwise following a to b to c ?
net = netb, netb = netc, netc = net;
o = b, b = c, c = o;
}

// connect delaunay network
delaunay_net_connect(neta, netb);
delaunay_net_connect(netb, netc);
delaunay_net_connect(netc, neta);

// Synthesize a triangle and insert it.
//triangles = dlist_new();
triangles = triangle_list_new();
tri = triangle_new_and_assign(a, b, c);
triangle_list_insert(tri, triangles);

//triangle_dump(tri);
// Synthesize three line and insert it.
lines = line_list_new();
line = line_new_and_assign(a, b);
line_list_insert(line, lines);
line = line_new_and_assign(b, c);
line_list_insert(line, lines);
line = line_new_and_assign(c, a);
line_list_insert(line, lines);

// main routine
for (x = lines->tail->next; x != lines->head; x = x->next) {
line = (line_t *)x->object;
//line_dump(line);
for (y = list->tail->next; y != list->head; y = y->next) {
netc = (delaunay_net_t *)y->object;
c = delaunay_net_get_coordinate(netc);
if (point_cmp(c, line->a) == 0 || point_cmp(c, line->b) == 0)
continue;
// counter-clockwise rotation of the triangle formed by a, b, and c
point_subtract(v1, line->b, line->a);
point_subtract(v2, c, line->a);
point_xproduct(v3, v1, v2);
if (v3->z < REAL_EPSILON) continue;
// Does the triangle exist in the pool of delaunay triangles before ?
tri = triangle_new_and_assign(line->a, c, line->b);
//if (is_triangle_in_list(tri, triangles)) {
if (triangle_list_query(tri, triangles) >= 0) {
triangle_destroy(tri);
tri = NULL;
continue;
}
// Does the unique triangle form the delaunay triangle ?
radius = triangle_circumcircle(org, tri);
for (z = list->tail->next; z != list->head; z = z->next) {
net = (delaunay_net_t *)z->object;
o = delaunay_net_get_coordinate(net);
if (point_cmp(o, line->a) == 0 ||
point_cmp(o, line->b) == 0 ||
point_cmp(o, c) == 0)
continue;
if (distance_between_points(o, org) < radius) {
triangle_destroy(tri);
tri = NULL;
break;
}
}

if (tri) {
/* put the trinangle into triangle's list,
* establish connection between delaunay networks, and
* allocate and put two lines into line's list */
triangle_list_insert(tri, triangles);
l = line_new_and_assign(line->a, c);
if (line_list_query(l, lines) < 0) {
line_list_insert(l, lines);
for (z = list->tail->next; z != list->head; z = z->next) {
net = (delaunay_net_t *)z->object;
o = delaunay_net_get_coordinate(net);
if (point_cmp(o, line->a) == 0) {
delaunay_net_connect(netc, net);
break;
}
}
} else line_destroy(l);

l = line_new_and_assign(c, line->b);
if (line_list_query(l, lines) < 0) {
line_list_insert(l, lines);
for (z = list->tail->next; z != list->head; z = z->next) {
net = (delaunay_net_t *)z->object;
o = delaunay_net_get_coordinate(net);
if (point_cmp(o, line->b) == 0) {
delaunay_net_connect(netc, net);
break;
}
}
} else line_destroy(l);
break;
}
}
}

// de-allocation lines and triangles
/*
while (triangles->count) {
tri = triangle_list_pop(triangles);
//triangle_dump(tri);
triangle_destroy(tri);
}
triangle_list_destroy(triangles);
*/
list->head->object = (void *)triangles;

while (lines->count) {
line = line_list_pop(lines);
//line_dump(line);
line_destroy(line);
}
line_list_destroy(lines);

point_destroy(v1);
point_destroy(v2);
point_destroy(v3);
point_destroy(org);

// return list;
}

Its reason is because of four points on the same circle. Delaunay Triangulation is determined by whether there is a point in circumference circle, but for this four points it is impossible to construct the unique delaunay triangulation.

The triangle with red color has properties like this;
three vertex: (333, 311), (333, 183), (308, 315)
center of circumference circle: (309.94, 247)
radius of circumference circle: 68.027688

The triangle with blue color has properties like this;
three vertex: (308, 315), (333, 311), (308, 179)
center of circumference circle: (309.94, 247)
radius of circumference circle: 68.027688

In this case it is necessary that additional algorithm is to avoid a crossing between lines.

The subroutine to detect whether two lines has intersecting point or not. It is referenced at http://paulbourke.net/geometry/lineline2d/

int intersect_point_between_lines(point_t *p, line_t *l, line_t *m)
{
real_t denom, ua, ub;

assert(p);
assert(l);
assert(m);

// Pl = l->a + Ua * (l->b - l->a)
// Pm = m->a + Ub * (m->b - m->a)
// for Pl = Pm
denom = (m->b->y - m->a->y) * (l->b->x - l->a->x) - (m->b->x - m->a->x) * (l->b->y - l->a->y);
ua = (m->b->x - m->a->x) * (l->a->y - m->a->y) - (m->b->y - m->a->y) * (l->a->x - m->a->x);
ub = (l->b->x - l->a->x) * (l->a->y - m->a->y) - (l->b->y - l->a->y) * (l->a->x - m->a->x);

if (abs(denom) < REAL_EPSILON) {
if (abs(ua) < REAL_EPSILON && abs(ub) < REAL_EPSILON)
return 0; // COINCIDENT
return -1; // PARALLEL
}

ua = ua / denom;
ub = ub / denom;

if (ua < 0 || ua > 1 || ub < 0 || ub > 1)
return -2; // Not intersecting

p->x = l->a->x + ua * (l->b->x - l->a->x);
p->y = l->a->y + ua * (l->b->y - l->a->y);

if ((abs(ua) < REAL_EPSILON || abs(ua - 1.0) < REAL_EPSILON) ||
(abs(ub) < REAL_EPSILON || abs(ub - 1.0) < REAL_EPSILON))
return 2; // ENDPOINT

return 1; // INTERSECTING
}

And the code of modified delaunay triangulation is here.

void delaunay_list_create_2d_network(delaunay_list_t *list, point_list_t *points, bitmap_t *roi)

{

dlink_t *x, *y, *z;

point_t *o, *a, *b, *c;

point_t *v1, *v2, *v3, *org;

delaunay_net_t *neta, *netb, *netc, *net;

//dlist_t *triangles;

triangle_list_t *triangles;

triangle_t *tri;

//dlist_t *lines;

line_list_t *lines;

line_t *line, *l, *m;

real_t da, db, dc, s;

real_t val, vmin;

real_t radius, area;

assert(list);

assert(points);

// allocate delaunay network list

//list = delaunay_list_new();

for (x = points->tail->next; x != points->head; x = x->next) {

net = delaunay_net_new((point_t *)x->object);

delaunay_list_insert(net, list);

}

// starting point

neta = delaunay_list_get_net(0, list);

a = delaunay_net_get_coordinate(neta);

// find minimum distance between points

for (netb = NULL, x = list->tail->next->next; x != list->head; x = x->next) {

net = (delaunay_net_t *)x->object;

o = delaunay_net_get_coordinate(net);

val = distance_between_points(a, o);

if (val < REAL_EPSILON) continue;

if (netb == NULL || val < vmin) { netb = net; vmin = val; }

}

b = delaunay_net_get_coordinate(netb);

// find minimum radius of circumscribed circle for triangle

for (netc = NULL, x = list->tail->next->next; x != list->head; x = x->next) {

net = (delaunay_net_t *)x->object;

if (net == netb) continue;

o = delaunay_net_get_coordinate(net);

// Heron's formula to find the area formed by three points

da = distance_between_points(a, b);

db = distance_between_points(b, o);

dc = distance_between_points(o, a);

s = (da + db + dc) / 2;

area = sqrt(s * (s - da) * (s - db) * (s - dc));

if (area < REAL_EPSILON) continue;

// calculate radius of circumscribed circle for triangle

radius = (da * db * dc) / (4 * area);

if (netc == NULL || radius < vmin) { netc = net; vmin = radius; }

}

c = delaunay_net_get_coordinate(netc);

v1 = point_new();

v2 = point_new();

v3 = point_new();

org = point_new();

// re-ordering

point_subtract(v1, b, a);

point_subtract(v2, c, a);

point_xproduct(v3, v1, v2);

if (v3->z > 0) { // counterclockwise following a to b to c ?

net = netb, netb = netc, netc = net;

o = b, b = c, c = o;

}

// connect delaunay network

delaunay_net_connect(neta, netb);

delaunay_net_connect(netb, netc);

delaunay_net_connect(netc, neta);

// Synthesize a triangle and insert it.

//triangles = dlist_new();

triangles = triangle_list_new();

tri = triangle_new_and_assign(a, b, c);

triangle_list_insert(tri, triangles);

//triangle_dump(tri);

// Synthesize three line and insert it.

lines = line_list_new();

line = line_new_and_assign(a, b);

line_list_insert(line, lines);

line = line_new_and_assign(b, c);

line_list_insert(line, lines);

line = line_new_and_assign(c, a);

line_list_insert(line, lines);

// main routine

for (x = lines->tail->next; x != lines->head; x = x->next) {

line = (line_t *)x->object;

//line_dump(line);

for (y = list->tail->next; y != list->head; y = y->next) {

netc = (delaunay_net_t *)y->object;

c = delaunay_net_get_coordinate(netc);

if (point_cmp(c, line->a) == 0 || point_cmp(c, line->b) == 0)

continue;

// counter-clockwise rotation of the triangle formed by a, b, and c

point_subtract(v1, line->b, line->a);

point_subtract(v2, c, line->a);

point_xproduct(v3, v1, v2);

if (v3->z < REAL_EPSILON) continue;

// Does any line crossing with line in the line list ?

l = line_new_and_assign(line->a, c);

m = line_new_and_assign(c, line->b);

for (z = lines->tail->next; z != lines->head; z = z->next) {

// Is there an intersecting ?

if ((intersect_point_between_lines(v3, l, (line_t *)z->object) == 1) ||

(intersect_point_between_lines(v3, m, (line_t *)z->object) == 1)) {

line_destroy(l); l = NULL;

line_destroy(m); m = NULL;

break;

}

if (l == NULL && m == NULL) continue;

// Does the triangle exist in the pool of delaunay triangles before ?

tri = triangle_new_and_assign(line->a, c, line->b);

//if (is_triangle_in_list(tri, triangles)) {

if (triangle_list_query(tri, triangles) >= 0) {

triangle_destroy(tri);

tri = NULL;

continue;

}

// Does the unique triangle form the delaunay triangle ?

radius = triangle_circumcircle(org, tri);

for (z = list->tail->next; z != list->head; z = z->next) {

net = (delaunay_net_t *)z->object;

o = delaunay_net_get_coordinate(net);

if (point_cmp(o, line->a) == 0 ||

point_cmp(o, line->b) == 0 ||

point_cmp(o, c) == 0)

continue;

if (distance_between_points(o, org) < radius) {

triangle_destroy(tri);

tri = NULL;

break;

}

if (tri) {

/* put the trinangle into triangle's list,

* establish connection between delaunay networks, and

* allocate and put two lines into line's list */

triangle_list_insert(tri, triangles);

//triangle_dump(tri);

//point_dump(org);

//printf("radius: %lf\n", radius);

//l = line_new_and_assign(line->a, c);

if (line_list_query(l, lines) < 0) {

line_list_insert(l, lines);

//line_dump(l);

for (z = list->tail->next; z != list->head; z = z->next) {

net = (delaunay_net_t *)z->object;

o = delaunay_net_get_coordinate(net);

if (point_cmp(o, line->a) == 0) {

delaunay_net_connect(netc, net);

break;

}

} else line_destroy(l);

//m = line_new_and_assign(c, line->b);

if (line_list_query(m, lines) < 0) {

line_list_insert(m, lines);

//line_dump(m);

for (z = list->tail->next; z != list->head; z = z->next) {

net = (delaunay_net_t *)z->object;

o = delaunay_net_get_coordinate(net);

if (point_cmp(o, line->b) == 0) {

delaunay_net_connect(netc, net);

break;

}

} else line_destroy(m);

//keyhit();

break;

}

// de-allocation lines and triangles

while (triangles->count) {

tri = triangle_list_pop(triangles);

//triangle_dump(tri);

triangle_destroy(tri);

}

triangle_list_destroy(triangles);

list->head->object = (void *)triangles;

while (lines->count) {

line = line_list_pop(lines);

//line_dump(line);

line_destroy(line);

}

line_list_destroy(lines);

point_destroy(v1);

point_destroy(v2);

point_destroy(v3);

point_destroy(org);

// return list;

}

At result the line crossing is disappeared.

Unfortunately, two delaunay triangles are deleted in the figure. In the figure dot is the mean center of delaunay triangle. As you see, you can find two triangles without mean center. It is remaining side effect of four points on the same circle.

So, it needs to change comparison between variables expressed in floating points.

Change
....

// Does the unique triangle form the delaunay triangle ?

radius = triangle_circumcircle(org, tri);

for (z = list->tail->next; z != list->head; z = z->next) {

net = (delaunay_net_t *)z->object;

o = delaunay_net_get_coordinate(net);

if (point_cmp(o, line->a) == 0 ||

point_cmp(o, line->b) == 0 ||

point_cmp(o, c) == 0)

continue;

if (distance_between_points(o, org) < radius) {

triangle_destroy(tri);

tri = NULL;

break;

}

}
....
to
....

// Does the unique triangle form the delaunay triangle ?

radius = triangle_circumcircle(org, tri);

for (z = list->tail->next; z != list->head; z = z->next) {

net = (delaunay_net_t *)z->object;

o = delaunay_net_get_coordinate(net);

if (point_cmp(o, line->a) == 0 ||

point_cmp(o, line->b) == 0 ||

point_cmp(o, c) == 0)

continue;

val = distance_between_points(o, org);

if ((val < radius) && (round(100*val) < round(100*radius)) {

triangle_destroy(tri);

tri = NULL;

break;

}

}
....

Then,

Watershed by immersion (Vincent-Soille watershed algorithm)

2011-09-18T03:44:00.000-07:00

procedure: Watershed-by-Immersion
Input: digital grey scale image G = (D;E; im).
Output: labelled watershed image lab on D.

#define INIT   1 // initial value of lab image
#define MASK   2 // initial value at each level
#define WSHED 0 // label of the watershed pixels
#define FICTITIOUS (-1, -1) // fictitious pixel 2 D

curlab <= 0; // curlab is the current label
fifo_init(queue);

for all p in D do
   lab[p] <= INIT;
   dist[p] <= 0; // dist is a work image of distances
end for

SORT pixels in increasing order of grey values (minimum hmin, maximum hmax)

// Start Flooding
for h = hmin to hmax do // Geodesic SKIZ of level h   1 inside level h

   for all p in D with im[p] = h do // mask all pixels at level h
             // these are directly accessible because of the sorting step
      lab[p] <= MASK;
      if p has a neighbour q with (lab[q] > 0 or lab[q] == WSHED) then
         // Initialize queue with neighbours at level h of current basins or watersheds
         dist[p] <= 1;
         fifo_add(p, queue);
      end if
   end for

   curdist <= 1;
   fifo_add(FICTITIOUS, queue);

   loop // extend basins
      p <= fifo_remove(queue);
      if p == FICTITIOUS then
         if fifo_empty(queue) then
            break;
         else
            fifo_add(FICTITIOUS, queue);
            curdist <= curdist + 1;
            p <= fifo_remove(queue);
         end if
      end if
      for all q in neighbours of p do // labelling p by inspecting neighbours
         if dist[q] < curdist and (lab[q] > 0 or lab[q] == WSHED) then
            // q belongs to an existing basin or to watersheds
            if lab[q] > 0 then
               if lab[p] == MASK or lab[p] == WSHED then
                  lab[p] <= lab[q];
               else if lab[p] != lab[q] then
                  lab[p] <= WSHED;
               end if
            else if lab[p] == MASK then
               lab[p] <= WSHED;
            end if
         else if lab[q] == MASK and dist[q] == 0 then // q is plateau pixel
            dist[q] <= curdist + 1;
            fifo_add(q, queue);
         end if
      end for
   end loop

   // detect and process new minima at level h
   for all p in D with im[p] = h do
      dist[p] <= 0; // reset distance to zero
      if lab[p] == MASK then // p is inside a new minimum
         curlab <= curlab + 1; // create new label
         fifo_add(p, queue);
         lab[p] <= curlab;
         while not fifo_empty(queue) do
            q <= fifo_remove(queue);
            for all r in neighbours of q do // inspect neighbours of q
               if lab[r] == MASK then
                  fifo_add(r, queue);
                  lab[r] <= curlab;
               end if
            end for
         end while
      end if
   end for
end for
// End Flooding

Reference

The Watershed Transform: De nitions, Algorithms and
Parallelization Strategies
Jos B.T.M. Roerdink and Arnold Meijster
Institute for Mathematics and Computing Science
University of Groningen
P.O. Box 800, 9700 AV Groningen, The Netherlands
Email: roe@cs.rug.nl,a.meijster@rc.rug.nl

Debian Wheezy의 safe-upgrade 해제한 최신 업데이트

2011-09-14T19:57:00.000-07:00

다른건 아무런 문제가 없지만 gnome-terminal과 xterm등등의 gnome 어플리케이션의 동작에 문제가 생겼다. 아무래도 gnome버전이 2.x에서 3.0.1로 옮아가는 가정에서 생기는 문제점인거 같다. 일명 slow key라고 하는 한 5번 누르면 하나 입력이 되는 키보드 먹통현상이 가장 큰 문제이다. 그리고 gnome의 기본 테마가 적용이 안된다. xorg에서 제공하는 기본 테마가 적용되서 다른 어플리케이션과 스크롤바와 메뉴등의 테마가 다르다. 수정되겠지 하면서 쓰고 있다. 그냥 safe-upgrade로 업테이트 해야 겠다. 빨리 gnome 3.x가 안정화 되길 바란다.

여러가지 씨름을 한후 ibus와 scim의 충돌로 인한 키보드 입력에 문제가 발생한 것을 알아 차리고 scim을 지웠다. 그 이후로는 ibus를 입력기로 사용하고 있다. 키보드 입력문제는 해결을 보았으나, 여전히 테마적용문제는 풀리지 않고 있다.

이 문제도 해결을 보았다. home 디렉토리 밑의 .gconf를 싹 지웠더니 설정이 제대로 되고 있다. 하지만 그동안의 설정이 사라졌다.

멋있는 방문자 지도 제공 사이트

2011-09-12T21:35:00.000-07:00

http://www.revolvermaps.com/?target=setup

지도와 해상도 기타등등을 선택하고 코드를 자신의 사이트에 추가하면 된다. 간단하다.

이러한 노력의 결과물을 많이 꽁짜로 제공하는데 왜 그런지 의문스럽다.

KS 규격 나사산 그리기(탭 깊이, 나사 맞춤 깊이)

2011-09-10T21:00:00.000-07:00

* M3

숫나사의 바깥 지름 d: 3mm
나사를 내기 위한 구멍의 지름 d1: 2.4mm
드릴 깊이 b: 6(강, 주강, 청동, 청동주물), 7.5(주철), 8.5(알미늄과 기타의 경합금제)
탭 깊이 a: 3(강, 주강, 청동, 청동주물), 4.5(주철), 5.5(알미늄가 기타의 경합금제)

* M6

숫나사의 바깥 지름 d: 6mm
나사를 내기 위한 구멍의 지름 d1: 5mm
드릴 깊이 b: 10(강, 주강, 청동, 청동주물), 13(주철), 15(알미늄과 기타의 경합금제)
탭 깊이 a: 6(강, 주강, 청동, 청동주물), 9(주철), 11(알미늄가 기타의 경합금제)

Untouchable HDD noise on Laptop as having age.

2011-09-05T22:02:00.000-07:00

I give a priority to change HDD as SDD for out of noise.

But my debian system is too heavy to reconfigure.

One of the disk copy between HDD is through "cat" commnand.

# cat /dev/sda1 > /dev/sda2

or

# dd if=/dev/sda1 of=backup , dd if=backup of=/dev/sda2

Other is backup image of HDD through "tar".

# cd /
# tar cvpjf backup.tar.bz2 --exclude=/proc --exclude=/lost+found --exclude=/backup.tar.bz2 --exclude=/mnt --exclude=/media --exclude=/sys /

Then copy this tar into new HDD and extract it.

# cd $(new_HDD)
# tar xvpjf backup.tar.bz2 -C /
# mkdir proc
# mkdir lost+found
# mkdir mnt
# mkdir sys
# mkdir media

3x3 bitmask list sorted by crossing number and pixel count

2011-09-05T02:34:00.000-07:00

Here is the source code for generating 3x3 bitmask list.

#include

int value_3x3[256][9+2];

int main(int argc, char *argv[])
{
int i, j, value;
int crossing = 0, count = 0;

for (i = 0; i < 256; i++) {
value_3x3[i][0] = '*';
for (j = 0; j < 8; j++) {
if (i & (1 << j)) value_3x3[i][j+1] = '*';
else value_3x3[i][j+1] = '=';
}
crossing = 0;
count = 0;
value = value_3x3[i][8];
for (j = 1; j <= 8; j++) {
if (value_3x3[i][j] != value) crossing++;
if (value_3x3[i][j] == value_3x3[i][0]) count++;
value = value_3x3[i][j];
}
value_3x3[i][9] = crossing;
value_3x3[i][10] = count;

}
for (crossing = 0; crossing <= 8; crossing += 2) {
for (count = 0; count <= 8; count++) {
for (i = 0; i < 256; i++) {
if (value_3x3[i][9] == crossing && value_3x3[i][10] == count) {
printf("%c%c%c\n", value_3x3[i][4], value_3x3[i][3], value_3x3[i][2]);
printf("%c%c%c %d %d\n", value_3x3[i][5], value_3x3[i][0], value_3x3[i][1], value_3x3[i][9], value_3x3[i][10]);
printf("%c%c%c\n", value_3x3[i][6], value_3x3[i][7], value_3x3[i][8]);
printf("\n");
}
}
}
}
return 0;
}

And then the sorted list is here.

3x3 bitmask, crossing number, pixels equal to center
-------------------------------------------------------------------------

To remove clicking configuration of hdparm in Asus 1215t debian system

2011-09-01T22:12:00.000-07:00

1. edit /etc/hdparm.conf

...

/dev/sda
{
apm = 255
spindown_time = 0
dma = on
}
...

2. edit /etc/apm/event.d/20hdparm

...
APMD_DRIVES=sda
...

power_performance ()
{
# Disable IDE hard disk spindown.
for DRIVE in $APMD_DRIVES; do
"${HDPARM}" -q -S 0 -B 254 "${DRIVE}" || true
done
}
...

3. edit /etc/acpi/fglrx-powermode.sh

...
for dev in /dev/sd? /dev/hd? ; do
if [ -b $dev] ; then

# Check for APM support; discard errors since not all drives
# support HDIO_GET_IDENTITY (-i).
if /sbin/hdparm -i $dev 2> /dev/null | grep -q 'AdvancedPM=yes' ; then
if [ -f /sys/class/power_supply/AC0/online ] && [ grep -q '1' /sys/class/power_supply/AC0/online ]; then
/sbin/hdparm -B 254 $dev
else
/sbin/hdparm -B 128 $dev
fi
fi
fi
done
...

[번역] hdparm - 완벽 가이드 (hdparm - complete guide)

2011-09-01T21:41:00.000-07:00

출처 : http://forums.debian.net/viewtopic.php?t=32794

이글의 목적은 당신의 하드 디스크의 load cycle(하드 디스크의
clicking 음이 반복되어 나타나게 되는 원인)을 조절하기 위한 매개변수를 설정하는 것이다. 이를 통해서 laptop이나 desktop에서 불필요한 하드의 헤드의 정차(하드 디스크의 끍는 소리는 이 헤드의 정차음이다)를 멈출 수 있다. 그러나 이렇게 설정하더라도 베터리로 동작할때 laptop이 헤드의 정차를 하도록 한다.

그레서 나는 사람들이 작성한 스크립트를 이용한 방법을 전혀 사용하지 않았다. 여러분은 이 화면에 있는 것들을 복사해서 사용해야 할 것이다. 대신 약 1년쯤 전에 여러분이 지지하는 방식의 방법을 찾았다. 나는 항상 모두 공유하고 싶지만 워낙에 게을러서. 지금은 작업을 밀어붙이고 있으며 진작에 이것을 공유하지 않은 것에 대해서 자책하고 있다. 나는 여러분이 이 방법을 온라인상에서 찾아 볼 수 있을지 어떨지는 잘 모르겠다. howto 문서와 떨어져 있기 때문에 그러나 여러분은 검색을 통해서 찾아 볼 수 있다.

그리고 문서 작성하면서 기억해 낸 것인데 한가지 주의할 점은 255의 apm 순위는 어떤 컴퓨터에서는 작동하지 않는다. 여러분이 이 값으로 254를 사용하지 않는다면 모든 것은 기본적으로 동일하다.

1. 여러분은 repos(패키지 저장소)로 부터 구할 수 있는 hdparm, sysvconfg, smartmontools를 설치해야한다.

2. hdparm 이 로그인시 시작되는지 확인하자. sysvconfg를 실행하고 hdparm을 enable시키면 로그인시 자동으로 시작된다.

3. 그럼 여기서 내가 다르게 접근한 약간의 것들이 필요하다. 여러분이 직접 수정해야할 몇개의 파일들이 있다. 로그인시 로드되거나 이미 존재하는 hdparm의 config 파일과 약간의 스크립트가 필요하다. 그래서 여러분은 정확한 정보를 그들에게 알려주어야한다.

3.1. /etc/hdparm.conf 수정

...

/dev/${your hard drive name}
{
apm = 254
spindown_time = 0
dma = on
}
...

3.2. 기타 파일 작성

# cd /etc/apm/event.d
# ls
# nano 20hdparm

그저 간단하다. 다음 몇 줄을 추가하자.

APMD_DRIVES=sda ##Again put YOUR device

그런 다음 스크립트가 동작하도록 편집하자.

power_performance ()
{
# Disable IDE hard disk spindown.
for DRIVE in $APMD_DRIVES; do
"${HDPARM}" -q -S 0 -B 255 "${DRIVE}" || true
done
}

기본적으로 여러분은 AC 전원일때 loading/unloading cycles을 얻을 수 없다.

3.3. 마지막 단계.

ACPI 전원 관리를 위한 것이다.

# cd /etc/acpi/ac.d
# ls
# nano 90-hdparm.sh ## or whatever ls shows you

...

for dev in /dev/sd? /dev/hd? ; do
if [ -b $dev ] ; then
# Check for APM support; discard errors since not all drives
# support HDIO_GET_IDENTITY (-i).
if hdparm -i $dev 2> /dev/null | grep -q 'AdvancedPM=yes' ; then
if [ $STATE = "BATTERY" ] ; then
hdparm -B 128 $dev
else
hdparm -B 255 $dev
fi
fi
fi
done
fi

...

3.4. reboot후 확인하기

# smartctl -a /dev/sda | grep Load_Cycle_Count

# smartctl -a -d ata -i /dev/sda | grep Load_Cycle_Count

이 명령은 94116이라는 다섯 자리의 숫자를 제공할 것이다. 몇 분후에 다시 해보자. 숫자는 같아야 한다.

그리고 몇일 후에 알아낸 것인데 내 컴퓨터의 경우 hdparm이 suspend후 resume되었을대 작동되지 않는다. 이경우 나는 다음과 같이 수동으로 사용한다.

# hdparm -B 254 -S 0 /dev/sda

Square Root of Complex Number

2011-08-30T22:48:00.000-07:00

As before, r = sqrt(a * a + b * b)

sqrt(a + J*b)
= sqrt(r) * (sqrt((1 + a / r) / 2) (+ or -) J * sqrt((1 - a / r) / 2))
= sqrt((r + a) / 2) (+ or -) J * sqrt((r - a) / 2)
= sqrt((sqrt(a * a + b * b) + a) / 2) (+ or -) J * sqrt((sqrt(a * a + b * b) - a) / 2)

caution ! (+ or -) should be chosen to be the same as the sign of b

example:

sqrt(3 - J * 2) = sqrt((sqrt(3 * 3 + 2 *2) + 3) / 2) - J * sqrt((sqrt(3 * 3 + 2 * 2) - 3) / 2)

정부, 삼성, LG의 독자 OS 개발에 대한 생각

2011-08-23T00:14:00.000-07:00

아마 예전에도 있었던 행태가 아닌가 기억한다. 보도에서 600억이라는 금액을 투자한다고 하는데 우선 드는 생각이 연구비 사냥꾼에 좋은 먹이감이 아닌가 생각이 든다. 누군가의 펀딩 욕심이거나 누군가의 실적올리기용 꼼수이 아닌가 한다.

아이폰에서 안드로이드폰으로 이어지는 OS 와 제조사의 결합이 이루어 졌다. 뭐 그렇다고 안드로이드가 개방정책을 철회하는것은 아니지만 (혹 그럴지도 모른다) 이 결합에 많이들 걱정하는 것 같다. 뭐 사용자 입장에서 나쁜 것 보다 좋은 점이 더 많다. 우후 죽순처럼 양산되는 완벽하지 못한 안드로이드 폰들을 사용하는 것보다 잘 정돈된 안드로이드폰이 훨씬 쓸만 하니깐. 그런 의미에서 삼성과 LG는 안드로이드 OS자체를 자신의 용도에 맞게 고쳐 쓸 수 있는 혹은 수정할 수 있는 인력의 확보가 더 중요할 지도 모른다. 한 발 더 나아가서 안드로이드 OS에서 파생되는 자체 OS를 개발 할 수도 있다. 안드로이드의 OS의 라이센스를 잘 모르겠지만 물리적으로 얼마든지 가능하다.

그렇지 않고 이 두 OS와 경쟁하기 위해서 지금 당장 처음부터 개발한다는 것은 무모하다. 일정 수준의 완성도를 갖기 위해서 두 회사가 거쳐온 기간과 자금, 지금 현재 유지, 보수되고 있는 노력들을 살펴보면 이해가 쉽다. 만약 안드로이드의 라이센스정책이 그러한 파생 OS를 허용하지 않거나 공개를 철회할 경우, 비슷한 부류의 GNU프로젝트에 투자를 한다던가 개발자를 고용해서 그 프로젝트에 기여하는 방식을 생각해 볼 수 있다. 물론 GPL 라이센스 형식으로 말이다. 그게 유지, 보수 및 개발에 훨씬 유리할 지도 모른다. 안드로이드의 태생도 이러한 방식의 소프트웨어를 기반으로 탄생한 것이었으니깐 말이다. 물론 사적이익 창출이 목표인 삼성이나 LG에서 못마땅 할 수 있지만 공공의 이익을 생각하는 정부 입장에서는 충분히 고려해 볼만한 가치가 있다. 반면 삼성이나 LG의 입장에서는 노키아에 팽당한 인텔의 MiGO프로젝트에 참여하는 방식도 충분히 고려해 볼만한 사항이다. RIM이나 HP를 염두해 두는 목소리도 많은데 내가 보기엔 이 둘과 손을 잡으면 지금 현재 상태의 구글과의 관계보다 더 어려울 수 있다. 이미 OS와 제조를 다 가지고 있기 때문이다.

사용자들은 아무런 걱정을 할 필요가 없다 좋은 것만 잘 골라서 쓰기만 하면된다. 개발자들도 마찬가지로 편한 환경을 골라서 개발하면 아무 문제 없다. 그저 이익이 줄어들어서 걱정하는 삼성과 LG가 문제이지. 그리고 안드로이드가 이상해 지더라도 걱정할 것 없다. 비슷한 목표를 가진 GNU프로젝트는 얼마 든지 많으니깐 당장 쓸정도가 아닌지 모르겠지만 스마트폰의 명세만 잘 공개 된다면 그 누군가의 그룹들이 뚝딱뚝딱하면서 쓸모 있게 만들테니깐 자원들은 풍부하다. 안드로이드도 많은 부분 이러한 자원들에 의존하고 있다. 그러니깐 대부분의 사람들은 삼성, LG 걱정할 필요없다.

OpenGL 프로그래밍에 대해서

2011-08-15T08:13:00.000-07:00

OpenGL은 하드웨어 독립적인 2D/3D 그래픽 라이브러리이다. SGI사의 IRIS GL라이브러리를 기초로 개발되었지만 개방형이다. 그래픽 라이브러리이기 때문에 기본적인 창관리와 입력처리에 대한 부분이 없다. 그래서 GLUT이라는 보조 도구를 이용해서 창을 관리하고 입력을 관리하는 방법을 제공한다. 이것도 역시 시스템 독립적이기 때문에 OpenGL과 GLUT의 조합이면 어떤 시스템에서도 동일한 결과를 갖는 프로그램을 개발 할 수 있다.

GLUT: OpenGL 유틸리티 툴킷을 줄인 말인 것 같다. 내부에 많은 상태 변수들을 관리하고 있단다. 간단히 초기화하고 창에 점과 선만 그릴 수 있었음 만족한다. 간단히 초기화 하는 법은 다음과 같다.

#include

void MyIdle(void)
{
/* 사용된 변수값을 다음 프레임의 상태로 수정해 주는부분*/
};

void MyDisplay(void)
{
/* 프레임에 실제 그림을 그리는 OpenGL명령어들이 들어간다. */
/* 프레임이 다 그려진 다음 버퍼를 교체한다. */
glutSwapBuffers();
};

void main(int argcp, char **argv)
{
/* GLUT상태를 초기화 한다. */
glutInit(&argcp, argv);
glutInitWindowSize(640, 480);
glutInitWindowPosition(0, 0);

/* 창을 연다. */
glutCreateWindow("My OpenGL Application");

/* 디스플레이모드를 선택한다.:
이중버퍼 & RGBA색상 */
glutInitDisplayMode(GLUT_RGBA | GLUT_DOUBLE);

/* 콜백함수들을 등록한다. */
glutDisplayFunc(MyDisplay)
glutIdleFunc(MyIdle)

/* 이벤트 처리엔진을 시작한다 */
glutMainLoop();
};

카메라 신호 분석하기

2011-08-07T09:27:00.000-07:00

카메라 센서 설정을 보게 되면 Active High혹은 Active Low라는 설정이 있다. 회사마다 다르겠지만, 아래의 두가지 경우를 설정하는 방식을 다양하게 제공하고 있을 꺼다.

설명하자면 바로 아래의 사진은 Active Low로 설정된 HREF신호이다. 초기 상태는 이와 같은 상태 일 것이다. 대부분이. HREF신호가 Rising edge로의 변경은 PCLK의 Low edge에서 발생한다. 즉 PCLK 의 rising edge에서 화소값을 읽어 들이면 된다.

이와는 반대로 Active High는 정반이다. 마찬가지로 HREF신호의 HIGH상태로의 전이는 PCLK의 Rising Edge에서 잃어난다. 그럼 동일한 방식으로 화소값은 Falling edge에서 읽어 들이면 된다.

문제는 무엇이냐 하면 둘 중 어떤 것으로 설정해도 화면은 나온다. 읽어 들이는 시점이 이전 상태를 읽어 들여서 문제지. 이상하게 좌측에 검은 줄이 생긴다면 의심해 볼 필요가 있다. Active Low 인지 High인지 읽어 들이는 시점은 정확한지.

비오는 날

2011-08-07T09:15:00.000-07:00

비가 와서 산책은 하고 싶고 산에는 못가도 언덕 비슷한 것은 가고 싶다.

비가 많이 와도 나무가 어느 정도 막아 주었다.

이 곳을 내려가면 산에 바짝 붙은 아담한 집이 있다. 보기엔 좋은데 왕래하는 사람이 많아져서 실상 좋은 집은 아닐꺼다.

밤깍는 기계

2011-08-07T09:10:00.000-07:00

농기계 전시회에서 구경한 것중 가장 유용한 것 아니었나 싶다. 간단하면서 가장 유용한 것이 진리다.

우기의 자전거 나들이

2011-08-07T08:50:00.000-07:00

장마가 너무나도 길었다. 잠시 소강상태 일때 간만에 열대야가 시작되었을때 방구석에 있는 녀석을 데리고 나갔다. 약간 질퍽하겠거니하고 생각했다. 너무 일찍 나왔다 갯벌 수준이었다. 오자마자 분해 할 수 있는 부분은 분해하고 휠먼저 딱았다.

그리고 사무실에 있는 크랭크 분해 공구들을 빌려서 휠 다음으로 더러워진 크랭크를 딱아 주었다.

지금도 비가 온다. 아무래도 여름내내 이럴거 같다. 느긋하게 생각하고 기다려야지.

메모리 이식 ㅋㅋ

2011-08-07T08:43:00.000-07:00

처음엔 메모리가 필요없다고 느꼈다. 심지어는 USB 전용 칩같은건 필요없다고 생각해서 MCU에 딸려 있는 USB를 그냥 사용할 목적으로 사용했었다. 모든 게 잘못된 생각이었다. 결국은 MCU가 CPLD로 바뀌고 USB전용 칩을 달아야만 했다.

하지만 그래도 여전히 부족했다. 그래서 추가적으로 SRAM을 달고 CPLD에 코딩이 추가되었다. 결국은 동작하였지만, 여전히 불안하다. 도망가고 싶을 정도로 걱정 스럽다. 그래도 시간이 지나 니깐 해보아야 겠구나 하는 곳들이 생겨난다. 그냥 들여다 보고 스트레스 받지 말고 멀리서 보자. 이젠 하드웨어는 해볼데로 해보았으니 소프트웨어를 점검해 봐야 겠다.

5각형으로 만든 돔 (Geodesic dome)

2011-07-29T06:07:00.000-07:00

조명 장치와 카메라를 여러대 달려고 구상중 아크릴판을 통짜로 가공하는 방식은 너무나도 번거롭고 모양의 변형이 많다.

그래서 구조적으로 견고하고 일정한 반경의 돔을 간단한 기초 작업으로 쉽게 만드는 방법을 찾던 중 pentagon(5각형)으로 암실 돔을 구성하면 어떨까? 하고 찾아 보았다. 외국에서는 과학 캠프나 실습차원에서 많이 하는 것 같다. 나도 예전에 서울 대공원에서 본 것 같은 기억이 있다. 가본지 너무 오래 되어서 아마도 거기 입장료가 좀 쎘었던거 같다.

아무튼 http://hilaroad.com/camp/projects/dome/domemath.html 에 가보면 원하는 직경의 Geodesic 돔을 구성할 수 있는 계산을 할 수 있게 해준다. 다음에는 이것을 번역해 보아야 겠다. 위와 같은 것을 레고 처럼 구성할 수 있다면 조명과 카메라를 자유 자재로 안정적인 위치에 고정 시킬 수 있을 거 같다.

Hila Science Camp

Dome Math

Geodesic Domes

돔을 만드는 사람들은 이러한 방식의 돔을 2V Icosa Alternate 라고 부른다. icosahedron의 각 삼각형 표면안에 딱 맞는 4개의 삼각형에 의해서 만들어 진다. icosahedron란 뭔가? icosahedron는 고대 그리스인들에 의해서 만들어진 5개의 구조물(solids) 중의 하나이다. 이것이 그들 5개이다:

어떻게 이름 지어졌을까? 우측 맨 끝이 icosahedron이다. 이것은 20의 삼각형 면을 가지고 있다.

여러분이 정말로 정말로 알고 싶다면, 다음의 링크를 따라가서 어떻게 삼각형이 icosahedron의 각각의 삼각형에 맞아 떨어지는지 알아 볼 수 있다. - Icosahedron Domes Explained!

맨 위에서 보여주는 것과 같은 돔을 만들기 위해서, 여러분은 40개의 삼각형, 30개의 isosceles 삼각형, 그리고 10 equilateral 삼각형을 만들어야 한다.

돔을 디자인해 봅시다!

1단계:

여러분이 만들고자하는 돔의 직경을 결정하자. 우선 작은 30cm정도의 직경을 갖는 작은 모델을 만들어 보자.

2단계:

필요한 삼각형의 크기를 결정하자. 우리가 만들어야할 삼각형은 각면에 접착용 탭이 필요하다. 다음과 같이 생겼다:

10개의equilateral 삼각형은 각 면이 "A" 길이로 같다. "A"를 결정하기 위해서 "Chord Factor" 라고 불리는 숫자를 사용한다.

"A"를 위한 chord factor는 0.61803. "A"의 길이를 결정하기 위해서 여러분이 디자인한 돔의 반경에 이 chord factor를 곱하자. 반경은 직격의 1/2에 해당한다는 것을 잊지 말자. 지금 디자인하는 돔은 직경이 30cm이고 반경은 15cm이다.

A = 돔의 반경 x 0.61803
A = 15 cm x 0.61803
A = 9.27045 cm
A=9.27 cm

30 isosceles 삼각형은 한 면이 "A"이고, 두 면은 "B"이다. "A"값은 9.27cm으로 알려져 있고, "B"를 계산하기 위해서 우리는 다른 chord factor를 사용한다. "B"를 위한 chord factor는 0.54653이다.

B = 돔의 반경 x 0.54653
B = 15 cm x 0.54653
B = 8.19795 cm
B = 8.20 cm

3단계:

이 돔을 구성하기 위해서 "A"와 "B"값으로 구성된 삼각형 요소들을 만들어야 한다. 단단한 cardboard위에 30cm 직경의 원을 그리고 다음 링크의 내용을 따라하도록 하자. instructions for constructing the 50 cm dome.어떤 크기의 돔이라도 만들 수 있다. 그저 반경에 chord factor만 곱하면 된다. 그리고 나서 조립하는 순서만 따라하면 된다.

이 페이지는 Hila Science Camp: http://hilaroad.com 에 의해서 만들어 졌습니다.

출처 : http://hilaroad.com/camp/projects/dome/domemath.html

돔 조립하기

Geodesic dome들은 우아하고 정교한 구조물이다. 가벼운 무게와 강인함으로 인해서 재료를 아주 효율적으로 사용하도록 만든다. 북극과 같은 먼 지역에서 사용하기에 이상적인 구조다. 언젠가 달이나 covering(?) 도시에서 돔 구조물을 볼지도 모른다.

이 사이트는 종이로 geodesic dome을 만드는 과정을 차례차례로 제공한다. 그 내용은 다른 크기의 dome을 디자인할 수 있는 것도 포함한다.

이 프로젝트는 많은 지방과 주에서의 과학과 기술 교과 과정의 구조와 매커니즘 부분을 지원한다. 교과 과정의 목표는 geodesic dome이나 shell 구조물을 구축하는 것을 종종 포함한다. geodesic dome은 독특한 구조에서 생기는 응측력과 팽창력을 보여준다. 여러분은 다양한 재질로된 geodesic dome을 구축할 수 있다. 지금의 dome은 종이로 구축될 것이다. 다른 크기의 dome을 구축하기 위해선 여기서 제시된 디자인 개념을 사용하고 Bristol 보드나 골판지, plywood 심지어는 얼음과 같은 재료로도 도전해 보자!

6, 7 학년의 학생에 대해서 어떤 교사들은 활기차고 화려한 조작으로 변화 시켜서 이 같은 활동을 미술 시간과 결합하기도 한다.

이 프로젝트의 수학적인 부분은 8학년과 9학년의 응용 수학에 아주 적당하다. 응용 수학과 "구조와 매커니즘"에 대한 여러분의 주와 지방의 교과 과정 문서를 확인해보자.

Geodesic Dome 구축하기

Geodesic dome은 Buckminister Fuller에 의해서 대중화 되고 삼각형들의 구조적 접목이 dome이나 shell 구조물의 고유의 강도에 더해 졌다.

다음은 50cm의 직경을 갖는 geodesic dome의 구축하는데 필요한 지시사항들이다. 여러분이 다른 크기의 dome을 만들고자 한다면 다음 링크를 따라하자. - DOME MATH

먼저 아래에 보여진 삼각형들의 정확한 그림을 생성해서 시작해보자. 여러분은 30개의 B-A-B 삼각형과 10개의 A-A-A 삼각형이 필요하다. 이 삼각형으로 만들어진 유닛들은 기본 구성 요소가 된다. 여러분은 이 삼각형들을 정확하기 만들어야 하기 때문에 제시된 수치들은 매우 중요하다. 콤파스와 자를 사용하자. A, B의 수치는 삼각형의 각 면의 길이가 되고 접착면을 위한 tab은 여분이 됨을 주의하자.

A-A-A 삼각혀은 equilateral 이라고 하고, B-A-B 삼각형은 isosceles라고 한다.

이것은 기하학에서 훌룡한 예제가 된다. 각 삼각형의 요소요소의 정확한 도면을 만들기 위해서 콤파스와 자, 연필로 작업하자.

각 삼각형의 도면을 그린 후 photo copy만으로 dome을 만드는데 충분하다. 어떤 종류의 종이도 적당하다.
두 가지 색상의 종이를 사용하자. 하나는 A-A-A 삼각형에, 다른 하나는 B-A-B 삼각형에 사용하자. 색상을 다르게 하는 것은 구축을 쉽게하고 여러분의 dome을 좀 더 극적이게 보이도록 만든다.

나는 이 삼각형들의 도면을 그리고 싶지 않다. 그져 dome을 만들고 싶다!
좋다. 이미 그려진 요소들을 위한 프린트 가능한 페이지에 대한 링크를 따라가라.

Link to template for A-A-A triangle: (PDF File)

Link to template for B-A-B triangle: (PDF File)

주의 사항: 프린트 한후 "A"와 "B"면이 정확한 길이로 프린트 되었는지 확인하라. 모든 프린트가 동일한 해상도 이지는 않다.

삼각형을 잘라내고 tab 라인을 조심스럽게 접자. 깔끔하게 접는 것이 무엇보다도 중요하다.
다섯개의 A-B-B삼각형을 조합해서 pentagon을 만들자. "B"면을 접착면을 활용해서 접착시키자. 접착용 풀은 이 프로젝트에 아주 적당하다.

pentagon의 외부 면이 삼각형의 "A" 면이 되도록하는 것을 명심하자.

완성된 pentagon은 위의 이미지와 같게 보여야 한다. 가운데가 쏫아 올라와야 한다. 다섯개의 pentagon을 더 만들자.

판지로 75cm 직경의 원은 만들어서 dome을 위한 바닥을 만들자.

그리고 그 바닥의 중앙에 50cm 직경의 원은 그리자.

바닥의 중앙의 25cm의 원을 잘라내자. 이것은 dome의 안쪽을 접근할 수 있도록 해줄것이다.

중요: 원 주위에 chords의 나열을 조심스럽게 그리자. 각 chord의 길이는 "A"이다.

50cm dome을 위해서 A값은 15.45 cm 이다. 원 주위에 10개의 chord가 완벽하게 들어 맞도록 위치시키자.

이 라인들은 각 삼각형이 어디에 위치해야하는 지를 알려준다.

두개의 pentagon 사이에 하나의 A-A-A 삼각형을 붙여서 두개의 pentagon을 붙이자.

또 다른 삼각형은 다음과 같이 추가하자.

이 완성된 부분의 바닥 부분에 접착제를 붙이고, chords의 한 부분에 각 삼각형의 한 면을 위치시키자.

계속해서 pentagon과 삼각형을 추가하고, 접착제로 이들은 고정시키자.

equilateral ( A - A - A ) 삼각형은 맨 꼭대기의 pentagon 사이사이에 붙이자.

매 꼭대기에 마지막 pentagon을 붙이자.

여러분의 dome은 완성되었다. dome에 관한 것을 더 배우고 싶다면 다음의 링크를 따라가자. Dome Math.

이 페이지는 Hila Science Camp: http://hilaroad.com 에 의해서 만들어 졌습니다.

출처 : http://hilaroad.com/camp/projects/dome/dome.html

그리고 좀 더 복잡하고 실제적인 내용은 다음 링크를 따라가 보자. 다른 내용도 유용한 것들이 많다. - http://www.byexample.com/articles/topics/geodesic_domes

Matrox Solios eV-CL + MIL-Lite 9.0 + 2 * VCC-F22V39ACL

2011-07-29T02:18:00.000-07:00

At first, it is necessary to change the configuration from Solios Camera Link Medium to Solios Camera Link Dual Base in Solios under Boards of left tree on Milconfig during unplugged the two camera. and then, reboot system.

In login time, the system detects Solios eV-CL as a new hardware named by Solios eV-CL Dual 0. And we can select the digitizer in Intellicam.

And the system totally detects three different new hardware like this:

In development part, for each camera a digitizer is able to allocated;

MdigAlloc(MilBoard[0].MilSystem, M_DEV0, DCF_NAME, M_DEFAULT, &MilDigitizer0);
MdigAlloc(MilBoard[0].MilSystem, M_DEV1, DCF_NAME, M_DEFAULT, &MilDigitizer1);

To control Matrox Solios Camera Link Com Port 0(COM3) and Port 1 (COM4) configuring via MsysControl is better than via Windows API.

MsysControl(MilBoard[0].MilSystem, M_UART_DATA_LENGTH+M_UART_NB(M_DEV0), 8);
MsysControl(MilBoard[0].MilSystem, M_UART_PARITY+M_UART_NB(M_DEV0), M_DISABLE);
MsysControl(MilBoard[0].MilSystem, M_UART_STOP_BITS+M_UART_NB(M_DEV0), 1);
MsysControl(MilBoard[0].MilSystem, M_UART_SPEED+M_UART_NB(M_DEV0), 9600);
MsysControl(MilBoard[0].MilSystem, M_UART_TIMEOUT+M_UART_NB(M_DEV0), M_INFINITE); // ms

MsysControl(MilBoard[0].MilSystem, M_UART_DATA_LENGTH+M_UART_NB(M_DEV1), 8);

MsysControl(MilBoard[0].MilSystem, M_UART_PARITY+M_UART_NB(M_DEV1), M_DISABLE);

MsysControl(MilBoard[0].MilSystem, M_UART_STOP_BITS+M_UART_NB(M_DEV1), 1);

MsysControl(MilBoard[0].MilSystem, M_UART_SPEED+M_UART_NB(M_DEV1), 9600);

MsysControl(MilBoard[0].MilSystem, M_UART_TIMEOUT+M_UART_NB(M_DEV1), M_INFINITE); // ms

In serial-communication of VCC-F22V39ACL, the format of serial frame is:

* Total 14bytes

* Function settings

Trilinear coordination 이란 ?

2011-07-28T01:28:00.000-07:00

삼각형 ABC에 대한 주어진 어떤 점 P의 trilinear coordinates라고 표현한다. trilinear coordinate를 표현하기 위해선 하나의 삼각형과 하나의 점이 주어져야한다.

그럼 trilinear coordinate는 주어진 점 P에서 삼각형의 세 변의 수직거리에 해당하는 α, β, γ의 비로 나타나게 된다. 여기서, α 는 점 P와 변 BC의 수직거리이고, β는 점 P와 변 CA의 수직거리, 그리고 γ는 점 P와 변 AB의 수직거리이다.

이와 같이 α : β : γ는 homogeneous coordinates 나 "trilinears"라고 부른다. 그리고 trilinear coordinate는 1835년에 Plücker에 의해서 소개되었다. 간단히 표현해보면 삼각형 ABC의 세 꼭지점에 해당하는 A, B, C는 1:0:0, 0:1:0 그리고 0:0:1의 trilinear coordinate로 표현된다.

참조 : http://mathworld.wolfram.com/TrilinearCoordinates.html

리눅스에서 파일 정보 캐시 만들고 빠르게 검색하기

2011-07-27T03:06:00.000-07:00

# updatedb

# locate $(your_file_name)

삼각형의 면적 구하기 (Heron's formula)

2011-07-27T02:54:00.000-07:00

Heron's formula에 의하면 삼각형 ABC의 넓이는 다음과 같다. 단, s는 (a+b+c)/2로 정의한다.

위와 같은 식은 다음과 같이 유도할 수 있다.

삼각형 ABC를 수직으로 분할하는 두 개의 삼각형에 대한 피타고라스의 정리는 다음과 같은 관계식을 제공한다.

그리고 두 수식으로부터 다음과 같은 관계식을 유도할 수 있다.

그럼 삼각형 ABC의 넓이는 c*h/2이고 이 넓이의 제곱값은 다음과 같이 전개할 수 있다.

그리고 s=(a+b+c)/2로 정의하면 Heron's의 방정식은 증명된다.

/* By Heron's formula, the area of the triangle
sqrt(s*(s-a)*(s-b)*(s-c)), s: (a+b+c)/2 */
real_t triangle_area(triangle_t *tri)
{
real_t s, a, b, c;

assert(tri);

a = distance_between_points(tri->b, tri->c);
b = distance_between_points(tri->c, tri->a);
c = distance_between_points(tri->a, tri->b);
s = (a + b + c) / 2;

return sqrt(s * (s - a) * (s - b) * (s - c))
}

랩탑 시스템에서 하드디스크 끍는 소리 제거하기

2011-07-27T01:29:00.000-07:00

알려져 있다시피 랩탑용 하드 디스크는 저전력을 위해서 디스크 암이 아무런 작업이 없을때 수시로 초기위치로 되돌아 갔다 트랙으로 갔다를 반복하게 "틱틱틱"하는 소음이 발생한다.

윈도우즈의 경우는 하드디스크 튜닝 프로그램을 항시 상주시켜서 APM(Advanced Power Management)모드를 한단계 아래인 "254"로 설정 시키면 소음이 제거된다고 보고 되고 있다.

리눅스의 경우는 hdparm이라는 명령어를 적절히 사용하면 동일한 효과를 볼 수 있다.

# hdparm -B 254 /dev/hda1

부팅하게 되면 하드 디스크의 APM은 128로 설정되어 있던 것이 위의 명령어 이후에는 254로 설정된다. 이것을 부팅시 254로 영구적으로 설정되도록 하기 위해서 /etc/hdparm.conf파일에 다음과 같은 내용을 추가하자.

/dev/your-hard-drive {
apm = 254
}

그리고 설정을 바로 적용하기 위해서 다음과 같은 명령어를 실행 시키도록 하자.

# update-rc.d hdparm defaults

한가지 더 손바주어야 할 부분은 랩탑의 경우 절전모드나 최대절전모드에서 되돌아 올 경우 하드디스크의 원래의 설정으로 되돌아 가버리기 때문에 /etc/pm/sleep.d/50_hdparm-pm 파일을 만들고 다음과 같은 내용을 추가하자.

#!/bin/sh
if [ -n "$1" ] && ([ "$1" = "resume" ] || [ "$1" = "thaw" ]); then
hdparm -B 254 /dev/your-hard-drive > /dev/null
fi

그리고 실행가능하도록 다음과 같은 명령어를 실행하자.

# chmod +x /etc/pm/sleep.d/50_hdparm-pm

별 효과가 없다면 다음과 같은 내용을 갖는 "99-hdd-spin-fix.sh" 스크립트를 생성하자.

#!/bin/sh
# Use a less aggressive hard disk power management to get rid of
# clicking noise when the drive is parking its heads
hdparm -B 254 /dev/sda

이 스크립트를 /etc/acpi/suspend.d/, /etc/acpi/resume.d/, /etc/acpi/start.d/ 밑에다가 복사하고 실행가능하도록 해주자.

# chmod +x /etc/acpi/?????.d/99-hdd-spin-fix.sh

참조 : http://www.thinkwiki.org/wiki/Problem_with_hard_drive_clicking#Possible_solutions_.28Linux.29

삼각형의 외접원 구하기

2011-07-25T10:03:00.000-07:00

삼각형의 세 꼭지점 A, B, C에 대한 외접원의 반지름 r과 중심점 O는 다음과 같이 구할 수 있다.

위의 수식에서는 삼각형의 꼭지점 A, B, C의 벡터와 외접원의 중심점 O의 벡터로 표현되었다.

/* tri : given triangle
* org : the center of circumscribed-circle
* return value : the radius of circumscribed-circle */
real_t triangle_circumcircle(point_t *org, triangle_t *tri)
{
real_t area, a, b, c;
real_t alpha, beta, gamma;
point_t *A, *B, *C;
point_t *d1, *d2, *d3;

assert(org);
assert(tri);

A = tri->a, B = tri->b, C = tri->c;
area = parallelogram_area_of_3points(A, B, C);
a = distance_between_points(B, C);
b = distance_between_points(C, A);
c = distance_between_points(A, B);

d1 = point_new();
d2 = point_new();

point_subtract(d1, A, B);
point_subtract(d2, A, C);
alpha = a * point_dotproduct(d1, d2) / (2 * sqr(area));

point_subtract(d1, B, A);
point_subtract(d2, B, C);
beta = b * point_dotproduct(d1, d2) / (2 * sqr(area));

point_subtract(d1, C, A);
point_subtract(d2, C, B);
gamma = c * point_dotproduct(d1, d2) / (2 * sqr(area));

point_destroy(d2);
point_destroy(d1);

org->x = alpha * A->x + beta * B->x + gamma * C->x;
org->y = alpha * A->y + beta * B->y + gamma * C->y;
org->z = alpha * A->z + beta * B->z + gamma * C->z;

return a*b*c/(2*area);
}

리눅스 CAS(Computer Algebra System)설치 (Maple같은 수식 전개방식의 수학 프로그램)

2011-07-24T03:26:00.000-07:00

삼각형의 외접원의 중심점과 반지름을 구하기 위해서 wiki를 검색해 보았다. circumscribed circle, 혹은 circumcenter등으로 너무나도 친절하게도 잘되어 있다. 그런데 하는 김에 3차원으로 되어 있는 것으로 프로그램을 작성해 보고 싶어 직접 수식을 구성해서 해보기로 했다.

이와 같이 간혹 수식을 종이에 전개해야할 일들이 종종 있다. 머리와 종이도 아낄겸 예전에 windows에서 사용하던 Maple을 리눅스에서 사용하기 위해서 찾아 보았다.

두개의 CAS(Computer Algebra System)이 검색되었다. MuPAD와 MAXIMA라는 프로그램이다. 우선은 MuPAD를 설치하기위해서 데비안 리눅스 패키지를 찾아 보았으나 없었다. 보아하니 상용으로 전환하는 수순이가 보다 Mathworks에 넘어가 MATLAB에서 툴박스로 사용할 수 있게 되어 있다. 예전부터 MATLAB에서 수식 전개와 같은 범용 계산이 가능했으면 했는데 이미 사들였던 것이다.

그래서 MAXIMA로 눈을 돌렸다. 텍스트 버전도 있고 GUI버전에다가 emacs에서 사용할 수도 있게 제공하는것같다. 무엇보다도 데비안 리눅스 배포본을 착실하게 제공하는 것같다. 우선은 설치해봤다.

Maxima Algebra System
Maxima backend for Cantor
A computer algebra system -- documentation
A computer algebra system -- emacs interface
A computer algebra system -- source code
A computer algebra system -- test suite
A computer algebra system -- base system
A computer algebra system -- extra code

여러 개가 다수 존재한다. 우선 Maxima Algebra System만 선택해서 설치하도록 하고 도움 문서와 base system 정도만 선택해서 설치해 보기로 했다. 그랬더니 documentation, source code, test suite, base system, extra code가 같이 설치 되어 버렸다. 우선은 써봐야 겠다.

프로그램의 교육에 설치가 되어있다. 우선은 삼각형의 외접원의 원점을 구하는 방정식으로 시험해 보기로 했다.

결과는 보았으나 가져다 쓰기가 영 좋지가 못하다. 이것도 감지덕지이긴 하지만.

Double linked list in c language

2011-07-18T10:03:00.000-07:00

* data type double linkage

struct dlink {
void *object;
struct dlink *prev;
struct dlink *next;
};

typedef struct dlink dlink_t;

- object : reference pointer indicating the address for any-type user-data using casting to void type pointer.
- prev : previous pointer for double-linked component
- next : next pointer for double-linked component

* data type double linked list

typedef struct {
int count;
dlink_t *head;
dlink_t *tail;
} dlist_t;

- count : the number of components contained by the list.
- head, tail : the start and end linkage for data- inport and outport.

----------------------------------- dlink.h -------------------------------------------

/* Double linked list

This program is free software; you can redistribute it and/or modify it

under the terms of the GNU Lesser General Public License as published by

the Free Software Foundation; either version 2 of the License, or

(at your option) any later version.

This program is distributed in the hope that it will be useful, but

WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY

or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public

License for more details.

You should have received a copy of the GNU Lesser General Public License

along with this program; if not, please visit www.gnu.org.

#ifndef __DLINK_H__
#define __DLINK_H__

struct dlink {
void *object;
void *private;
struct dlink *prev;
struct dlink *next;
};

typedef struct dlink dlink_t;

dlink_t *dlink_new(void);
void dlink_prepend(dlink_t *obj, dlink_t *at);
void dlink_append(dlink_t *obj, dlink_t *at);
void dlink_cutoff(dlink_t *obj);
void dlink_substitute(dlink_t *dummy, dlink_t *at);
void dlink_exchange(dlink_t *a, dlink_t *b);
void dlink_destroy(dlink_t *dlink);

typedef struct {
int reference;
int count;
void *private;
dlink_t *head;
dlink_t *tail;
} dlist_t;

dlist_t *dlist_new(void);
void dlist_push(dlink_t *link, dlist_t *list);
dlink_t *dlist_pop(dlist_t *list);
void dlist_insert(dlink_t *link, dlist_t *list);
dlink_t *dlist_extract(dlist_t *list);
dlink_t *dlist_glimpse(int index, dlist_t *list);
dlink_t *dlist_pick(int index, dlist_t *list);
void dlist_put(dlink_t *link, int i, dlist_t *list);
void dlist_exchange(int i, int j, dlist_t *list);
void dlist_destroy(dlist_t *list);

#endif /* __DLINK_H__ */
----------------------------------- dlink.c -------------------------------------------

/* Double linked list

This program is free software; you can redistribute it and/or modify it

under the terms of the GNU Lesser General Public License as published by

the Free Software Foundation; either version 2 of the License, or

(at your option) any later version.

This program is distributed in the hope that it will be useful, but

WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY

or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public

License for more details.

You should have received a copy of the GNU Lesser General Public License

along with this program; if not, please visit www.gnu.org.

#include
#include
#include
#include "dlink.h"

dlink_t *dlink_new(void)
{
dlink_t *dlink;

dlink = malloc(sizeof(dlink_t));
assert(dlink);
memset(dlink, 0, sizeof(dlink_t));

return dlink;
}

void dlink_prepend(dlink_t *obj, dlink_t *at)
{
assert(obj);
assert(at);

obj->prev = at->prev;
obj->next = at;

at->prev->next = obj;
at->prev = obj;
}

void dlink_append(dlink_t *obj, dlink_t *at)
{
assert(obj);
assert(at);

obj->prev = at;
obj->next = at->next;

at->next->prev = obj;
at->next = obj;
}

void dlink_cutoff(dlink_t *obj)
{
assert(obj);

obj->prev->next = obj->next;
obj->next->prev = obj->prev;
}

void dlink_substitute(dlink_t *dummy, dlink_t *at)
{
assert(dummy);
assert(at);

at->prev->next = dummy;
at->next->prev = dummy;
dummy->prev = at->prev;
dummy->next = at->next;
}

void dlink_exchange(dlink_t *a, dlink_t *b)
{
dlink_t *dummy;

assert(a);
assert(b);

dummy = dlink_new();

dlink_substitute(dummy, a);
dlink_substitute(a, b);
dlink_substitute(b, dummy);

dlink_destroy(dummy);
}

void dlink_destroy(dlink_t *dlink)
{
assert(dlink);
free(dlink);
}

dlist_t *dlist_new(void)
{
dlist_t *dlist;

dlist = malloc(sizeof(dlist_t)+2*sizeof(dlink_t));
assert(dlist);

memset(dlist, 0, sizeof(dlist_t)+2*sizeof(dlink_t));

dlist->head = (dlink_t *)(dlist+1);
dlist->tail = dlist->head+1;

dlist->head->next = NULL;
dlist->head->prev = dlist->tail;

dlist->tail->next = dlist->head;
dlist->tail->prev = NULL;

return dlist;
}

void dlist_push(dlink_t *link, dlist_t *list)
{
assert(list);
assert(link);

dlink_append(link, list->tail);
//dlink_prepend(link, list->head);
list->count++;
}

dlink_t *dlist_pop(dlist_t *list)
{
dlink_t *link;
assert(list);

//link = list->head->prev;
link = list->tail->next;
dlink_cutoff(link);
list->count--;

return link;
}

void dlist_insert(dlink_t *link, dlist_t *list)
{
assert(list);
assert(link);

dlink_prepend(link, list->head);
list->count++;
}

dlink_t *dlist_extract(dlist_t *list)
{
dlink_t *link;

assert(list);

if (list->count == 0) return NULL;

//link = list->tail->next;
link = list->head->prev;
dlink_cutoff(link);
list->count--;

return link;
}

dlink_t *dlist_glimpse(int index, dlist_t *list)
{
dlink_t *link;

assert(list);
assert(index >= 0 && index < list->count);

link = list->tail->next;
while (index--) link = link->next;

return link;
}

dlink_t *dlist_pick(int index, dlist_t *list)
{
dlink_t *link;

assert(list);
assert(index >= 0 && index < list->count);

link = list->tail->next;
while (index--) link = link->next;
dlink_cutoff(link);
list->count--;

return link;
}

void dlist_put(dlink_t *link, int i, dlist_t *list)
{
dlink_t *ilink;

assert(link);
assert(i >= 0 && i < list->count);

ilink = list->tail->next;
while (i--) ilink = ilink->next;
dlink_prepend(link, ilink);
}

void dlist_exchange(int i, int j, dlist_t *list)
{
dlink_t *ilink, *jlink;

assert(list);
assert(i >= 0 && i < list->count);
assert(j >= 0 && j < list->count);

ilink = list->tail->next;
while (i--) ilink = ilink->next;

jlink = list->tail->next;
while (j--) jlink = jlink->next;

dlink_exchange(ilink, jlink);
}

void dlist_destroy(dlist_t *list)
{
dlink_t *link;

assert(list);

if (list->reference == 0) {
while (list->count > 0) {
link = dlist_extract(list);
dlink_destroy(link);
}
free(list);
} else list->reference--;
}

Use GdkPixbuf in GTK+

2011-07-16T09:32:00.000-07:00

At first it needs to create GtkImage on GtkViewport. Then, we can use GdkPixbuf-embedded GtkImage.

static void gdk_pixbuf_clear(GdkPixbuf *pixbuf)
{
int x, y;
int width, height, rowstride, n_channels;
guchar *pixels, *p;

n_channels = gdk_pixbuf_get_n_channels(pixbuf);
g_assert(gdk_pixbuf_get_colorspace(pixbuf) == GDK_COLORSPACE_RGB);
g_assert(gdk_pixbuf_get_bits_per_sample(pixbuf) == 8);
g_assert(gdk_pixbuf_get_has_alpha(pixbuf));
g_assert(n_channels == 4);

width = gdk_pixbuf_get_width(pixbuf);
height = gdk_pixbuf_get_height(pixbuf);
rowstride = gdk_pixbuf_get_rowstride(pixbuf);
pixels = gdk_pixbuf_get_pixels(pixbuf);

for (y = 0; y < height; y++) {
for (x = 0; x < width; x++) {
p = pixels + y*rowstride + x*n_channels;
*(p+0) = 255; // red channel
*(p+1) = 0; // green channel
*(p+2) = 0; // blue channel
*(p+3) = 255; // alpha channel
}
}
}

....

image = gtk_image_new();

....

pixbuf = gdk_pixbuf_new(GDK_COLORSPACE_RGB, TRUE, 8, 640, 480);
g_assert(pixbuf);

gdk_pixbuf_clear(pixbuf);

gtk_image_set_from_pixbuf(image, pixbuf);

...

FabFi - 무료, 공개된 오픈 소스 장거리 wireless ethernet 전송 시스템(Reflector + Router)

2011-07-14T19:57:00.000-07:00

[출처] http://fabfi.fablab.af

FabLab의 공용으로 성장중인 시스템이다. 제목에서 보듯이 장거리 wireless ethernet 신호를 전송해 줄 목적으로 된 소스가 공개된 시스템이다. 서울에서는 필요없을 수도 있겠지만 wireless 네트워크가 미약한 지방의 경우 활용 용도가 분명히 있을 것이다. 그래도 reflector는 충분히 wireless 전파가 미약한 곳에서 RF 전파를 포집하는데 많은 도움이 될 수 있는 안테나 증폭기로 사용될 수 있을 듯 싶다. 찾아 보면 reflector의 도면 파일이 있다. 현재 케냐와 아프가니스탄 같은 기반시설이잘 갖추어 지지 않은 곳에서 구현되고 있다.

현제 wiki에서 제공하는 내용은 맨 처음으로 구성된 아프가니 스탄의 Jalalabad의 FabFi를 기반으로 설명하고 있다. 이 시스템은 uplink와 Open courseware 미러링을 제공하고 링크할 수 있는 거리는 3,880m이다.

FabFi의 안테나는 reflector와 router로 구성된다. reflector는 직접 개개인이 만들어야 하고, router는 수정된 기존 상품이다. reflector 디자인은 4인치에서 18인치까지로 범위가 다양한다. (plywood, signboard, 아크릴로 만든다.) radio frequency(RF)에 대한 반사 재료로 chicken wire(?), 스텐레스 짠 그물망, 그리고 강철로 된 윈도우 스크린를 사용한다. 특히 똑똑한 아프가니스탄 사람들은 오래된 USAID 쿠킹 오일 캔을 reflector로 사용하기 시작했다.

router는 쉽게 판매되는 가정용 wireless AP에 수정된 router-radio-antenna를 사용한다. 버전 3.0에서는 FabFi는 공식적으로 Linksys WRT54GL과 ASUS WL-520GU를 지원한다. 비록 다른 많은 router들도 그에 걸맞는 image를 생성하면 쉽게 사용할 수 도 있다. FabFi는 OpenWRTOpenWRT 8.09.2 (Kamikaze) firmware가 구동되는 상황에서 구동된다.

FabFi 기술 요약

개요

FabFi은 확장가능한 장거리 P2P 시스템이고 mesh hybrid-wireless broadband transmission 기반시설이다. 개념은 간단하다. 간단하고 지능적이고 서로 연결가능한 network가 불안정한 상황에서도 믿을 수 있는 network 환경을 만들어 낼 수 있다는 개념이다. FabFi는 low-cost 장치를 아주 먼거리에 대해서 직접적으로 통신할 수 있도록 하는 간단한 물리 개념과 아주 다양한 상황들에 적용할 수 있는 탄력적인 구성을 활용한다.

극한의 상황에 대해서 RF(Radio Frequency) 에너지를 방향 전환시킬 수 있는 wire mesh surface를 갖는 제작된 RF reflector위에 상용 wireless routers를 장착하여 사용한다. Reflector의 증폭율은 reflector의 크기와 사용된 재료에 의존한다. 그러나 현재의 디자인에서는 15dBi정도로 높게 측정되어진다.

FabFi 시스템에서 단일 wireless link는 wireless router가 장착된 두 대의 reflector로 구성된다. 이와 유사하게, 두 대의 router도 wired 연결로 link될 수 있다. 단일 router는 동시에 wireless 연결과 wired연결 둘 다로 link될 수도 있다. 이 시스템은 수많은 방법으로 구성된 개개의 link에 대한 구성도 가능하다. (매우 먼 거리나 작은 지역에 많은 사용자를 수용할수 있는 링크 같은 것들 말이다. 이 링크의 중요한 요소는 "meshing"이라는 것이다. mesh network는 구조화 되지 않는 ad-hoc같은 방식으로 하나 이상의 다른 근방의 디바이스에 접속할 수 있도록 허용하는 시스템이다. Mesh network는 소프트웨어가 network의 topology(구조)를 감지해서 실시간으로 자동적으로 데이터의 흐름을 결정하기 때문에 강인하고 단순하다. 전통적인 mesh network들은 단일 라디오와 wireless 만으로된 연결과 모든 방향을 커버할 수 있는 안테나에 의존하기 때문에 구조적으로 제한되어져 있다. 가능한 상황에 따라서 직접적인 wireless link나 유선 전송을 이용해서 Fabfi 시스템은 매우 큰 규모의 한곳에 고정되서 사용하는 mesh network를 구성할 수 있도록 최적화 되어 있다. 높은 데역폭의 좁은 영역의 uplink에 의존하지 않고 지역내에서 강인한 디지털 통신이 가능한 규모에서는 가격이 비싸고 많은 지역에서 사용 불가능하다. 좀 더 자세한 설명을 해주는 에니매이션을 확인해보자. (2MB의 용량으로 다소 시간이 걸릴수 있다.)

Reflector의 작동 방식

FabFi의 reflector는 포물선의 준선에 수직한 전파가 포물선의 표면에 접촉할 때 이 전파들을 포물선의 초점에 모이게 하는 포물선 모양(Y=cX^2)의 구조를 사용한다.(이것에 대한 좀더 자세한 것은 Mathworld 를 보라...) RF 전파를 반사 시킬 수 있는 window screen, chicken wire같은 것을 3차원 의 포물선 모양을 이루는 frame에 붙이고 정확하게 계산된 추점에 wireless router를 붙인다. 그럼 전송 지점의 router에서 오는 RF 에너지가 수신부의 router에서 RF 에너지가 효율적으로 모아진다.

Reflector 만들기

FabFi 시스템의 본질적인 부분은 지역특색에 따라 사용가능한 어떤 재료로도 구현할 수 있는 유연성이 있다는 것이다. 필요한 모든 것들은 2D 디자인 파일을 프린트하고 여러분이 찾을 수 있는 재료를 조각으로 만들 수 있는 능력이면 된다. 여러분이 Fab Lab을 가지고 있다면. 나무나 금속, 아크릴로 직접 reflector를 만들기 위해서 CNC 나무 제단기나 레이져 커터를 사용할 수 있다. 그러나 진흙으로 빚거나 돌을 깍거나 얼음 블럭을 끌로 파거나 해서 만들지 못할 이유는 없다. 2009년 1월에 초기 배치시 Jalalabad에서는 세 가지의 다른 reflector 디자인이 구현되었다: 큰 4인치의 나무로 된것과 2인치의 나무 버전 그리고 18인치의 아크릴 버전이 구현되었다. 반사면의 재질은 chicken wire와 스테인레스 강철로 짜여진 그물망과 window screen이 될 수 있다.

다음으로 실제 현장에서 필요한 것들은 봉합된 방수가 가능한 뺏다 낄 수 있는 조립품의 수정된 reflector 디자인을 개발하는 것이다. 이 새로운 디자인 2010년 여름에 개시 되었다.

그러나 멀리 볼 것도 없이 network 사용자들은 이미 지역적으로 구할 수 있는 고철 더미로 만든 reflector를 디자인하고 만들고 있다. 아직 중요한 부분에 정교한 장업이 필요하지만, 이 reflector들은 기술적인 면에서 지역 인력이 자체적으로 개발되는 기술이 전파되고 있다는 분명한 신호이다. 이것들은 모두 $3달러 이하의 가격이다.!

Router와 FabFi 펌웨어는 OpenWRT라 불리는 open source 3rd party 펌웨어를 사용한다. OpenWRT는 리눅스 기반의 유연성을 장점으로 가지고 있어 FabFi 장치들은 다양한 범위의 모니터링과 자가 진단 도구들을 운용할 수 있다. 현재의 시스템은 실시간 network 모니터링과 local web을 저장해 집중된 접근을 제어하는 기능과 사용자 관리, 요금 부과를 위한 사용량 추적을 지원한다. 이 모든 것들은 $50에서 $100 달러의 장치에서 구현된다. 자동 구성 기능은 FabFi 1.0이 발표된 이후로 지속적으로 개선되어 왔다. 지금은 다중 fabfi상에서 다중 router를 지원하고 있다. 그리고 802.11n의 속도로 network를 구성할 수 있는 능력이 있다. 개발 부서에서 사용하는 전원인 안정적인 전원은 계속적인 도전거리다. 편리하게도 지금 지원하는 모든 장치들은 12VDC에서 동작하고 자동차나 작은 엔진 베터리에서 직접 전원으로 사용 될 수 있다. 자동차 배터리와 두 대의 비싸지 않은 충전기는 Jalalbad network에서 두 개의 큰 distribution 허브에서 작동하는 믿을 만한 UPS 장치로 작동한다. 도시의 전원없이 거의 이틀동안 router 를 동작시킬수 있는 전원이다. 케냐에서 주 전원과 태양열 전원을 모두 지원하는 network에 어느 node에나 부착할 수 있고 방수 기능이 있는 UPS를 지원하는 "node in a box"를 고안했다. 다음 개발 버전은 router에 연결된 일반 UPS나 제작된 UPS에 장착하는 100-240VAC switching power를 이용해서 설치할 수 있도록 하는 bare-bone 타입의 12V-12V UPS가 될 것이다. 바람이나 다른 그 지역에서 생산되는 전력으로 동작하는 충전회로들은 동시에 진행되는 FabLab의 프로젝트이다.

Fab Future는 미적인 부분도 고려되었을 뿐만 아니라, Fabfi는 일상적으로 모래 폭퐁이 부는 여름의 Jalalabad에서 130F의기온인 혹독한 아프가니스탄 기후에서도 놀랄정도로 안정적인 것이 증명되어 왔다. 2010년도 초여름에 fabfi 시스템을 client 장치에 직접 wireless 접속하도록 제공하고 통신망 규모의 wifi ISP업체를 운영해 왔다. 2년 이상의 운영하는 동안 한편으로는 여전히 하드웨어적인 실폐가 발견되기도 한다. 앞으로의 가장 큰 과제는 uplink 대역폭 문제이다. 많은 나라들은 범용적인 접속 문제의 해결방안으로써 "mobile broadband"를 밀고 있지만, 대부분의 지역에서의 실제적인 현실은 mobile 장치들만으로는 어떤 보조없이 감당할만한 충분한 가격이나 충분한 효율을 제공하지 못한다. 케냐에서 fabfi는 mobile 연결이 사용자당 데이타 허용량이 감소해서 단지 접속의 도구로 전락한 상황에서 통신하는데 가치있는 서비스를 제공한다.

이것에 관한것이다 기타 다른 프로젝트의 진행 과정은 여기 FabFi Blog에서 살펴 볼 수 있다.

reflectror 만들기

http://fabfi.fablab.af/distribution.html에 가서 알맞은 크기의 디자인을 다운로드 받아서 나무를 자르고 조립한다. 원하면 페인트를 칠하거나 실링을 할 수 도 있다. 그리고 조립된 프레임에 철로된 그물망이나 스크린을 고정시킨다. 고정면이 구부러지거나 주름이 생기지 않도록 단단하게 해야한다.

라우터 구성하기

안테나 설치하기

reflector에 router를 고정시킨다. 비가 올 경우를 대비해서 router주변을 비닐 봉다리로 둘러 싸거나 방수가 되는 플라스틱 박스를 사용한다. router에 전원을 인가하고 LAN포트에 인터넷이 연결될 수 있는 ethernet 케이블을 연결시킨다. 다른 FabFi의 방향으로 안테나를 조절한다. 컴퓨터의 LAN포트에 연결해서 연결상태를 점검한다. 안테나를 끈을 이용해서 단단히 고정한다.