Skip to main content
Search type Search syntax
Tags [tag]
Exact "words here"
Author user:1234
user:me (yours)
Score score:3 (3+)
score:0 (none)
Answers answers:3 (3+)
answers:0 (none)
isaccepted:yes
hasaccepted:no
inquestion:1234
Views views:250
Code code:"if (foo != bar)"
Sections title:apples
body:"apples oranges"
URL url:"*.example.com"
Saves in:saves
Status closed:yes
duplicate:no
migrated:no
wiki:no
Types is:question
is:answer
Exclude -[tag]
-apples
For more details on advanced search visit our help page
Results tagged with
Search options not deleted user 242802

Algorithms which solve mathematical problems by means of numerical approximation (as opposed to symbolic computation).

2 votes
2 answers
407 views

Implementing multidimensional integral for a custom function in C++

I am not an expert with C++, but I am trying to implement a 4-dimensional integral using GSL numerical integration approach. The code below shows the whole algorithm. Although it seems correct what I …
CafféSospeso's user avatar
0 votes

Implementing multidimensional integral for a custom function in C++

Following the suggestion of @G.Sliepen, I have edited the PDFfunction to reduce divisions and replace cos(x) with an expansion of two complex exponentials, such that: cos(x) = 1/2 * ( cexp(I * x) + ce …
CafféSospeso's user avatar
7 votes
7 answers
1k views

A probability distribution function, to be called repeatedly during numerical integration

I am trying to speed up as much a possible this function in C++. As I explained in another post, Implementing multidimensional integral for a custom function in C++, this function will be used inside …
CafféSospeso's user avatar
3 votes
1 answer
284 views

Implementation of Multidimensional numerical integration in C++ and R

I'm trying to perform a 4-dimensional numerical integration in R using a function I wrote in C++ code which is then sourced in R using the Rcpp package. Below there is my code: // [[Rcpp::depends(Rcpp …
CafféSospeso's user avatar
0 votes
1 answer
258 views

Parallelise numerical integration with OpenMP in C++

I would like to parallelise with OpenMP a one-dimensional integral using the integrate() function implemented in the Boost library. I found a rather obscure example (for my C++ skills) in the Boost w …
CafféSospeso's user avatar