Потоки, увеличить скорость работы программы

using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
using System.IO;
using System.Threading;
 
namespace Ekz
{
    class 
    class Program
    {
        static int N;//глобальная переменная для хранения количества точек
        static int countThreads, countThreadsEnd;//глобальная переменная для хранения количества потоков и количества завершенных потоков
        static double[,] matrix;//глобальная переменная для хранения координат точек
        static int count;//глобальная переменная для хранения количества треугольников
        static double S;//глобальная переменная для хранения заданной площади
        static double E;//глобальная переменная для хранения заданной точности
        static Semaphore hSemaphore;//глобальный семафор
 
        static void Main(string[] args)
        {
            int i;
            count = 0;
            S = Convert.ToDouble(Console.ReadLine());//Считывание из файла значения площади
            E = Convert.ToDouble(Console.ReadLine());//Считывание из файла значения точности
            N = Convert.ToInt32(Console.ReadLine());//Считывание из файла значения количества точек
            matrix = new double[N, 2];//создание матрицы, для считывания в нее координат точек
            string s = "";//объявление строки
            for (i = 0; i < N; i++)
            {
                s = Console.ReadLine();//считывание координат одной точки в строку
                matrix[i, 0] = Convert.ToDouble(s.Split(' ')[0]);//занесение первой координаты из строки в матрицу
                matrix[i, 1] = Convert.ToDouble(s.Split(' ')[1]);//занесение второй координаты из строки в матрицу
            }
        
 
            hSemaphore = new Semaphore(0, 1);//создание семафора
            countThreads = 10;//задание количества потоков
            countThreadsEnd = countThreads;//задание количества завершенных потоков
 
            for (i = 0; i < countThreads; i++)
            {
                //WaitCallback start = delegate { Program.Ckeck(k); };//делегат для передачи в метод потока параметров
                object m = i;
                ThreadPool.QueueUserWorkItem(Ckeck, m);//добавление потока в пул потоков
            }
            hSemaphore.WaitOne();//ожидание, пока освободится семафор методом Release
            Console.WriteLine(count);//вывод количества треугольников подходящих под условие
            Console.ReadKey();
        }
 
        static void Ckeck(object ind)//Метод выполняемый в потоке
        {
            int index=(int) ind;
            int dx = N / countThreads;//интервал точек, который будет проверяться в потоке
            int a = dx * index;//начало интервала
            int b = dx * (index + 1);//конец интервала
 
            if (index == countThreads - 1)//если это последний поток
                b = N;//то конец интервала будет равняться оличеству точек
            int i, j, k;
            for (i = a; i < b; i++)
                for (j = i + 1; j < N - 1; j++)
                    for (k = j + 1; k < N; k++)
                    {
                        double ss = triangleArea(matrix[i, 0], matrix[i, 1], matrix[j, 0], matrix[j, 1], matrix[k, 0], matrix[k, 1]);//определение площади треугольник
                        if (Math.Abs(S - ss) <= E)//проверка, если площадь найденного треугольника равна заданной площади S с заданной точностью E
                            Interlocked.Increment(ref count);//увеличение значения переменной count на еденицу, с помощью атомарной операции
                    }
            int l = Interlocked.Decrement(ref countThreadsEnd);//уменьшение значения переменной количества завершенных потоков countThreadsEnd на еденицу, с помощью атомарной операции
            if (l == 0)//если это был последний поток
                hSemaphore.Release();//то освобождаем семафор
        }
 
        static double triangleArea(double x1, double y1, double x2, double y2, double x3, double y3)//метод для определения площади треугольника
        {
            double a, b, c, p;
            c = Math.Sqrt(Math.Pow(y1 - y2, 2) + Math.Pow(x1 - x2, 2));//длина стороны c
            a = Math.Sqrt(Math.Pow(y2 - y3, 2) + Math.Pow(x2 - x3, 2));//длина стороны a
            b = Math.Sqrt(Math.Pow(y1 - y3, 2) + Math.Pow(x1 - x3, 2));//длина стороны b
 
            p = (a + b + c) / 2;//определение полупериметра
            return (Math.Sqrt(p * (p - a) * (p - b) * (p - c)));//возврат значения площади
        }
    }
}