evsl_direct_f90.c

Go to the documentation of this file.

#include <stdio.h>
#include <stdlib.h>
#include <stdint.h>
#include <string.h>
#include "def.h"
#include "struct.h"
#include "internal_proto.h"
#include "evsl_direct.h"
void EVSLFORT(setup_bsol_direct)(uintptr_t *Bf90, 
                                 uintptr_t *Bsoldataf90) {
  /* cast csr pointer of B */
  csrMat *B = (csrMat *) (*Bf90);
  void *Bsol;
  /* setup B sol and LT sol*/
  SetupBSolDirect(B, &Bsol);
  SetBSol(BSolDirect, Bsol);
  SetLTSol(LTSolDirect, Bsol); 
  /* cast pointer for output */
  *Bsoldataf90 = (uintptr_t) Bsol;
}

void EVSLFORT(free_bsol_direct)(uintptr_t *Bsolf90) {
  /* cast pointer */
  void *Bsol = (void *) (*Bsolf90);
  FreeBSolDirectData(Bsol);
}

void EVSLFORT(setup_asigmabsol_direct)(uintptr_t *Af90,
                                       int *flagB,
                                       uintptr_t *Bf90,
                                       uintptr_t *ratf90, 
                                       uintptr_t *solshiftf90) {
  /* cast csr pointer of A */
  csrMat *A = (csrMat *) (*Af90);
  /* cast csr pointer of B */
  csrMat *B = (*flagB) ? (csrMat *) (*Bf90) : NULL;
  /* cast pointer */
  ratparams *rat = (ratparams *) (*ratf90);
  /* allocate and setup solshiftdata */
  void **solshiftdata = (void **) malloc(rat->num*sizeof(void *));
  SetupASIGMABSolDirect(A, B, rat->num, rat->zk, solshiftdata);
  /* cast pointer for output */
  *solshiftf90 = (uintptr_t) solshiftdata;
}

void EVSLFORT(set_asigmabsol_direct)(uintptr_t *ratf90, 
                                     uintptr_t *solshiftf90) {
  /* cast pointers */
  ratparams *rat = (ratparams *) (*ratf90);
  void **solshiftdata = (void **) (*solshiftf90);

  SetASigmaBSol(rat, NULL, ASIGMABSolDirect, solshiftdata);
}

void EVSLFORT(free_asigmabsol_direct)(uintptr_t *ratf90, 
                                      uintptr_t *solshiftf90) {
  /* cast pointers */
  ratparams *rat = (ratparams *) (*ratf90);
  void **solshiftdata = (void **) (*solshiftf90);
  FreeASIGMABSolDirect(rat->num, solshiftdata);
  free(solshiftdata);
}