Carmen.h

Go to the documentation of this file.

/***************************************************************************
                          Carmen.h  -  description
                             -------------------
    begin                : Thu Jun 7 2001
    copyright            : (C) 2001 by Olivier Roussel & Alexei Tsigulin
    email                : roussel@ict.uni-karlsruhe.de, lpsoft@mail.ru
 ***************************************************************************/

/***************************************************************************
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU General Public License as published by  *
 *   the Free Software Foundation; either version 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 ***************************************************************************/

#include "carmen.pre"


#if defined PARMPI
#include <mpi.h>
#endif

#include "PreProcessor.h"


#include <stdio.h>
#include <stdlib.h>
#include <iostream>
#include <math.h>

using namespace std;


#include "Vector.h"
#include "Matrix.h"
#include "Timer.h"
#include "Parameters.h"
#include "PrintGrid.h"
#include "TimeAverageGrid.h"
#include "Cell.h"
#include "Node.h"
#include "FineMesh.h"


#define Max(x,y) (((x) > (y)) ? (x):(y))

#define Max3(x,y,z) (Max((x),Max((y),(z))))

#define Min(x,y) (((x) < (y)) ? (x):(y))

#define Min3(x,y,z) (Min((x),Min((y),(z))))

#define power2(x)  ((x)*(x))

#define power3(x)  ((x)*(x)*(x))

#define Abs(x) ( ((x) < 0)? -(x):(x) )

/*
______________________________________________________________________________________________

External functions
______________________________________________________________________________________________

*/

void AdaptTimeStep();

int BC(int i, int AxisNo, int N=(1<<ScaleNb));

int BoundaryRegion(const Vector& X);

void Backup(Node* Root);

void Backup(FineMesh* Root);

double CPUTimeRef(int iterations, int scales);

real ComputedTolerance(const int ScaleNo);

int DigitNumber(int arg);

int FileWrite(FILE *f, const char* format, real arg);

Vector Flux(Cell& Cell1, Cell& Cell2, Cell& Cell3, Cell& Cell4, int AxisNo);

void GetBoundaryCells(Cell& Cell1, Cell& Cell2, Cell& Cell3, Cell& Cell4,
            Cell& C1, Cell& C2, Cell& C3, Cell& C4, int AxisNo);

Vector InitAverage(real x, real y=0., real z=0.);

real InitResistivity(real x, real y=0., real z=0.);

void InitParameters();

void InitTimeStep();

void InitTree(Node* Root);

Vector Limiter(const Vector u, const Vector v);

real Limiter(const real x);

real MinAbs(const real a, const real b);

real NormMaxQuantities(const Vector& V);

void Performance(const char* FileName);

void PrintIntegral(const char* FileName);

void RefreshTree(Node* Root);

void Remesh(Node* Root);

Vector FluxX(const Vector& Avg);

Vector FluxY(const Vector& Avg);

Vector FluxZ(const Vector& Avg);

Vector SchemeHLL(const Cell& Cell1, const Cell& Cell2, const Cell& Cell3, const Cell& Cell4, const int AxisNo);

Vector SchemeHLLD(const Cell& Cell1, const Cell& Cell2, const Cell& Cell3, const Cell& Cell4, const int AxisNo);

Matrix stateUstar(const Vector& AvgL, const Vector& AvgR, const real prel, const real prer, real& slopeLeft, real& slopeRight,real& slopeM, real& slopeLeftStar, real& slopeRightStar, int AxisNo);

void fluxCorrection(Vector& Flux,const Vector& AvgL, const Vector& AvgR, int AxisNo);

void ShowTime(Timer arg);

int Sign(const real a);

Vector ArtificialViscosity(const Vector& Cell1, const Vector& Cell2, real dx, int AxisNo);

Vector ResistiveTerms(Cell& Cell1, Cell& Cell2, Cell& Cell3, Cell& Cell4, int AxisNo);

Vector Source(Cell& UserCell);

real Step(real x);

void TimeEvolution(FineMesh* Root);

void TimeEvolution(Node* Root);

void View(FineMesh* Root, const char* AverageFileName);

void View(Node* Root, const char* TreeFileName, const char* MeshFileName, const char* AverageFileName);

void ViewEvery(FineMesh* Root, int arg);

void ViewEvery(Node* Root, int arg);

void ViewIteration(FineMesh* Root);

void ViewIteration(Node* Root);


//Parallel
void CreateMPIType(FineMesh *Root);

void CPUExchange(FineMesh *Root,int);

void FreeMPIType();

void CreateMPITopology();

void CreateMPILinks();

void ReduceIntegralValues();