sma_memory.C File Reference

Go to the source code of this file.

Functions
void	uexternaldb_ (int *LOP, int *, Real [], Real *, int *KSTEP, int *KINC)

Function Documentation

uexternaldb_()

void uexternaldb_	(	int *	LOP,
		int *	,
		Real	[],
		Real *	,
		int *	KSTEP,
		int *	KINC
	)

Definition at line 21 of file sma_memory.C.

{
  const bool threaded = getnumthreads_() > 1;

  switch (*LOP)
  {
    // beginning of the simulation
    case 0:
    {
      Real * fa = threaded ? SMALocalFloatArrayCreate(7, 5, 1.1) : SMAFloatArrayCreate(7, 5, 1.1);
      Real fasum = 0.0;
      for (int i = 0; i < 5; ++i)
        fasum += fa[i];

      // test wrong handle error based on kstep setting
      Real * fa2 = threaded ? SMALocalFloatArrayAccess(*KSTEP == 1 ? 2 : 7)
                            : SMAFloatArrayAccess(*KSTEP == 1 ? 2 : 7);
      if (fa != fa2)
        mooseError("Mismatching pointers");

      int * ia = threaded ? SMALocalIntArrayCreate(67, 39, 9) : SMAIntArrayCreate(67, 39, 9);
      int iasum = 0;
      for (int i = 0; i < 39; ++i)
        iasum += ia[i];

      // test wrong handle error based on kstep setting
      int * ia2 = threaded ? SMALocalIntArrayAccess(*KSTEP == 2 ? 13 : 67)
                           : SMAIntArrayAccess(*KSTEP == 2 ? 13 : 67);
      if (ia != ia2)
        mooseError("Mismatching pointers");

      if (get_thread_id_() == 0)
        Moose::out << *LOP << "lop " << fasum << ' '
                   << (threaded ? SMALocalFloatArraySize(7) : SMAFloatArraySize(7)) << ' ' << iasum
                   << ' ' << (threaded ? SMALocalIntArraySize(67) : SMAIntArraySize(67)) << '\n';
      break;
    }

    // beginning of the timestep
    case 1:
    {
      Real * fa = threaded ? SMALocalFloatArrayAccess(7) : SMAFloatArrayAccess(7);
      for (int i = 0; i < (threaded ? SMALocalFloatArraySize(7) : SMAFloatArraySize(7)); ++i)
        fa[i] += i + *KINC;

      int * ia = threaded ? SMALocalIntArrayAccess(67) : SMAIntArrayAccess(67);
      for (int i = 0; i < (threaded ? SMALocalIntArraySize(67) : SMAIntArraySize(67)); ++i)
        ia[i] += i + *KINC;
      break;
    }

    // end of the simulation
    case 3:
    {
      Real * fa = threaded ? SMALocalFloatArrayAccess(7) : SMAFloatArrayAccess(7);
      Real fasum = 0.0;
      for (int i = 0; i < (threaded ? SMALocalFloatArraySize(7) : SMAFloatArraySize(7)); ++i)
        fasum += fa[i] + 0;
      if (threaded)
        SMALocalFloatArrayDelete(7);
      else
        SMAFloatArrayDelete(7);

      int iasum = 0;
      int * ia = threaded ? SMALocalIntArrayAccess(67) : SMAIntArrayAccess(67);
      for (int i = 0; i < (threaded ? SMALocalIntArraySize(67) : SMAIntArraySize(67)); ++i)
        iasum += ia[i];

      if (threaded)
        SMALocalIntArrayDelete(67);
      else
        SMAIntArrayDelete(67);

      if (get_thread_id_() == 0)
        Moose::out << *LOP << "lop " << fasum << ' ' << iasum << ' ' << '\n';
      break;
    }
  }
}