数値計算

Size: px

Start display at page:

Download "数値計算"

ふじよしなかきむら
4 years ago
Views:

1 プログラム作成から実行まで数値計算垣谷公徳 17 号館 3 階電子メール : kimi@ee.ous.ac.jp Source program hello.c printf("hello\n"); コンパイラ Library libc.a 0011_printf _getc _scanf コンパイル Object module hello.o _main_ _printf_0 リンカリンク Executable module a.out プログラム作成実行手順原始プログラム (Source program) の作成原始プログラムをコンパイルし目的プログラム (Object module) に変換目的プログラムを連係編集 (link) し実行可能プログラム (Executable module/load(able) moduleに変換リンク ( 連係編集 ) リンケージエディタ ( リンカ ) によって目的プログラムのサブルーチン関数のアドレス定義を行いライブラリなどと結合させて実行可能プログラムを作成する実行可能プログラムをメモリ上にロードし実行するプログラム作成用語集原始プログラム (source program/source code/source file) プログラム言語によって記述されたもの目的プログラム (object module/object code/object file) コンパイラによって変換されたものアドレスなどが未定義のため実行不可実行可能プログラム (executable file/load module/command) リンカ ( リンケージエディタ ) によりプログラムの実行に必要なサブルーチンや関数などを含めた計算機で実行可能なもの 3 4 プログラミング言語数値計算に用いられる主なプログラミング言語 FORTRAN 66/77/90/95 c (c++, java) python etc. FORTRAN 言語の特徴最古の高水準言語汎用構造化言語 (FORTRAN77 以降 ) FORTRAN90 以降はオブジェクト指向数値計算に適した言語仕様過去の豊富な資産 ( ライブラリ ) 数値計算の分野では圧倒的なシェア 5 6

2 c 言語の特徴低水準な処理も可能な高水準言語汎用構造化言語オブジェクト指向言語 (c++, java) の基礎ソフトウェア技術者の教養的知識 c 言語数値計算のための c 言語 (1) 2013 年 4 月 16 日最新のライブラリ ( 過去の資産は使えない ) 数値計算には不向き (FORTRAN に比べて ) 7 8 Linux 等でコンパイル実行する場合 $ gcc this_file.c -lm 用意されているキーワード ( 予約語 ) 特別な意味を持つ記号 9 10 ユーザ定義名前付きオブジェクトユーザ定義リテラル 11 システムに依存するリテラル標準入出力ライブラリ標準数学ライブラリ 12

3 プリプロセッサにより1 行 ( レコード ) 毎に処理される # キーワードのフィールドから始める 13 コンパイラで処理される ( ホワイトスペース ) ; { /* */ で区切られたバイト列行の概念なし 14 文とブロックブロック double x, y; printf( Hello, world.\n ); x = 0.1; y = sin(x); 文 printf( SIN(%lf) = %lf\n, x, y); return 0; # プリプロセッサ命令 # 定数マクロ定義 int main( コマンドライン引き数 ) { 変数の型変数名 ; /* 変数宣言 */ 実行すべき処理 ; return 0; プログラミング言語の一般論データ型 ( 定数と変数配列 ) 代入基本演算 ( 四則演算 ) 入出力条件分岐繰返処理関数外部手続きキーワード ( 予約語 ) データタイプ char int float double enum struct union typedef データタイプ修飾子 signed unsigned short long const void データ記憶クラス auto static extern volatile register 制御構造 if else do while for switch case break default return continue goto 演算子, 記号 sizeof + - * / % &, ; : [ ] ( ) { # 17 18

4 変数データ ( ビット列 ) を入れることができる箱名前がついている大きさが決まっている実行時に用意される ( 実行中ではない ) 'a' c 1026 n pi 文字型 [ 通常 8bit] 変数 char 整数型 [ 通常 16bit or 32bit] int 拡大整数型 [ 通常 32bit 以上 ] long int, long 実数型 ( 単精度浮動小数点数 )[ 通常 double と同一 ] float 倍精度実数型 ( 倍精度浮動小数点数 )[ 通常 64bit] double 変数 Local Rule.1: 変数は全て32bit 符号付き整数 (int or long int) と 64bit 実数 (double) を使う 32bit 符号付き整数 ~ bit 実数有効桁数約 16 桁の実数 21 変数宣言データ型変数名 ; データ型変数名, 変数名,...; データ型変数名 = 初期値 ; 整数型倍精度実数型 int a; int a, b; int n=10; double x; double x, y, z; double pi = ; 22 整数 ( 拡大整数 ) 定数 0, -1, 1234,( L) 倍精度実数型 ( 倍精度浮動小数点数 ) 0.0, 1., 123.4, 1.234e e-3 = = 文字列定数 ( 文字列リテラル ) "Hello world\n" 23 記号定数 ( マクロ置換 ) 表に名前裏に値 ( 定数 ) の書かれたカードプリプロセッサの機能を利用する厳密な意味ではユーザ定義名前付きオブジェクトではない #define 記号名置換すべきテキスト #define N 1000 #define PI PI double pi2; pi2 = 2.*PI; 24

5 四則演算と代入変数 = 変数 ( 定数 ) 演算子変数 ( 定数 ) 和 y = a + b y = a + b 差 y = a b y = a - b 積 y = a b y = a * b 商 y = a b y = a / b 剰余 y = a mod b y = a % b 四則演算と代入整数 = 整数 + 整数 ; 整数 = 整数 - 整数 ; 整数 = 整数 * 整数 ; 整数 = 整数 / 整数 ; 整数 = 整数 % 整数 ; 実数 = 実数 + 実数 ; 実数 = 実数 - 実数 ; 実数 = 実数 * 実数 ; 実数 = 実数 / 実数 ; Local Rule.2: 代入文の右辺と左辺の型を明示的に一致させる自動的な型変換を用いず, 常に明示的にキャストを行う int n; double a, x; x = a / n; int n; double a, x; x = a / (double) n; 配列配列名と数字 ( 添字 ) で指定される変数の集まり double a[100]; double b[10][5]; a[0] = 1.0; a[1] = 2.0; a[99] = 5.0; b[0][0] = 1.2; b[9][4] = 0.0; 'H' c[0] 'e' c[1] 'l' c[2] 'l' c[3] 'o' c[4] m[0][0] m[0][1] m[0][2] m[1][0] m[1][1] m[1][2] m[2][0] m[2][1] m[2][2] m[3][0] m[3][1] m[3][2] \0 c[5] 一次元配列とベクトル double b[3]; b[0] = 1.0; b[1] = 2.0; b[2] = 1.0; double b[3]={1.0, 2.0, 3.0; double b[3]= {1.0, 2.0, 3.0; B = 一次元配列と数列 #define N 100 二次元配列と行列 double a[2][2]; a[0][0] = 1.0; int n; double a, r, c[n]; c[n] = a*pow(r, n-1); c[0] = a; c[n] = r*c[n-1]; C n = ar n 1 { C0 = 0 C n = rc n 1 a[0][1] = -1.0; a[1][0] = -1.0; a[1][1] = 2.0; double a[2][2]={{1.0, -1.0, {-1.0, 2.0; double a[2][2]= {{1.0, -1.0, A = ( ) {-1.0, 2.0 ; 29 30

6 XX00 x[0] XX06 x[1] XX00 x 配列に関する注意 int x[2][3] = { {0, 1, 2, {4, 8, 16 ; XX00 x[0][0] XX02 x[0][1] XX04 x[0][2] XX06 XX08 XX0A x[1][1] x[1][2] x[1][0] x[1][1] x[1][2] XX0C Local rule: 配列の要素のアドレスを直接意識するようなプログラムにはしない必ずすべての ( 二次元なら2つ三次元なら3つ ) の添え字を指定する 31 入出力入力関数 getchar, getc, fgetc ; 一文字入力 gets, fgets, fread ; 文字列 ( データ列 ) 入力 scanf, fscanf, sscanf ; 書式付入力出力関数 putchar, putc, fputc ; 一文字出力 puts, fputs, fwrite ; 文字列 ( データ列 ) 出力 printf, fprintf, sprintf; 書式付出力 Local Rule.3: 入力は行わず, データはプログラム中に直接書く. 出力は printf のみ. 32 printf 関数 int printf (const char *, 変数, 変数,...); 書式制御文字列リテラル %\" 以外はそのまま出力 "%+ 数字 + 文字 " %d 整数を出力 %f 固定小数点実数を出力 %e 浮動小数点実数を出力 "%%" % 自体を出力その他 %i %o %u, %c, %s %x, %X %E, %g, %G 書式制御文字列 int n = 578; printf("%d\n", n); printf("%8d\n", n); printf("%08d\n", n); printf("%2d\n", n); 書式制御文字列 double x = ; 書式制御文字列 double x = ; printf("%f\n", x); printf("%9.3\n", x); printf("%e\n", x); e printf("%15.7e\n", x); e printf("%10.1e\n", x); e printf("%6.3d\n", x); printf("%6.1d\n", x); printf("%27.20e\n", x); e

7 printf 関数 int printf (const char *, 変数, 変数,...); 書式制御文字列リテラル Local Rule.4 原則 %d 整数 %f 実数 %15.7e 実数を使用する 37

数値計算

数値計算垣谷公徳 17 号館 3 階電子メール : kimi@ee.ous.ac.jp プログラミング言語の一般論データ型 ( 定数と変数配列 ) 代入基本演算 ( 四則演算 ) 入出力分岐繰返処理関数外部手続き 1 2 入力関数入出力 getchar, getc, fgetc ; 一文字入力 gets, fgets, fread ; 文字列 ( データ列 ) 入力 scanf,