|
1 |
cmake ../ -DCMAKE_C_COMPILER=icc -DCMAKE_CXX_COMPILER=icpc -DCMAKE_Fortran_COMPILER=ifort -DCMAKE_INSTALL_PREFIX=/packages/opencv/3.4.2-intel2018 -DWITH_TBB=ON -DWITH_V4L=OFF -DWITH_QT=OFF -DWITH_OPENGL=OFF -DWITH_OPENCL=OFF -DWITH_GTK=OFF -DBUILD_TIFF=ON -DBUILD_JAVA=OFF -DWITH_GPHOTO2=OFF -DCMAKE_BUILD_TYPE=Release -DWITH_OPENCL=OFF -DWITH_IPP=OFF -DENABLE_CXX11=ON -DBUILD_opencv_video=OFF -DBUILD_TESTS=OFF -DBUILD_opencv_ts=OFF |
Format wyjściowy biblioteki deal.ii
|
1 2 3 4 5 6 |
(0,0) 6.73636e-09 (0,1) 2.76187e-09 (0,2) -3.27721e-09 (0,3) -3.6236e-11 (0,4) 1.13877e-08 (0,5) -5.36419e-10 |
Matlab wymaga numeracji od 1 i linii ze spacjami zaraz przecinków:
|
1 2 3 4 5 6 |
1 1 6.73636e-09 1 2 2.76187e-09 1 3 -3.27721e-09 1 4 -3.6236e-11 1 5 1.13877e-08 1 6 -5.36419e-10 |
Skrypt dokonujący konwersji:
|
1 2 3 4 5 6 7 8 9 |
#!/usr/bin/python import sys import re with open(sys.argv[1]+".new", 'a+') as w: with open(sys.argv[1]) as file: for line in file: m = re.match(r'\((.*?),(.*?)\)(.*?)$', line) w.write('{0} {1} {2}\n'.format(int(m.group(1))+1, int(m.group(2))+1, m.group(3))) |
Uruchomienie:
|
1 |
./matrix_convert.py A.dat |
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 |
#!/bin/bash # prune include files one at a time, recompile, and put them back if it doesn't compile # arguments are list of files to check removeinclude() { file=$1 header=$2 perl -i -p -e 's+([ \t]*#include[ \t][ \t]*[\"\<]'$2'[\"\>])+//REMOVEINCLUDE $1+' $1 } replaceinclude() { file=$1 perl -i -p -e 's+//REMOVEINCLUDE ++' $1 } for file in $* do #echo $file includes=`grep "^[ \t]*#include" $file | awk '{print $2;}' | sed 's/[\"\<\>]//g'` echo $includes for i in $includes do touch $file # just to be sure it recompiles removeinclude $file $i if make -j10 >/dev/null 2>&1; then grep -v REMOVEINCLUDE $file > tmp && mv tmp $file echo removed $i from $file else replaceinclude $file echo $i was needed in $file fi done done |
Skrypt usuwa kolejne #include a następnie rekompiluje kod by sprawdzić czy nadal działa. Jako argument skrypt przyjmuje katalog z plikami do sprawdzenia.
Narzędzie do statycznego analizowania kodu C++.
|
1 |
./cppcheck /home/kbzowski/femafat_new/ --std=c++11 --report-progress --enable=all -I /packages/deal.ii/8.2pre-gcc482/include/ -I /packages/boost/1.55.0-gcc482/ --xml --xml-version=2 -i /home/kbzowski/femafat_new/CMakeFiles 2> /home/kbzowski/femafat_new/cppcheck.xml |
-I – Include directory -i – Ignore directory –xml – Wyjście jako xml (potrzebne do raportów HTML) –enable=all – Sprawdza co się da Generowanie raportów HTML:
|
1 |
cppcheck-htmlreport --file=/home/kbzowski/femafat_new/cppcheck.xml --report-dir=/home/kbzowski/femafat_new/Reports --source-dir=/home/kbzowski/femafat_new |
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 |
#include <math.h> #include <iostream> using namespace std; extern "C" void uhard(double *, double *, double *, double *, double *, double *, double *, double *, int *, int *, int *, int *, int *, int *, char *, int *, double *, int *, double *, double *, int*, double *); /** User subroutine UHARD: is called at all material calculation points of elements for which the material definition includes user-defined isotropic hardening or cyclic hardening for metal plasticity; - can be used to define a material's isotropic yield behavior; - can be used to define the size of the yield surface in a combined hardening model; - can include material behavior dependent on field variables or state variables; - and requires, when appropriate, that the values of the derivatives of the yield stress (or yield surface size in combined hardening models) be defined with respect to the strain, strain rate, and temperature. */ void uhard(double *syield, double *hard, double *eqplas, double *eqplasrt, double *time, double *dtime, double *temp, double *dtemp, int *noel, int *npt, int *layer, int *kspt, int *kstep, int *kinc, char *cmname, int *nstatev, double *statev, int *numfieldv, double *predef, double *dpred, int *numprops, double *props) { // USER CODING TO DEFINE: // SYIELD - (sigma_0) - Yield stress for isotropic plasticity. Yield surface size for combined hardening. // hard[0] - Variation of SYIELD with respect to the equivalent plastic strain // hard[1] - Variation of SYIELD with respect to the equivalent plastic strain rate, This quantity is not used with the combined hardening model. // hard[2] - Variation of SYIELD with respect to temperature, This quantity is required only in adiabatic, fully coupled temperature-displacement, and thermal-electrical-structural analyses. // statev@nstatev - Array containing the user - defined solution - dependent state variables at this point.These are supplied as values at the beginning of the increment or as values updated by other user subroutines and must be returned as values at the end of the increment. // VARIABLES PASSED IN FOR INFORMATION: // eqplas - Equivalent plastic strain // eqplasrt - Equivalent plastic strain rate // time[0] - Value of step time at the beginning of the current increment. // time[1] - Value of total time at the beginning of the current increment. // dtime - Time increment // temp - Temperature at the beginning of the increment. // dtemp - Increment of temperature. // noel - Element number. // npt - Integration point number. // layer - Layer number(for composite shells and layered solids). // kspt - Section point number within the current layer. // kstep - Step number. // kinc - Increment number. // cmname - User - specified material name, left justified. // nstatv - User - specified number of solution - dependent state variables associated with this material // numfieldv - Number of field variables. // predef@numfieldv - Array of interpolated values of predefined field variables at this material point at the start of the increment based on the values read in at the nodes(initial values at the beginning of the analysis and current values during the analysis). // dpred@numfieldv - Array of increments of predefined field variables at this material point for this increment; this includes any values updated by user subroutine USDFLD. // numprops - Number of hardening properties entered for this user - defined hardening definition. // props@numprops - Array of hardening properties entered for this user - defined hardening definition. double a = 800; double n = 0.2; const double zero = 0.0; double eps = 0; if(*eqplas == zero) { eps = 0.01; } // syield - Naprezenie - obliczone na podstawie odksztalcenia eqpla *syield = a * powf(eps, n); // hard[0] - wartosc pochodnej po odksztalceniu hard[0] = a * n * powf(eps, n - 1); // hard[1] - wartosc pochodnej po predkosci odksztalcenia // hard[2] - wartosc pochodnej po temperaturze return; } |
HPC Pack 2012 R2 MS-MPI Redistributable Package HPC Pack 2012 Client Utilities Redistributable Package HPC Pack 2012 SDK with Service Pack 1 Wyedytuj user-config.jam:
|
1 2 3 4 |
using mpi : : <include>"c:/Program Files/Microsoft MPI/Inc" <library-path>"c:/Program Files/Microsoft MPI/Lib/i386" <find-shared-library>msmpi ; |
string w postaci: (0,0), (0.01,0), (0.01,0.01) zamieniony na tablicę obiektów typu punkt (x,y)
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 |
void TokenizeVectorOfPoints(std::vector<Point<2>>& pointsVector){ boost::regex re("(\\((.*?),(.*?)\\).*?)+"); std::string s = _project->heatTransferProbePoints; boost::sregex_iterator it(s.begin(), s.end(), re); boost::sregex_iterator end; for (; it != end; ++it) { double x, y; try { x = boost::lexical_cast<double>((*it)[2]); y = boost::lexical_cast<double>((*it)[3]); } catch (boost::bad_lexical_cast const&) { x = 0; y = 0; } pointsVector.emplace_back(Point<2>(x,y)); } } |
|
1 2 3 |
void testFunction() __attribute__((optimize(0))) { //code } |
Zainstalować Intel Parallel Studio XE Do \tools\build\v2\user-config.jam dodać using mpi : : <include>”c:/Program Files (x86)/Intel/MPI/4.1.0.023/ia32/include” <library-path>”c:/Program Files (x86)/Intel/MPI/4.1.0.023/ia32/lib” <find-shared-library>impi ;
Do pliku: {BOOST_ROOT}\tools\build\v2\user-config.jam należy dopisać: using mpi : : <include>”c:/Program Files (x86)/MPICH2/include” <library-path>”c:/Program Files (x86)/MPICH2/lib” <find-shared-library>mpi <find-shared-library>mpe <find-shared-library>cxx ;
