#include <iostream>
#include <vector>
#include <cstdlib>
//#include <cblas.h> // Include the CBLAS header

extern "C" {
double ddot_(const int* N, const double *X, const int* incX,
                  const double *Y, const int* incY);
}

int main() {
    // Define the size of the vectors
    const int N = 5;

    // Initialize two vectors
    std::vector<double> A = {1.0, 2.0, 3.0, 4.0, 5.0};
    std::vector<double> B = {5.0, 4.0, 3.0, 2.0, 1.0};

    // Calculate the dot product using CBLAS
    int t = 1;
    double dot_product = ddot_(&N, A.data(), &t, B.data(), &t);

    // Print the result
    std::cout << "Dot product of A and B: " << dot_product << std::endl;

    return 0;
}
//This should return a value of 35