VPF::Перемножение квадратных матриц с вещ. числами

spk57

Дата 20.11.2014, 01:30 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 3
Регистрация: 4.7.2014

Репутация: нет
Всего: нет

Что-то никак не могу найти ошибку:

Код


#define _CRT_SECURE_NO_WARNINGS
#include <iostream>
#include <fstream>
#include <string>
#include <iomanip>

using namespace std;

int main(int argc, char *argv[])
{
    setlocale(LC_ALL, "rus");
    if (argc < 3)
    {
        cout << "Не заданы параметры" << endl;
        return 1;
    }

    FILE *firstFile, *secondFile, *resultFile;
    firstFile = fopen(argv[1], "r");
    secondFile = fopen(argv[2], "r");
    resultFile = fopen("result.txt", "w");

    /*Определение размера матриц*/
    char testStr[100000];
    int size1 = 0, size2 = 0;
    fgets(testStr, 100000, firstFile);
    for (int i = 0; i < strlen(testStr); i++)
    {
        if (testStr[i] == ' ') size1++;
    }
    fgets(testStr, 100000, secondFile);
    for (int i = 0; i < strlen(testStr); i++)
    {
        if (testStr[i] == ' ') size2++;
    }

    if (size1 != size2)
    {
        cout << "Матрицы разных размеров" << endl;
        return 1;
    }

    /*Считывание первой матрицы*/
    //создание матрицы
    float**firstMatrix = new float*[size1];

    //считывание матрицы
    const char separator[] = " ";
    char *ptr = NULL;
    char str[100000];
    for (int i = 0; i < size1; ++i)
    {
        fgets(str, 100000, firstFile);
        firstMatrix[i] = new float[size1];
        ptr = strtok(str, separator);
        for (int j = 0; j < size1; j++)
        {
            if (ptr == NULL)
                continue;
            firstMatrix[i][j] = atof(ptr);
            ptr = strtok(NULL, separator);
        }
    }

    /*Считывание второй матрицы*/
    //создание матрицы
    float**secondMatrix = new float*[size1];

    //считывание матрицы
    ptr = nullptr;
    for (int i = 0; i < size1; ++i)
    {
        fgets(str, 100000, secondFile);
        secondMatrix[i] = new float[size1];
        ptr = strtok(str, separator);
        for (int j = 0; j < size1; j++)
        {
            if (ptr == NULL)
                continue;
            secondMatrix[i][j] = atof(ptr);
            ptr = strtok(NULL, separator);
        }
    }


    cout << "Матрицы считаны" << endl;

    //результурующая матрицы
    float **resultMatrix = new float*[size1];
    for (int i = 0; i < size1; ++i)
        resultMatrix[i] = new float[size1];


    /*Перемножение матриц и вывод результата в файл*/
    cout << "Перемножаю матрицы" << endl;
    clock_t t1 = clock();

    for (int i = 0; i < size1; i++)
    {
        for (int j = 0; j < size1; j++)
        {
            resultMatrix[i][j] = 0;
            for (int k = 0; k < size1; k++)
            {
                resultMatrix[i][j] += firstMatrix[i][k] * secondMatrix[k][j];
            }
            fprintf(resultFile, "%1.2f ", resultMatrix[i][j]);
        }
        fprintf(resultFile, "\n");
    }
    fclose(firstFile);
    fclose(secondFile);
    fclose(resultFile);

    cout << "Подсчет окончен. Результат сохранен в result.txt" << endl;
    cout << "Время выполнения: " << (float)(clock() - t1) / CLOCKS_PER_SEC << endl;

    return 0;

Матрица 1:
0,3 0,2 3,9 0,0 4,6 8,0 9,7 3,5 2,3 6,0
9,2 1,8 4,6 9,5 9,9 7,3 6,8 9,9 3,0 4,6
5,4 4,0 7,0 8,5 7,5 5,4 1,6 5,2 2,8 0,7
0,6 0,4 2,8 4,0 2,3 6,8 3,4 3,3 7,3 4,0
9,7 1,8 9,1 6,5 7,9 8,4 0,3 2,7 9,5 8,5
7,6 8,7 1,7 1,4 3,3 4,6 0,1 9,3 3,0 5,0
8,2 1,8 9,9 9,9 8,0 6,7 1,8 6,9 3,4 1,2
4,8 2,4 3,7 2,0 1,5 8,3 9,9 7,2 9,3 8,9
1,9 4,7 4,9 7,3 3,7 4,7 0,6 1,0 3,2 1,4
6,1 7,4 3,5 2,1 4,4 2,6 3,7 7,8 6,7 2,6

Матрица 2:
5,8 8,5 5,0 9,3 5,4 9,9 7,3 0,7 7,6 1,5
7,8 6,7 5,9 7,1 1,7 0,8 2,2 4,8 8,4 4,9
1,0 7,2 4,6 2,9 8,1 1,4 1,3 2,7 6,9 7,3
8,6 8,0 7,9 8,1 2,3 2,1 1,3 3,8 1,2 6,0
6,7 6,2 2,0 5,4 5,3 5,3 3,8 4,2 1,2 5,9
0,7 3,6 8,6 1,1 0,2 1,7 4,6 4,1 1,1 6,2
4,0 0,5 5,6 0,7 4,1 5,3 1,7 1,3 8,3 7,1
0,2 7,3 4,8 9,6 2,6 1,9 9,6 4,0 0,6 1,5
3,5 1,6 7,9 8,7 0,3 9,9 0,0 0,1 2,2 0,8
7,5 4,6 7,2 4,8 1,4 6,3 9,6 0,7 8,5 7,5

Результат (явно не то):
305,91 -1985369216,00 -1985369088,00 -1985369088,00 -1985369344,00 -1985369216,00 -1985369216,00 -1985369344,00 -1985369216,00 -1985369216,00
234,89 -302121216,00 -302121152,00 -302121216,00 -302121312,00 -302121248,00 -302121280,00 -302121280,00 -302121248,00 -302121184,00
181,75 -1726408192,00 -1726408192,00 -1726408192,00 -1726408192,00 -1726408192,00 -1726408192,00 -1726408192,00 -1726408192,00 -1726408192,00
382,91 -3668617472,00 -3668617472,00 -3668617472,00 -3668617472,00 -3668617472,00 -3668617472,00 -3668617472,00 -3668617472,00 -3668617472,00
205,26 -2158010112,00 -2158010112,00 -2158010112,00 -2158010368,00 -2158010112,00 -2158010368,00 -2158010368,00 -2158010112,00 -2158010112,00

Guinness

Дата 20.11.2014, 08:17 (ссылка)

(нет голосов)

Загрузка ...

Опытный

Профиль
Группа: Участник
Сообщений: 310
Регистрация: 21.6.2009
Где: Зеленоград

Репутация: 2
Всего: 10

Рекомендую прогнать свою программу на более вменяемых(легко проверяемых) тестовых данных. Возьмите две матрицы:

Цитата

1 1 1
1 1 1
1 1 1

Я думаю, многое сразу станет ясным. И на будущее, разбивайте всё на функции, так проще отлаживаться и проверять корректность работы программы. В данном случае, можно выделить по крайней мере 3:
- подсчёт размера матрицы, в котором у Вас собственно и ошибка;
- заполнение матрицы;
- перемножение матриц.

1 Пользователей читают эту тему (1 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема \| C/C++: Для новичков \| Следующая тема »