行列(ベクトル)

[C言語] ベクトル同士を足し算する方法

2025-01-18更新日: 2025-01-18

C言語でベクトル同士を足し算するには、通常、配列を用いて各要素を個別に加算します。

例えば、2つのベクトルvec1とvec2がある場合、それぞれの対応する要素を加算し、新しいベクトルresultに格納します。

この操作はループを用いて実装され、各インデックスに対してresult[i] = vec1[i] + vec2[i]のように計算します。

この方法は、ベクトルの次元が固定されている場合に特に有効です。

目次から探す

ベクトルの足し算の実装
ベクトルの足し算の応用例
ベクトル演算の最適化
まとめ

ベクトルの足し算の実装

配列を用いたベクトルの足し算

配列を用いたベクトルの足し算は、最も基本的な方法です。

ここでは、2つの同じ長さのベクトルを要素ごとに足し合わせ、新しいベクトルを作成します。

#include <stdio.h>
#define VECTOR_SIZE 3
void addVectors(int vector1[], int vector2[], int result[]) {
    for (int i = 0; i < VECTOR_SIZE; i++) {
        result[i] = vector1[i] + vector2[i];
    }
}
int main() {
    int vector1[VECTOR_SIZE] = {1, 2, 3};
    int vector2[VECTOR_SIZE] = {4, 5, 6};
    int result[VECTOR_SIZE];
    // ベクトルの足し算を実行
    addVectors(vector1, vector2, result);
    // 結果を表示
    printf("Resultant Vector: ");
    for (int i = 0; i < VECTOR_SIZE; i++) {
        printf("%d ", result[i]);
    }
    printf("\n");
    return 0;
}

Resultant Vector: 5 7 9

このプログラムでは、vector1とvector2の各要素を足し合わせてresultに格納しています。

VECTOR_SIZEを変更することで、異なる次元のベクトルにも対応できます。

構造体を用いたベクトルの足し算

構造体を用いることで、ベクトルの要素をより直感的に扱うことができます。

以下の例では、構造体を使って2次元ベクトルの足し算を行います。

#include <stdio.h>
typedef struct {
    int x;
    int y;
} Vector2D;
Vector2D addVectors(Vector2D v1, Vector2D v2) {
    Vector2D result;
    result.x = v1.x + v2.x;
    result.y = v1.y + v2.y;
    return result;
}
int main() {
    Vector2D vector1 = {1, 2};
    Vector2D vector2 = {3, 4};
    // ベクトルの足し算を実行
    Vector2D result = addVectors(vector1, vector2);
    // 結果を表示
    printf("Resultant Vector: (%d, %d)\n", result.x, result.y);
    return 0;
}

Resultant Vector: (4, 6)

このプログラムでは、Vector2Dという構造体を定義し、2次元ベクトルの足し算を行っています。

構造体を使うことで、ベクトルの各成分に名前を付けて管理しやすくなります。

動的メモリを用いたベクトルの足し算

動的メモリを用いることで、実行時にベクトルのサイズを決定することができます。

以下の例では、動的メモリを使ってベクトルの足し算を行います。

#include <stdio.h>
#include <stdlib.h>
void addVectors(int *vector1, int *vector2, int *result, int size) {
    for (int i = 0; i < size; i++) {
        result[i] = vector1[i] + vector2[i];
    }
}
int main() {
    int size;
    printf("Enter the size of the vectors: ");
    scanf("%d", &size);
    // 動的メモリの割り当て
    int *vector1 = (int *)malloc(size * sizeof(int));
    int *vector2 = (int *)malloc(size * sizeof(int));
    int *result = (int *)malloc(size * sizeof(int));
    // ベクトルの入力
    printf("Enter elements of vector1: ");
    for (int i = 0; i < size; i++) {
        scanf("%d", &vector1[i]);
    }
    printf("Enter elements of vector2: ");
    for (int i = 0; i < size; i++) {
        scanf("%d", &vector2[i]);
    }
    // ベクトルの足し算を実行
    addVectors(vector1, vector2, result, size);
    // 結果を表示
    printf("Resultant Vector: ");
    for (int i = 0; i < size; i++) {
        printf("%d ", result[i]);
    }
    printf("\n");
    // 動的メモリの解放
    free(vector1);
    free(vector2);
    free(result);
    return 0;
}

Enter the size of the vectors: 3
Enter elements of vector1: 1 2 3
Enter elements of vector2: 4 5 6
Resultant Vector: 5 7 9

このプログラムでは、mallocを使って動的にメモリを割り当て、ユーザーからベクトルのサイズと要素を入力してもらいます。

動的メモリを使うことで、プログラムの柔軟性が向上します。

ベクトルの足し算の応用例

3次元ベクトルの足し算

3次元ベクトルの足し算は、2次元ベクトルの足し算と同様に、各成分を足し合わせることで実現できます。

以下の例では、3次元ベクトルの足し算を行います。

#include <stdio.h>
typedef struct {
    int x;
    int y;
    int z;
} Vector3D;
Vector3D addVectors(Vector3D v1, Vector3D v2) {
    Vector3D result;
    result.x = v1.x + v2.x;
    result.y = v1.y + v2.y;
    result.z = v1.z + v2.z;
    return result;
}
int main() {
    Vector3D vector1 = {1, 2, 3};
    Vector3D vector2 = {4, 5, 6};
    // ベクトルの足し算を実行
    Vector3D result = addVectors(vector1, vector2);
    // 結果を表示
    printf("Resultant Vector: (%d, %d, %d)\n", result.x, result.y, result.z);
    return 0;
}

Resultant Vector: (5, 7, 9)

このプログラムでは、Vector3Dという構造体を用いて3次元ベクトルの足し算を行っています。

各成分を個別に足し合わせることで、3次元空間でのベクトル演算が可能です。

ベクトルのスカラー倍との組み合わせ

ベクトルのスカラー倍は、ベクトルの各成分にスカラー値を掛ける操作です。

これをベクトルの足し算と組み合わせることで、より複雑なベクトル演算が可能になります。

#include <stdio.h>
typedef struct {
    int x;
    int y;
} Vector2D;
Vector2D scalarMultiply(Vector2D v, int scalar) {
    Vector2D result;
    result.x = v.x * scalar;
    result.y = v.y * scalar;
    return result;
}
Vector2D addVectors(Vector2D v1, Vector2D v2) {
    Vector2D result;
    result.x = v1.x + v2.x;
    result.y = v1.y + v2.y;
    return result;
}
int main() {
    Vector2D vector1 = {1, 2};
    int scalar = 3;
    // スカラー倍を実行
    Vector2D scaledVector = scalarMultiply(vector1, scalar);
    // ベクトルの足し算を実行
    Vector2D result = addVectors(vector1, scaledVector);
    // 結果を表示
    printf("Resultant Vector: (%d, %d)\n", result.x, result.y);
    return 0;
}

Resultant Vector: (4, 8)

このプログラムでは、vector1をスカラー倍した結果をscaledVectorに格納し、それをvector1と足し合わせています。

スカラー倍と足し算を組み合わせることで、ベクトルの拡大や縮小を伴う演算が可能です。

ベクトルの足し算を用いた物理シミュレーション

ベクトルの足し算は、物理シミュレーションにおいて重要な役割を果たします。

例えば、物体の位置や速度をベクトルで表現し、時間の経過とともに更新することで、物体の運動をシミュレートできます。

#include <stdio.h>
typedef struct {
    float x;
    float y;
} Vector2D;
Vector2D addVectors(Vector2D v1, Vector2D v2) {
    Vector2D result;
    result.x = v1.x + v2.x;
    result.y = v1.y + v2.y;
    return result;
}
int main() {
    Vector2D position = {0.0f, 0.0f}; // 初期位置
    Vector2D velocity = {1.0f, 1.5f}; // 速度ベクトル
    float timeStep = 1.0f; // 時間のステップ
    // 物体の位置を更新
    for (int i = 0; i < 5; i++) {
        position = addVectors(position, velocity);
        printf("Time %d: Position = (%.2f, %.2f)\n", i + 1, position.x, position.y);
    }
    return 0;
}

Time 1: Position = (1.00, 1.50)
Time 2: Position = (2.00, 3.00)
Time 3: Position = (3.00, 4.50)
Time 4: Position = (4.00, 6.00)
Time 5: Position = (5.00, 7.50)

このプログラムでは、物体の初期位置と速度をベクトルで表現し、時間のステップごとに位置を更新しています。

ベクトルの足し算を用いることで、物体の運動を簡潔にシミュレートすることができます。

ベクトル演算の最適化

ベクトル演算の最適化は、計算速度を向上させるために重要です。

特に大規模なデータセットを扱う場合、効率的な演算が求められます。

ここでは、いくつかの最適化手法を紹介します。

ループアンローリングによる最適化

ループアンローリングは、ループの繰り返し回数を減らすことで、ループのオーバーヘッドを削減する手法です。

これにより、処理速度が向上することがあります。

#include <stdio.h>
#define VECTOR_SIZE 8
void addVectors(int vector1[], int vector2[], int result[]) {
    for (int i = 0; i < VECTOR_SIZE; i += 2) {
        result[i] = vector1[i] + vector2[i];
        result[i + 1] = vector1[i + 1] + vector2[i + 1];
    }
}
int main() {
    int vector1[VECTOR_SIZE] = {1, 2, 3, 4, 5, 6, 7, 8};
    int vector2[VECTOR_SIZE] = {8, 7, 6, 5, 4, 3, 2, 1};
    int result[VECTOR_SIZE];
    // ベクトルの足し算を実行
    addVectors(vector1, vector2, result);
    // 結果を表示
    printf("Resultant Vector: ");
    for (int i = 0; i < VECTOR_SIZE; i++) {
        printf("%d ", result[i]);
    }
    printf("\n");
    return 0;
}

Resultant Vector: 9 9 9 9 9 9 9 9

このプログラムでは、ループアンローリングを用いて、1回のループで2つの要素を処理しています。

これにより、ループのオーバーヘッドを削減し、処理速度を向上させています。

SIMD命令を用いた最適化

SIMD(Single Instruction, Multiple Data)命令は、1つの命令で複数のデータを同時に処理することができるため、ベクトル演算の高速化に有効です。

以下は、SIMD命令を用いた最適化の例です。

#include <stdio.h>
#include <immintrin.h> // SIMD命令を使用するためのヘッダー
#define VECTOR_SIZE 8
void addVectorsSIMD(int vector1[], int vector2[], int result[]) {
    __m256i vec1 = _mm256_loadu_si256((__m256i*)vector1);
    __m256i vec2 = _mm256_loadu_si256((__m256i*)vector2);
    __m256i vecResult = _mm256_add_epi32(vec1, vec2);
    _mm256_storeu_si256((__m256i*)result, vecResult);
}
int main() {
    int vector1[VECTOR_SIZE] = {1, 2, 3, 4, 5, 6, 7, 8};
    int vector2[VECTOR_SIZE] = {8, 7, 6, 5, 4, 3, 2, 1};
    int result[VECTOR_SIZE];
    // ベクトルの足し算を実行
    addVectorsSIMD(vector1, vector2, result);
    // 結果を表示
    printf("Resultant Vector: ");
    for (int i = 0; i < VECTOR_SIZE; i++) {
        printf("%d ", result[i]);
    }
    printf("\n");
    return 0;
}

Resultant Vector: 9 9 9 9 9 9 9 9

このプログラムでは、SIMD命令を使用してベクトルの足し算を行っています。

_mm256_add_epi32は、256ビットの整数ベクトルを要素ごとに加算する命令です。

これにより、複数のデータを同時に処理し、演算を高速化しています。

メモリアクセスの最適化

メモリアクセスの最適化は、キャッシュの効率的な利用を促進し、処理速度を向上させるために重要です。

データの配置やアクセスパターンを工夫することで、キャッシュミスを減らすことができます。

#include <stdio.h>
#define VECTOR_SIZE 8
void addVectors(int vector1[], int vector2[], int result[]) {
    for (int i = 0; i < VECTOR_SIZE; i++) {
        result[i] = vector1[i] + vector2[i];
    }
}
int main() {
    int vector1[VECTOR_SIZE] = {1, 2, 3, 4, 5, 6, 7, 8};
    int vector2[VECTOR_SIZE] = {8, 7, 6, 5, 4, 3, 2, 1};
    int result[VECTOR_SIZE];
    // ベクトルの足し算を実行
    addVectors(vector1, vector2, result);
    // 結果を表示
    printf("Resultant Vector: ");
    for (int i = 0; i < VECTOR_SIZE; i++) {
        printf("%d ", result[i]);
    }
    printf("\n");
    return 0;
}