Skip to content

Instantly share code, notes, and snippets.

@hemmer
hemmer / gist:6864729
Last active December 24, 2015 21:19
Mathematica factoring
Here is an example input:
In[89]:= AA = 1/(24 \[Eta] \[Tau]C \[Tau]Q) (-12 \[Eta] \[Tau]Q - 12 a^2 GC \[Tau]C \[Tau]Q - 3 \[Zeta] \[Tau]C \[Tau]Q + 5 \[Zeta] \[Xi] \[Tau]C \[Tau]Q + Sqrt[(48 \[Zeta] \[Eta] (-3 + 5 \[Xi]) \[Tau]C + (12 \[Eta] + (12 a^2 GC + \[Zeta] (3 - 5 \[Xi])) \[Tau]C)^2) \[Tau]Q^2])
This method doesn't work:
In[90]:= Solve[AA == 0, \[Zeta]]
During evaluation of In[90]:= $RecursionLimit::reclim: Recursion depth of 20 exceeded. >>
@hemmer
hemmer / Makefile
Last active January 31, 2020 18:24
Example of FFTW transforms with various indexing methods. See printDataFourier for examples.
program_NAME := fourier
CC = cc
program_C_SRCS := $(wildcard *.c) $(wildcard */*.c)
program_C_OBJS := ${program_C_SRCS:.c=.o}
program_OBJS := $(program_C_OBJS) $(program_CXX_OBJS)
program_INCLUDE_DIRS :=
program_LIBRARY_DIRS :=
program_LIBRARIES := pthread fftw3 m gsl gslcblas
program_FLAGS := -Wall -Wextra -g -std=c99 -Wshadow
#program_FLAGS := -Wall -Wextra -O3 -std=c99 -Wshadow