NormMaxQuantities.cpp File Reference

Compute the Linf norm of a vector containing the physical quantities divided by their characteristic value. More...

#include "Carmen.h"

Include dependency graph for NormMaxQuantities.cpp:

Functions
real	NormMaxQuantities (const Vector &V)
	Returns the Max-norm of the vector where every quantity is divided by its characteristic value. More...

Detailed Description

Compute the Linf norm of a vector containing the physical quantities divided by their characteristic value.

Author

Anna Karina Fontes Gomes

Date

January-2017

Function Documentation

real NormMaxQuantities

(

const Vector &

V

)

Returns the Max-norm of the vector where every quantity is divided by its characteristic value.

Parameters

V

Vector

Returns

double

{
    Vector W(QuantityNb);
    int AxisNo=1;
    real MomentumMax=0.;
    real MagMax=0.;



    W.setZero();

/*
    // Density
    W.setValue(1, V.value(1)/QuantityMax.value(1));

    // Momentum
    W.setValue(2 , V.value(2)/QuantityMax.value(2) );
    W.setValue(3 , V.value(3)/QuantityMax.value(3) );
    W.setValue(4 , V.value(4)/QuantityMax.value(4) );

    // Energy
    W.setValue(5 , V.value(5)/QuantityMax.value(5) );

    // psi
   // W.setValue(6, V.value(6)/QuantityMax.value(6));

    // Magnetic Field
    W.setValue(7 , V.value(7)/QuantityMax.value(7) );
    W.setValue(8 , V.value(8)/QuantityMax.value(8) );
    W.setValue(9 , V.value(9)/QuantityMax.value(9) );
*/

  // --- Compute Linf norm --

    W.setValue(1, (V.value(1))/QuantityMax.value(1));
    W.setValue(5, (V.value(5))/QuantityMax.value(5));
    //W.setValue(6, (V.value(6))/QuantityMax.value(6));


    for (AxisNo = 1; AxisNo <= Dimension; AxisNo++)
    {
        MomentumMax = Max( MomentumMax, QuantityMax.value(AxisNo+1) );
        MagMax      = Max( MagMax,      QuantityMax.value(AxisNo+6) );
        W.setValue(2, W.value(2) + V.value(AxisNo+1)*V.value(AxisNo+1));
        W.setValue(7, W.value(7) + V.value(AxisNo+6)*V.value(AxisNo+6));
    }

    if(Dimension==2){
        W.setValue(4, (V.value(4))/QuantityMax.value(4));
        W.setValue(9, (V.value(9))/QuantityMax.value(9));
    }

    W.setValue(2, sqrt(W.value(2))/MomentumMax);
    W.setValue(7, sqrt(W.value(7))/MagMax );

    if(IterationNo==0)return NMax(V);
    return NMax(W);
}