// MPI_Broadcast - sends data from root to all processes
MPI_Comm_rank(MPI_COMM_WORLD, &rank);
MPI_Comm_size(MPI_COMM_WORLD, &size);

double data[1000];
if (rank == 0) {
    for (int i = 0; i < 1000; i++) data[i] = i;  // Root initializes
}
// All processes receive the data
MPI_Bcast(data, 1000, MPI_DOUBLE, 0, MPI_COMM_WORLD);

#pragma omp parallel for num_threads(8) schedule(static)
for (int i = 0; i < N; i++) {
    // Loop iterations are distributed statically
    result[i] = compute(a[i], b[i]);
}

#pragma omp parallel sections
{
    #pragma omp section
    { task1(); }  // Thread 1 executes
    #pragma omp section
    { task2(); }  // Thread 2 executes
}

// Initialize MPI
MPI_Init(&argc, &argv);
MPI_Comm_rank(MPI_COMM_WORLD, &rank);
MPI_Comm_size(MPI_COMM_WORLD, &size);

// Set OpenMP threads per MPI rank
omp_set_num_threads(4);

#pragma omp parallel for
for (int i = 0; i < local_n; i++) {
    local_data[i] = compute(i);  // Each rank computes its portion
}

// MPI communication across nodes
MPI_Allreduce(local_data, global_data, local_n, 
              MPI_DOUBLE, MPI_SUM, MPI_COMM_WORLD);