Foadsf

OS := $(shell uname)
OPTIONS:=

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

all: main.c

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

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

//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

//following this: https://openfoamwiki.net/index.php/Main_ContribExamples/AxiSymmetric
//  /cygdrive/c/Users/sojoodif/surfdrive/sync/MyPhD/simulation/OpenFOAM/20180131_2

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

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
{