Juan M. Bello-Rivas jmbr
jmbr
/ poor-mans-spy.cpp
Last active
August 29, 2015 14:08
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| /* | |
| * Example: spy("foobar.ppm", A); | |
| */ | |
| int spy(std::string const& ppmfile, arma::mat const& M) { | |
| ofstream ofs(ppmfile); | |
| if (!ofs.is_open()) | |
| return -1; | |
| ofs << "P3\n" << M.n_rows << " " << M.n_cols << "\n1\n"; |
jmbr
/ gist:667599
Created
November 8, 2010 11:18
Plotting a coarse-grained model of (part of) a protein with PLPlot. This library apparently doesn't feature interactive rotation capabilities (version 5.9.7).
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| /* gcc -Wall -std=c99 plplotline3.c -o plplotline3 -lplplotd */ | |
| #include <stddef.h> | |
| #include <plplot/plplot.h> | |
| int main(void) | |
| { | |
| const size_t num_atoms = 10; | |
| const double positions[][3] = {{ 13.935, 18.529, 29.843 }, { 13.088, 19.661, 26.283 }, { 12.726, 17.033, 23.612 }, { 12.179, 17.659, 19.887 }, { 10.253, 15.79, 17.221 }, { 11.082, 16.103, 13.475 }, { 8.009, 15.163, 11.389 }, { 8.628, 13.975, 7.913 }, { 5.213, 12.642, 6.966 }, { 3.589, 12.601, 3.497 }}; |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| /** | |
| * @file bitset.c | |
| * @brief Implementation of a bit set data structure. | |
| */ | |
| #include <stdlib.h> | |
| #include <assert.h> | |
| #include "bitset.h" |
jmbr
/ cross_product_and_triple_scalar_product.c
Created
November 8, 2010 18:53
Cross product and triple scalar product using GSL.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| #include <gsl/gsl_vector.h> | |
| #include <gsl/gsl_blas.h> | |
| void cross_product(const gsl_vector *u, const gsl_vector *v, gsl_vector *product) | |
| { | |
| double p1 = gsl_vector_get(u, 1)*gsl_vector_get(v, 2) | |
| - gsl_vector_get(u, 2)*gsl_vector_get(v, 1); | |
| double p2 = gsl_vector_get(u, 2)*gsl_vector_get(v, 0) | |
| - gsl_vector_get(u, 0)*gsl_vector_get(v, 2); |
jmbr
/ test_stack_allocation_of_vectors.c
Created
November 17, 2010 07:34
Use of stack allocated vectors with the GNU Scientific Library.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| /* | |
| Compile with: | |
| gcc -Wall -std=c99 test_stack_allocation_of_vectors.c -o test_stack_allocation_of_vectors `gsl-config --libs` | |
| */ | |
| #include <stdio.h> | |
| #include <stddef.h> | |
| #include <string.h> | |
| #include <assert.h> |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| #include <cstdlib> | |
| #include <vector> | |
| #include <iostream> | |
| #include <ESBTL/default.h> | |
| #include <ESBTL/selected_atom_iterator.h> | |
| #include <ESBTL/atom_selectors.h> | |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| #include <cmath> | |
| inline unsigned next_power_of_two(unsigned x) { | |
| return unsigned(exp2f(ceilf(log2f(float(x))))); | |
| } |
jmbr
/ abort_unless.h
Created
August 30, 2013 17:56
Unlike the assert() macro, abort_unless() is not disabled by defining the preprocessor macro NDEBUG.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| #define abort_unless(expr) do { \ | |
| if (!(expr)) { \ | |
| fprintf(stderr, "%s:%u (%s): Assertion `%s' failed.\n", \ | |
| __FILE__, __LINE__, __func__, __STRING(expr)); \ | |
| fflush(stderr); \ | |
| abort(); \ | |
| } \ | |
| } while (0) |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| template <typename T = double> | |
| vector<T> linspace(T a, T b, size_t N) { | |
| T h = (b - a) / static_cast<T>(N-1); | |
| vector<T> xs(N); | |
| typename vector<T>::iterator x; | |
| T val; | |
| for (x = xs.begin(), val = a; x != xs.end(); ++x, val += h) | |
| *x = val; | |
| return xs; | |
| } |
jmbr
/ unrolled.s
Created
March 16, 2018 20:32
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| .text | |
| .align 4,0x90 | |
| .globl __Z21matrix_vector_productPKdPd | |
| __Z21matrix_vector_productPKdPd: | |
| LFB18: | |
| movsd (%rdi), %xmm0 | |
| movsd 8(%rdi), %xmm4 | |
| movapd %xmm0, %xmm5 | |
| movsd 16(%rdi), %xmm3 | |
| mulsd %xmm0, %xmm5 |
OlderNewer