userset.cpp

Go to the documentation of this file.

#include "userset.h"
#include <cmath>
 
void SetU(int const nx, int const ny, double u[])
{
    int dx    = nx + 1;
    for (int j = 0; j <= ny; ++j)
    {
        int k = j * dx;
        for (int i = 0; i <= nx; ++i, ++k)
        {
            u[k] = 0.0;
        }
    }
}
 
void SetF(int const nx, int const ny, double f[])
{
    int dx    = nx + 1;
    for (int j = 0; j <= ny; ++j)
    {
        int k = j * dx;
        for (int i = 0; i <= nx; ++i, ++k)
        {
            f[k] = 1.0;
        }
    }
}
 
 
double FunctF(double const x , double const y)
{
// return  std::sin(3.14159*1*x)*std::sin(3.14159*1*y);
//  return 16.0*1024. ;
// return (double)1.0 ;
    return x * x * std::sin(2.5 * 3.14159 * y);
}
 
double FunctU(const double /* x */, double const /* y */)
{
    return 1.0 ;
}