Найти расстояние между отрезками

#include <cstdlib>
#include <iostream>
#include <cmath>
#include <algorithm>
 
/** 
 *  1. реализован основной алгоритм типа метода секущих Ньютона с минимальным интерфейсом для ввода параметров
 *  2. для удобства восприятия код сконцентрирован в одном файле
 */
using namespace std;
 
/** точка */
class Point {
public:
 
    Point() {
        x = 0;
        y = 0;
        z = 0;
    }
 
    Point(double inX, double inY, double inZ) {
        x = inX;
        y = inY;
        z = inZ;
    }
 
    void init(double inX, double inY, double inZ) {
        x = inX;
        y = inY;
        z = inZ;
 
    }
    double x;
    double y;
    double z;
};
 
/** отрезок */
class Line {
public:
    Point* A; //начало отрезка
    Point* B; //конец отрезка
    Point* M; //середина отрезка
 
    ~Line() {
        delete A;
        delete B;
        delete M;
    }
};
 
/** паттерн фабрика */
class InitFabriq {
public:
 
    static Point* createPoint(double x, double y, double z) {
        Point* p = new Point(x, y, z);
        return p;
    }
 
    static Line* createLine(Point* A, Point* B) {
        Point& refA = *A;
        Point& refB = *B;
        Point* M = InitFabriq::createPoint((refA.x + refB.x) / 2, (refA.y + refB.y) / 2, (refA.z + refB.z) / 2);
        Line* l = new Line();
        l->A = A;
        l->B = B;
        l->M = M;
        return l;
    }
 
    static Point* copyPoint(Point* src, Point* dst) {
        dst->x = src->x;
        dst->y = src->y;
        dst->z = src->z;
        return dst;
    }
 
    static Point* mediumPoint(Point* A, Point* B, Point* M) {
        Point& refA = *A;
        Point& refB = *B;
        M->x = (refA.x + refB.x) / 2;
        M->y = (refA.y + refB.y) / 2;
        M->z = (refA.z + refB.z) / 2;
        return M;
    }
};
 
/** нечто типа метода секущих Ньютона - сложность log(n) */
class Algoritm {
public:
 
    /** расстояние */
    static double distance(Point& A, Point& B) {
        return sqrt((A.x - B.x) *(A.x - B.x) + (A.y - B.y)*(A.y - B.y)+(A.z - B.z)*(A.z - B.z));
    }
 
    /** основной алгоритм */
    static double MediumLogScan(Line* line1, Line* line2, double inEps) {
 
        double dist[3];
 
        double eps = numeric_limits<double>::infinity();
        Point* focus = line1->M;
 
        int count = 0;
        while (eps > inEps) {
            dist[0] = Algoritm::distance(*focus, *line2->A);
            dist[1] = Algoritm::distance(*focus, *line2->M);
            dist[2] = Algoritm::distance(*focus, *line2->B);
            cout << "distances:" << dist[0] << "  " << dist[1] << "  " << dist[2] << endl;
            if (dist[0] > dist[2]) {
                InitFabriq::copyPoint(line2->M, line2->A);
                eps = abs(dist[1] - dist[0]);
            } else {
                InitFabriq::copyPoint(line2->M, line2->B);
                eps = abs(dist[1] - dist[2]);
            }
            line2->M = InitFabriq::mediumPoint(line2->A, line2->B, line2->M);
            count++;
        }
        cout << "distance=" << dist[1] << " iterations=" << count << endl;
        return dist[1];
    }
 
};
 
 
int main(int argc, char** argv) {
 
    cout << "Упрощенный интерфейс ввода параметров" << endl;
    cout << "Пример ввода точки пространства через пробел: 0.1 0.1 0.5" << endl;
    
    double a, b, c;
    cout << "Введите первую точку первого отрезка: " << endl;
    cin >> a >> b >> c;
    Point* A1 = InitFabriq::createPoint(a, b, c);
    cout << "Введите вторую точку первого отрезка: " << endl;
    cin >> a >> b >> c;
    Point* B1 = InitFabriq::createPoint(a, b, c);
    Line* line1 = InitFabriq::createLine(A1, B1);
 
    cout << "Введите первую точку второго отрезка: " << endl;
    cin >> a >> b >> c;
    Point* A2 = InitFabriq::createPoint(a, b, c);
    cout << "Введите вторую точку второго отрезка: " << endl;
    cin >> a >> b >> c;
    Point* B2 = InitFabriq::createPoint(a, b, c);
    Line* line2 = InitFabriq::createLine(A2, B2);
    
    cout << "Введите порядок точности(например такой 0.0001) " << endl;
    double eps;
    cin >> eps;
  
    Algoritm::MediumLogScan(line1, line2, eps);
    double result = Algoritm::MediumLogScan(line2, line1, eps);
    cout<<"Минимальное расстояние = "<<result<<endl;
 
 
    delete line1;
    delete line2;
 
    return 0;
}