Home - News ( United States | United Kingdom | Italy | Germany ) - Football scores

Showing content from http://www.arrayfire.org/docs/pde_2swe_8cpp-example.htm below:

ArrayFire: pde/swe.cpp

#include <math.h>

#include <stdio.h>

#include <stdlib.h>

#include <string.h>

float mx = max * 0.5;

float mn = -max * 0.5;

return (a - mn) / (mx - mn);

}

static void swe(bool console) {

const unsigned Lx = 1600, nx = Lx + 1;

const unsigned Ly = 1600, ny = Ly + 1;

const float dx = Lx / (nx - 1);

const float dy = Ly / (ny - 1);

ZERO = constant(0, nx, ny);

um = ZERO, vm = ZERO;

unsigned io = (unsigned)floor(Lx / 6.0f), jo = (unsigned)floor(Ly / 6.0f),

k = 15;

x = tile(range(nx), 1, ny);

y = tile(range(

(1, ny), 1), nx, 1);

0.01f * exp((-((x - io) * (x - io) + (y - jo) * (y - jo))) / (k * k));

float m_eta = max<float>(etam);

float dt = 0.5;

float h_diff_kernel[] = {9.81f * (dt / dx), 0, -9.81f * (dt / dx)};

float h_lap_kernel[] = {0, 1, 0, 1, -4, 1, 0, 1, 0};

h_diff_kernel_arr(3, h_diff_kernel);

h_lap_kernel_arr(3, 3, h_lap_kernel);

if (!console) {

win =

(1536, 768,

);

}

unsigned iter = 0;

unsigned random_interval = 30;

(!win->

()) {

if (iter > 2000) {

etam = 0.01f *

((-((x - io) * (x - io) + (y - jo) * (y - jo))) /

(k * k));

m_eta = max<float>(etam);

eta = etam;

iter = 0;

}

if (iter % 100 == 0 || iter % 130 == 0 || iter % random_interval == 0) {

io = (unsigned)

(rand() % Lx),

jo = (

)

(rand() % Ly);

random_interval = rand() % 200 + 1;

eta += 0.01f *

((-((x - io) * (x - io) + (y - jo) * (y - jo))) /

(k * k));

}

etap = 2 * eta - etam + (2 * dt * dt) / (dx * dy) * e;

etam = eta;

eta = etap;

m_eta = max<float>(etam);

if (!console) {

(*win)(0, 1).setAxesLimits(0, hist_out.

(), 0,

max<float>(hist_out));

(*win)(0, 0).image(normalize(eta, m_eta));

(*win)(0, 1).hist(hist_out, 0, 1,

"Normalized Pressure Distribution");

"Pressure at left boundary");

(*win)(1, 1).plot(

"Gradients versus Magnitude at left boundary");

win->show();

} else

iter++;

}

}

int main(int argc, char* argv[]) {

int device = argc > 1 ? atoi(argv[1]) : 0;

bool console = argc > 2 ? argv[2][0] == '-' : false;

try {

printf("Simulation of shallow water equations\n");

swe(console);

fprintf(stderr,

, e.

());

throw;

}

return 0;

}

Window object to render af::arrays.

A multi dimensional data container.

dim4 dims() const

Get dimensions of the array.

dim_t elements() const

Get the total number of elements across all dimensions of the array.

Generic object that represents size and shape.

An ArrayFire exception class.

virtual const char * what() const

Returns an error message for the exception in a string format.

@ AF_COLORMAP_BLUE

Blue hue map.

AFAPI array exp(const array &in)

C++ Interface to evaluate the exponential.

AFAPI array floor(const array &in)

C++ Interface to floor numbers.

array::array_proxy col(int index)

Returns a reference to a col.

AFAPI array range(const dim4 &dims, const int seq_dim=-1, const dtype ty=f32)

C++ Interface to generate an array with [0, n-1] values along the seq_dim dimension and tiled across ...

array & eval(array &a)

Evaluate an expression (nonblocking).

AFAPI void setDevice(const int device)

Sets the current device.

bool close()

Check if window is marked for close.

void grid(const int rows, const int cols)

Setup grid layout for multiview mode in a window.

AFAPI array histogram(const array &in, const unsigned nbins, const double minval, const double maxval)

C++ Interface for histogram.

AFAPI array flat(const array &in)

C++ Interface to flatten an array.

AFAPI array convolve(const array &signal, const array &filter, const convMode mode=AF_CONV_DEFAULT, const convDomain domain=AF_CONV_AUTO)

C++ Interface for convolution any(one through three) dimensional signals.

RetroSearch is an open source project built by @garambo | Open a GitHub Issue

Search and Browse the WWW like it's 1997 | Search results from DuckDuckGo

HTML: 3.2 | Encoding: UTF-8 | Version: 0.7.4