Foadsf

OS := $(shell uname)
OPTIONS:=

ifeq ($(OS),Darwin)
	OPTIONS += -framework OpenCL
else
	OPTIONS += -l OpenCL
endif

all: main.c

Foadsf

cmake_minimum_required(VERSION 3.1)

project(OpenCL_Example)

find_package(OpenCL REQUIRED)
include_directories(${OpenCL_INCLUDE_DIRS})
link_directories(${OpenCL_LIBRARY})

add_executable(main main.c)
target_include_directories(main PUBLIC ${CMAKE_CURRENT_SOURCE_DIR})

Foadsf

#!/bin/sh -f

export CC="mpicc -fPIC"
export CXX="mpicxx"
export FC="mpif90"
export F77="mpif90"

export ELMER_INSTALL="/some/random/folder/elmerfem/install/"
export ELMER_HOME="/some/random/folder/elmerfem/install/"
export MPI_HOME="/usr/local/Cellar/open-mpi/3.1.1"

Foadsf

from functools import reduce
from operator import mul
import sympy as sp
import numpy as np

# Both of these function will be vectorized via np.vectorize.
# This will not speed up the computation at all (behind np.vectorize, there still are for loops)
# but this avoid nested for loop, making that easier to read (and deal with arbitrary number of dims).

def build_X(a, xs):

Foadsf

import numpy as np
from copy import deepcopy

def list_dim(a):
    if not type(a) == list:
        return []
    return [len(a)] + list_dim(a[0])

class Monomial:
    def __init__(self, coefficient, powers):

Foadsf

Foadsf

import numpy as np

def crop(A,D1,D2):
    return A[tuple(slice(D1[i], D2[i]) for i in range(A.ndim))]

def sumall(A):
    sum1=A
    for k in range(A.ndim):
        sum1 = np.sum(sum1,axis=0)
    return sum1

Foadsf

Foadsf

//following this: https://openfoamwiki.net/index.php/Main_ContribExamples/AxiSymmetric

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  5                                     |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile

Foadsf

// test blockMeshDict file
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  5                                     |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{