highs-sys 1.14.2

Rust binding for the HiGHS linear programming solver. See http://highs.dev.
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88

// Example 2. Application Using C and cuBLAS: 0-based indexing
//-----------------------------------------------------------
#include <cuda_runtime.h>
#include <math.h>
#include <stdio.h>
#include <stdlib.h>

#include <iostream>

#include "cublas_v2.h"
#define N 5

int main(void) {
  cudaError_t cudaStat;
  cublasStatus_t stat;
  cublasHandle_t handle;
  int i;
  float* devPtrA;
  float* a = 0;
  a = (float*)malloc(N * sizeof(*a));
  if (!a) {
    printf("host memory allocation failed");
    return EXIT_FAILURE;
  }

  for (i = 0; i < N; i++) a[i] = (float)(i * N);
  for (i = 0; i < N; i++) printf("%7.0f", a[(i)]);
  std::cout << std::endl;

  cudaStat = cudaMalloc((void**)&devPtrA, N * sizeof(*a));
  if (cudaStat != cudaSuccess) {
    printf("device memory allocation failed");
    free(a);
    return EXIT_FAILURE;
  }
  stat = cublasCreate(&handle);
  if (stat != CUBLAS_STATUS_SUCCESS) {
    printf("CUBLAS initialization failed\n");
    free(a);
    cudaFree(devPtrA);
    return EXIT_FAILURE;
  }

  stat = cublasSetVector(N, sizeof(*a), a, 1, devPtrA, 1);
  if (stat != CUBLAS_STATUS_SUCCESS) {
    printf("data download failed");
    free(a);
    cudaFree(devPtrA);
    cublasDestroy(handle);
    return EXIT_FAILURE;
  }

  float r;

  cublasStatus_t status = cublasSnrm2(handle, N, devPtrA, 1, &r);
  if (stat != CUBLAS_STATUS_SUCCESS) {
    printf("error at dnorm");
    free(a);
    cudaFree(devPtrA);
    cublasDestroy(handle);
    return EXIT_FAILURE;
  }

  std::cout << "Snorm: " << r << std::endl << std::endl;

  float* b;
  b = (float*)malloc(N * sizeof(*b));

  stat = cublasGetVector(N, sizeof(*a), devPtrA, 1, b, 1);
  if (stat != CUBLAS_STATUS_SUCCESS) {
    printf("data upload failed");
    free(a);
    cudaFree(devPtrA);
    cublasDestroy(handle);
    return EXIT_FAILURE;
  }

  cudaFree(devPtrA);
  cublasDestroy(handle);

  for (i = 0; i < N; i++) {
    printf("%7.0f", b[(i)]);
  }
  printf("\n");

  free(a);
  return EXIT_SUCCESS;
}