45 eigen_assert(dimension == 1 || dimension == 2 &&
"Transformation dimension must be less than 3.");
47 DFTI_DESCRIPTOR_HANDLE
res =
nullptr;
50 "DftiCreateDescriptor failed.")
51 handl.
reset(
res, [](DFTI_DESCRIPTOR_HANDLE handle) { DftiFreeDescriptor(&handle); });
52 if (forward_domain == DFTI_REAL) {
54 RUN_OR_ASSERT(DftiSetValue(handl.get(), DFTI_CONJUGATE_EVEN_STORAGE, DFTI_COMPLEX_COMPLEX),
55 "DftiSetValue failed.")
58 RUN_OR_ASSERT(DftiCreateDescriptor(&
res, precision, DFTI_COMPLEX, dimension,
sizes),
"DftiCreateDescriptor failed.")
59 handl.
reset(
res, [](DFTI_DESCRIPTOR_HANDLE handle) { DftiFreeDescriptor(&handle); });
62 RUN_OR_ASSERT(DftiSetValue(handl.get(), DFTI_PLACEMENT, DFTI_NOT_INPLACE),
"DftiSetValue failed.")
63 RUN_OR_ASSERT(DftiCommitDescriptor(handl.
get()), "DftiCommitDescriptor failed.")
#define eigen_assert(x)
Definition: Macros.h:910
cout<< "Here is the matrix m:"<< endl<< m<< endl;Matrix< ptrdiff_t, 3, 1 > res
Definition: PartialRedux_count.cpp:3
std::vector< Array2i > sizes
Definition: dense_solvers.cpp:12
#define RUN_OR_ASSERT(EXPR, ERROR_MSG)
Definition: ei_imklfft_impl.h:20
void reset(Vector< double > &unknowns)
Reset unknowns to some nontrivial garbage.
Definition: spring_contact.cc:53
Container::iterator get(Container &c, Position position)
Definition: stdlist_overload.cpp:29