Scientific Computing

Add missing LaTeX fonts

Missing LaTeX fonts can be added via TeXLive or MikTeX. Using the symbolic fonts is as easy as:

\documentclass[a4paper,12pt]{article}

\usepackage{fontawesome5}

\begin{document}
  \faGithub This is a GitHub logo.
\end{document}

This may need to use XeLaTeX. We use fontawesome5 instead of obsolete fontawesome, which is version 4. FontAwesome5 is in TeXLive 2018.

TeXLive is popular across operating systems. Linux users can use system package managers to install groups of TeXLive packages. Advanced Linux TeXLive users can independently setup TeXLive to get the latest packages individually, to save hundreds of megabytes of install space.


For marvosym:

Windows/Linux tlmgr:

tlmgr install marvosym

Linux: apt:

apt install texlive-fonts-recommended

For fontawesome5:

FontAwesome gives popular emoji and website icons.

Windows/Linux: tlmgr:

tlmgr install fontawesome5

Linux: apt:

apt install texlive-fonts-extra

MikTeX is another LaTeX distro on Windows. MiKTeX will automatically install packages as needed. If that isn’t happening, check:

Start → MiKTeX Console → Settings → General → “Always install missing packages on the fly”

Otherwise, manually search for package name and install:

Start → MiKTeX Console → Packages

Install OpenCV in Python

The unofficial OpenCV PyPI wheels work with pip install methods:

pip install opencv-python

For ARM / Raspberry Pi:

pip install opencv-python

also works for certain ARM platforms like the Raspberry Pi.

OpenCV is trivial and fast to install on a Raspberry Pi via pip as described above.

For the latest extended functionality that hasn’t yet been incorporated into the core package, OpenCV including the Extra contributed modules may be obtained by:

pip install opencv-contrib-python

If ... is not supported on this platform error be sure it’s using desired Python install. One may have to manually specify the path for the pip command e.g.

$Env:SystemDrive\anaconda3\scripts\python -m pip install opencv-python

We made several test scripts to try out the OpenCV install. Compiling OpenCV yourself allows customizing and optimizing OpenCV for your computer (e.g. using GPU, TBB, OpenCL, etc.).

The conda install opencv and conda install -c conda-forge opencv methods for OpenCV continue to be BROKEN for video/image reading and display. Use pip install above instead.

Convert Periscope video for YouTube

If YouTube won’t accept a video upload, or the video never completes “Processing” on YouTube, try re-encoding the video with FFmpeg. Sometimes lossy conversion is necessary to achieve the YouTube recommended upload settings.

Periscope video downloads use MPEG TS container with .ts file extension. These .ts files can be played back in VLC or similar to confirm content. A lossless conversion to YouTube is possible with:

ffmpeg -i pscp.ts -bsf:a aac_adtstoasc -codec copy -max_muxing_queue_size 1000 pscp.mp4

The FFmpeg option -max_muxing_queue_size is an arbitrary parameter. The default queue size is sometimes not large enough:

Too many packets buffered for output stream 0:1.
  • Do not use .mkv extension for YouTube uploads, it will fail to process.
  • certain formats like AV1 or FFV1 will fail to finish processing the upload
  • after completing Uploading, video should start Processing in less than 5 minutes. If not, it will probably never actually work.

Uploads that don’t finish conversion processing after upload are not viewable on YouTube.

Generate vectors of datetime in Python

Generating a range of datetime data is a common data analysis and simulation task. Here we show examples of generating datetime vectors for Python datetime and numpy.datetime64

Python datetime

Python datetime is often used as timezone-naïve with UTC as the assumed timezone. This custom avoids ambiguities when working with Pandas and Numpy, which are foundational for Python data science.

Generate a range of Python datetime like:

from __future__ import annotation
from datetime import datetime, timedelta


def datetime_range(start: datetime, end: datetime, step: timedelta) -> list[datetime]:
    """like range() for datetime"""
    return [start + i * step for i in range((end - start) // step)]


dt = datetime_range(datetime(2019, 12, 1), datetime(2020, 4, 1), timedelta(days=1))

Numpy datetime64

Numpy datetime64 generates a range of times like:

dt = numpy.arange('2019-12-01', '2020-04-01', dtype='datetime64[D]')

Pandas has the date_range function to generate time vectors.

Matlab pcolor datetime plots

Matlab datetime works much like Python datetime, and is generally recommended. Matlab plotting functions generally support “datetime” class. It is sometimes necessary to manipulate Matlab plots involving datetime for the desired result.

dt = datetime(2019, 12, 1):datetime(2020, 3, 1);

dat = rand([100, length(dt)]);
y = 1:size(dat, 1);

pcolor(dt, y, dat)

If the plot doesn’t have the desired datetime axis formatting, try datetick.

datetick('x', 'yyyy-mm-dd', 'keepticks')

Generate vector of datetime. Python matplotlib supports datetime.datetime and numpy.datetime64 in most plot functions.

Fix Matlab graphics error after driver update

After updating an operating system graphics driver while Matlab was open, an error may occur upon plotting in Matlab like:

MATLAB has experienced a low-level graphics error, and may not have drawn correctly. Read about what you can do to prevent this issue at Resolving Low-Level Graphics Issues then restart MATLAB.

The solution to this issue typically is, as the message suggests, just restarting Matlab. A full reboot is not always necessary, but try that if the error persists.

Skip test with return code 77

Skipping a test with Meson or CMake by returning error code 77 is a de facto practice. Sometimes it is only feasible to know a test should be skipped by attempting to run that test.

Meson build system accepts return code 77 by default as a signal to skip the test. Configure CMake to skip return code:

add_test(NAME MyTest COMMAND mytest)  # arbitrary

set_property(TEST MyTest PROPERTY SKIP_RETURN_CODE 77)

Check temperature of Raspberry Pi CPU

The vcgencmd utility allows reading a few dozen hardware measurements on the Raspberry Pi boards.

CPU Temperature is checked by:

vcgencmd measure_temp

Typical temperatures in office environment, with case:

Pi Model heatsink usage temp [C]
2B yes light-moderate 40..45
4B no light-moderate 65..70
4B no YouTube 720p60 85

A red thermometer icon GPU-superimposed on the Raspberry Pi display output means the Raspberry Pi is overheating and is throttling the CPU and GPU to avoid self-destruction.

Raspberry Pi 0, 1, 2, 3 temperature thresholds:

CPU temp. [C] icon throttle
< 80 none none
80 - 85
half full thermometer red
CPU
> 85
full thermometer red
CPU & GPU

Raspberry Pi 4 temperature thresholds:

CPU temp. [C] icon throttle
80 - 85 none CPU: 1000 MHz
85 - 90
half full thermometer red
CPU
> 90
full thermometer red
CPU & GPU

The current Raspberry Pi CPU clock speed is obtained from

vcgencmd measure_clock arm

The output is in units of Hertz. The Raspberry Pi CPU clock speed and power consumption is dynamic:

clock speed [MHz] Raspi 2 Raspi 3 Raspi 3+ Raspi 4
idle 600 600 600 600
100% one or more cores 900 1200 1400 1500

Log temperature measurements with crontab -e. This can periodically log temperature and CPU frequency, e.g. add a line like:

@hourly  vcgencmd measure_temp | /usr/bin/logger

logger writes the measured parameters into the system log.


Related: Monitor Raspberry Pi DC input voltage

Create a blank/orphan Git branch

Many basic Git use cases involve a main branch with feature branches periodically merged into the main branch. For certain purposes, totally distinct branch without a common history can exist in the same Git repo. One of the most common uses of this is for documentation. For example, GitHub will build a website from the gh-pages branch.

Setup blank Git branch

Do NOT force push during this procedure, you may accidentally erase years of work!

This example assumes you want to create a gh-pages empty branch for documentation on GitHub, but will of course work for other purposes too.

From the repo directory create a blank Git branch:

git switch --orphan gh-pages

Remove the unneeded files from this branch (by default, all existing files are staged from the previous branch)

git rm --cached -r .

git clean -id

This leaves the .git/ directory, which should not be disturbed.

Copy documentation files

What happens next depends on if your documentation files were already added to another branch (tracked) or were not added to Git (untracked). Assume wanted files for the blank gh-pages branch are in docs/ on main branch.

copy over the files to gh-pages

git checkout main -- docs/

git commit -am "moved docs"

Upload files

  1. push the documentation

    git push -u origin gh-pages
  2. Enable the documentation builds from github.invalid/username/myrepo/settings → GitHub Pages under Source select gh-pages.

  3. In a few minutes, the webpages at username.github.io/myrepo/ should be active.

Once everything is working, the old docs/ folder isn’t needed.

Python f-string benchmarks

Python 3.6 f-strings have been shown to be the fastest string formatting method in microbenchmarks by Python core dev Raymond Hettinger in relative speed factors:

  • f-string: 1.0
  • concatenate string: 1.33
  • join sequence of strings: 1.73
  • %s formatting operator: 2.40
  • .format() method: 3.62
  • Template() method: 11.48

The reason for this speedy performance was described by Python core dev Serhiy Storchaka.