Mercury13

Ошибка у тебя в scanf.

scanf("%lf",&x);

%lg тоже работает — по крайней мере, на MinGW.

Мало известно, но у scanf немного не такой формат подстановок, чем у printf.

void readResponses(string address, string (&rResponse)[15][4]) {}

Двухмерный массив — НЕ указатели на указатель, это тебе не Ява. Это массив одномерных массивов.

Нормального способа передачи массивов фиксированной длины в Си нет (только через ук-ль на первый элемент), в Си++ — через ссылку.

Первое возвращает по ссылке, второе по значению. Первое просто возвращает адрес, второе делает копию. Как определить, что нам нужно?
Где у нас результат: в локальной переменной или нет?

Перед нами операция «унарный плюс», которая, скорее всего, ничего не делает и return *this. Этот самый *this точно не локальный. Потому возвращаем по ссылке, String&.

А теперь возьмём бинарное сложение. Мы создаём новую строку (как локальную переменную, разумеется), заносим в неё сумму… а чтобы она не исчезла раньше, чем её подберёт вызывавшая подпрограмма, вернём по значению (String).

Для виртуальных функций обычно перестраховываются и возвращают по значению (String). Чхают на копирование, зато если вдруг потомку потребуется вернуть что-то локальное, он без проблем вернёт.

class Father {
private:
  String fName;
public:
  virtual String& name() const { return fName; }  // пока всё нормально, но…
};

class Son : public Father {
public:
  String& name() const
     { return Father::name() + "'s son"; }   // ошибка, возвращаем локальную переменную!
};

Как я понял, твоя матрица 3×3 — это однородные координаты в 2D? Я бы поступил так.
1. Убедиться, что элементы 3-1 и 3-2 нули (иначе — это не аффинное преобразование).
2. Элемент 3-3 превратить в единицу, соответственно увеличив остальные (на что — читай, что такое однородные координаты).
3. Элементы 1-3 и 2-3 — перенос. Отрежем их, получается матрица 2×2.
4. То, что осталось, должно быть вида (c, s), (-s, c). Если с какой-то погрешностью это не так и 2-норма строк не единица (тоже с какой-то погрешностью) — это не поворот (т.е. может быть масштабирование или наклон). Остаётся взять atan2(c, s) — получается угол.

Качество кода — на уровне студенческой лабораторной. Причём очень хреновой.

Программа решает четыре системы.
y = f0(t)
y' = f1(t, y)
y'' = f2(t, y, y')
y''' = f3(t, y, y'')

1. Если уж говорить о методе Эйлера — в программе ошибка.
y[1] = y[1] + delt * y[2];
Тут нужно старое значение y[2], а не новое.

2. Функцию eiler стоило бы обозвать eulerStep, заодно убрав лишние параметры.

3. Автор не знает такой концепции, как процедурный тип aka callback, поэтому работа с правой частью у него вышла вот через такую задницу.

4. В промышленном коде метод решения ОДУ я бы вынес в отдельную процедуру, с callback’ами на правую часть и на потребителя результата. Причём внешнее (доступное пользователю библиотеки) именование должно быть такое, чтобы её пользователь мог ею воспользоваться, даже если он не читал статью, по которой библиотеку писали. Так что смотрите, какие имена являются стандартизированными, понятными по любому учебнику, а какие стоит прояснить. Внутреннее именование — это уж думайте сами, научный код обычно сложен, и его постоянно приходится сличать со статьёй. Поэтому имена как в статье — для научного кода не WTF.

5. Не слишком удачно разделена ответственность между функцией правой части и функцией шага по Эйлеру. Это привело к ограничению — программа может решать только ОДУ y{n} = f(t, y, y', ..., y{n−1}).

Учи понятие «единица компиляции».

Два варианта.
a) Несколько единиц компиляции (основной вариант, когда идёт активная разработка проекта). Тогда в хедере каждую переменную отмечаешь модификатором extern — «переменная есть, но где-то в другом месте». В одной из единиц компиляции объявляешь переменные, уже без extern.

Насколько мне известно, защита от двойного включения для extern-переменных не важна, но, например, для классов или функций ой как пригодится.

///// unit.h /////
#ifndef UNIT_H
#define UNIT_H

extern int myVar;

#endif

///// unit.cpp /////
#include "unit.h"

int myVar;

///// main.cpp /////
#include "unit.h"

int main()
{
   myVar = 1;
}

б) Одна большущая единица компиляции (в такой вид часто преобразуют библиотеки перед выпуском). Тогда достаточно оградить каждый хедер от двойного включения, и всё в порядке.

///// unit.h /////
#ifndef UNIT_H
#define UNIT_H

int myVar;

#endif

///// main.cpp /////
#include "unit.h"

int main()
{
   myVar = 1;
}

ЗЫ. Когда я писал ответ, кода ещё не было. А тут он взял и пришёл. Для переменных (SYMBOL_WIDTH, LINE_HEIGHT) решение то самое. А для констант вариантов несколько.

Путь сишный.
#define LEFTPANEL 200
ВОЗМОЖНО для литеральных констант, и особенно классно работает с текстовыми строками (например, "[" TEXT "]").
НЕДОСТАТКИ: это препроцессор, и при пересечении имён будет нехорошо.

Путь перечислимый.

enum {
  LEFTPANEL = 200
};

ВОЗМОЖНО для констант, влезающих в int.
НЕДОСТАТКИ: возможны приколы с кроссплатформенностью, если на отдельных платформах константа вылезет за грань int.

Путь переменный.

extern const int LEFTPANEL;
....
const int LEFTPANEL = 200;

ВОЗМОЖНО для любых констант.
НЕДОСТАТКИ: нет вычисления при компиляции. Только линкер знает, что это 200.

Путь C++11.
constexpr int LEFTPANEL = 200;
ВОЗМОЖНО для литеральных констант.
НЕДОСТАТКИ: только C++11.

Путь «на рывок»
static const int LEFTPANEL = 200;
ВОЗМОЖНО для всех литеральных типов, но при этом может создавать дублирующиеся переменные: так, в Borland будет дублироваться всё сверх int — и, соответственно, не допускать предкомпилированных заголовков, ибо создаётся код.
НЕДОСТАТКИ: семантика зависит от компилятора и типа, ради чего и сделали constexpr.