- GiNaC-list - ginac.de

Re: [GiNaC-list] problems with cstring and string
by Tania Robens 09 Feb '17

09 Feb '17

ok i will try this i did install cln locally and configured with something like export LD_LIBRARY_PATH=$LD_LIBRARY_PATH:$CLN/lib ./configure --prefix=/home/trobens/programs/ginac/install CLN_CFLAGS=-I$CLN/include/cln/ CLN_LIBS=-L$CLN/lib/ CPPFLAGS=-I$CLN/include/ LDFLAGS=-L$CLN/lib/ LIBS=-lcln mb this causes the problem ?? On Thu, 9 Feb 2017, Mario Prausa wrote: > > > On 09.02.2017 03:51, Tania Robens wrote: >> >> In file included from >> /home/trobens/programs/cln/install/include/cln/string.h:9:0, >> from >> /opt/rh/devtoolset-3/root/usr/include/c++/4.9.2/cstring:42, >> from class_info.h:27, >> from registrar.h:26, >> from basic.h:29, >> from ex.h:26, >> from expair.h:26, >> from expairseq.h:26, >> from add.h:26, >> from add.cpp:23: >> > > This looks suspicious. The ginac file "class_info.h" includes the > standard C++ library file "cstring" which includes the CLN file > "string.h". In the standard C libray there is also a file by the name > "string.h" which is what "cstring" wants to include. > > I think the problem is the order of your include pathes. The pathes to > the standard libraries should go before any other libraries (in that > case CLN). > > Mario > ------------------------------------------------------------------------- Dr. habil. Tania Robens Michigan State University Department of Physics and Astronomy 567 Wilson Road East Lansing, MI 48824 USA

3 3

problems with cstring and string
by Tania Robens 09 Feb '17

09 Feb '17

Hi I have tried to implement the newest CLN version (cln-1.3.4) and ginac version (ginac-1.7.2). i am using gcc (GCC) 4.9.2 20150212 (Red Hat 4.9.2-6) (but also had the same problem with newer compilers) uname -a renders Linux 2.6.32-642.11.1.el6.x86_64 #1 SMP Tue Nov 15 14:13:21 CST 2016 x86_64 x86_64 x86_64 GNU/Linux CLN works, but for ginac i get checking readline/readline.h usability... no checking readline/readline.h presence... yes configure: WARNING: readline/readline.h: present but cannot be compiled configure: WARNING: readline/readline.h: check for missing prerequisite headers? configure: WARNING: readline/readline.h: see the Autoconf documentation configure: WARNING: readline/readline.h: section "Present But Cannot Be Compiled" configure: WARNING: readline/readline.h: proceeding with the compiler's result configure: WARNING: ## ---------------------------------- ## configure: WARNING: ## Report this to ginac-list(a)ginac.de ## configure: WARNING: ## ---------------------------------- ## checking for readline/readline.h... no checking readline/history.h usability... no checking readline/history.h presence... yes configure: WARNING: readline/history.h: present but cannot be compiled configure: WARNING: readline/history.h: check for missing prerequisite headers? configure: WARNING: readline/history.h: see the Autoconf documentation configure: WARNING: readline/history.h: section "Present But Cannot Be Compiled" configure: WARNING: readline/history.h: proceeding with the compiler's result configure: WARNING: ## ---------------------------------- ## configure: WARNING: ## Report this to ginac-list(a)ginac.de ## configure: WARNING: ## ---------------------------------- ## checking for readline/history.h... no and, when i ignore this and do make anyways, i have In file included from /home/trobens/programs/cln/install/include/cln/string.h:9:0, from /opt/rh/devtoolset-3/root/usr/include/c++/4.9.2/cstring:42, from class_info.h:27, from registrar.h:26, from basic.h:29, from ex.h:26, from expair.h:26, from expairseq.h:26, from add.h:26, from add.cpp:23: /opt/rh/devtoolset-3/root/usr/include/c++/4.9.2/cstring:75:11: error: ‘::memchr’ has not been declared using ::memchr; ^ /opt/rh/devtoolset-3/root/usr/include/c++/4.9.2/cstring:76:11: error: ‘::memcmp’ has not been declared using ::memcmp; ...etc ... which then leads to a stop in the compilation. i have tried various ways around this, if someone has a quick solution please let me know the same mistake occurs with earlier ginac versions and, as stated, with newer gccs (up to 6.1) ?? Thanks Best Tania Robens ------------------------------------------------------------------------- Dr. habil. Tania Robens Michigan State University Department of Physics and Astronomy 567 Wilson Road East Lansing, MI 48824 USA

3 3

About non commutative transformations in GiNaC
by abpetrov 26 Jan '17

26 Jan '17

Hello, I wrote a program to do calculations with non commutative variables in GiNaC. But when I made attempt to apply noncommutative rule I got error: terminate called after throwing an instance of 'std::logic_error' what(): add::eval(): sum of non-commutative objects has non-zero numeric term Below is simple example of program with program output. How to change the class ncsymbol so that it was possible to apply rules like a*ap==ap*a+1? Best regards Petrov A.B. -------------------------------------------------------------------- Output for program: Starting: ./test2 e1: a*ap*a terminate called after throwing an instance of 'std::logic_error' what(): add::eval(): sum of non-commutative objects has non-zero numeric term *** Crashed with return code: 0 *** ------------------------------------------------------------------ Program: #include <iostream> using namespace std; #include<ginac/ginac.h> using namespace GiNaC; class ncsymbol : public symbol { GINAC_DECLARE_REGISTERED_CLASS(ncsymbol, symbol) public: ncsymbol(const string &s); protected: unsigned return_type() const { return return_types::noncommutative_composite; } }; GINAC_IMPLEMENT_REGISTERED_CLASS(ncsymbol, symbol) ncsymbol::ncsymbol() { } ncsymbol::ncsymbol(const string &s) : inherited(s) { } int ncsymbol::compare_same_type(const basic &other) const { GINAC_ASSERT(is_a<ncsymbol>(other)); const ncsymbol *o = static_cast<const ncsymbol *>(&other); if (serial==o->serial) return 0; return serial < o->serial ? -1 : 1; } void calc_example(){ ncsymbol a("a"),ap("ap"); ex e1 = a*ap*a; cout << "e1: " << e1 << endl; e1 = e1.subs(a*ap==ap*a+1); cout << "e1 subs: " << e1 << endl; } int main(int argc, char** argv) { calc_example(); return 0; }

2 4

Automatic transformation of expressions in GiNaC
by abpetrov 12 Jan '17

12 Jan '17

Hello, I wrote some program to do some calculations with GiNaC. But I can see program makes some unwanted transformations of expressions automatically. Below is simple example. Can I avoid that transformations? Wanted output is sqrt(m/2). Best regards Petrov A.B. -------------------------------------------------------------------- Output for program: Starting: ./test2 1/2*sqrt(2)*sqrt(m) *** Exited normally *** ------------------------------------------------------------------ Program: #include <iostream> using namespace std; #include<ginac/ginac.h> using namespace GiNaC; void calc_example(){ possymbol m("m"); ex e = sqrt(m/2); cout << e << endl; } int main(int argc, char** argv) { calc_example(); return 0; }

2 2

Release 1.7.2
by Richard B. Kreckel 10 Jan '17

10 Jan '17

Hi, A new GiNaC release is out and available. These are the changes: * Fix memory leak in operators returning relationals, introduced in 1.7.0. * Make macro GINAC_DECLARE_REGISTERED_CLASS() work outside namespace GiNaC. * Fix some build issues. Happy hacking -richy. -- Richard B. Kreckel <http://in.terlu.de/~kreckel/>

1 0

Memory leak in dynallocate for relational
by Thomas Luthe 09 Jan '17

09 Jan '17

Hi, valgrind reports a memory leak caused by dynallocate<GiNaC::relational, const GiNaC::ex&, const GiNaC::ex&, GiNaC::relational::operators> (basic.h:336). I have attached a minimal example as well as the corresponding valgrind output I get. The executable was produced with g++ 5.4 -O0 and GiNaC 1.7.1. Repeated use of operators == or != for GiNaC::ex does indeed lead to a noticeable increase of memory usage in applications, whereas ex::is_equal does not have the same effect. Cheers, Thomas

2 1

Simplification in GiNaC
by abpetrov 07 Jan '17

07 Jan '17

Hello, I recently wrote a simple program to learn some calculations with GiNaC. This works, But I got some expression which automatically don't become simple. The question is: What is the simplest way to get the expected answer - Ax here? Best regards Petrov A.B. -------------------------------------------------------------------- Output for program: Starting: ./test2 sqrt(m*x^2*A)*sqrt(m^(-1)*A) *** Exited normally *** ------------------------------------------------------------------ Program: #include <iostream> using namespace std; #include<ginac/ginac.h> using namespace GiNaC; void calc_example(){ realsymbol A("A"),m("m"),x("x"); ex e = sqrt(A*m*pow(x,2))*sqrt(A/m); cout << e << endl; } int main(int argc, char** argv) { calc_example(); return 0; }

2 2

Various syntax errors while compiling appln. with Ginac1.7.1
by harish sukhwani 21 Oct '16

21 Oct '16

Hi Everyone, I am new to GiNaC. In our research group. we have a tool written predominantly in C (and some c++ with GiNaC). I have inherited this recently, and I am trying to cross-compile using mingw32 on Ubuntu 14.04 with Win32 target. Note that the compiling application for local target with Ginac1.5 installed works fine. I compile the only single file using GiNaC using i686-w64-mingw32-g++ -c -o symbolic-diff.o symbolic-diff.cpp I see the following error (and a LOT more errors: In file included from /usr/i686-w64-mingw32/include/ginac/basic.h:27:0, from /usr/i686-w64-mingw32/include/ginac/ginac.h:28, from symbolic-diff.cpp:17: (line corresponding to #include <ginac/ginac.h>) /usr/i686-w64-mingw32/include/ginac/ptr.h:37:13: error: expected ';' at end of member declaration refcounted() noexcept : refcount(0) {} ^ /usr/i686-w64-mingw32/include/ginac/ptr.h:37:15: error: 'noexcept' does not name a type refcounted() noexcept : refcount(0) {} I am not sure where to start. ..... FYI, here are my installation steps for CLN1.3.3 and GiNaC 1.7.1 HOST=i686-w64-mingw32 CXX=${HOST}-g++ CC=${HOST}-gcc CPPFLAGS=-I/usr/i686-w64-mingw32/include CXXFLAGS="$CPPFLAGS -O3" CFLAGS=$CXXFLAGS LDFLAGS=-L/usr/i686-w64-mingw32/lib PREFIX=/usr/i686-w64-mingw32 export CC CXX CPPFLAGS CXXFLAGS CFLAGS LDFLAGS HOST ** For CLN ** ./configure --host=$HOST --enable-shared --prefix=$PREFIX && make && make install CLN_CFLAGS=$CPPFLAGS CLN_LIBS="$LDFLAGS -lcln" export CLN_CFLAGS CLN_LIBS ** For GiNaC** ./configure --host=$HOST --enable-shared --prefix=$PREFIX && make && make install Thanks Harish

2 1

Patch: Replace mktemp with mkstemp
by Till Hofmann 08 Oct '16

08 Oct '16

Hi all, while building 1.7.1 I noticed that mktemp is still used. I've attached a patch that replaces mktemp with mkstemp. Note that the behavior is slightly different: The file is opened twice because mkstemp directly opens the file, too. I'm not aware of any other solution to safely create an ofstream for a tmp file. If there is, please let me know. Kind regards, Till

2 1

removed function in 1.7.1 without soname bump
by Till Hofmann 06 Oct '16

06 Oct '16

Hi all, I'm a Fedora co-maintainer for ginac. We're currently working on updating ginac to the new release 1.7.1. During this process, we realized that there is a function that was removed in 1.7.1, but the soname was not changed. Was this intentional? Usually, removing functions requires an soname bump so other software doesn't accidentally link against the wrong version. The removed function is 'method virtual GiNaC::ex GiNaC::pointer_to_map_function_1arg<const GiNaC::ex&>::operator()(const GiNaC::ex&)' I've attached the output of abidiff to this email. For reference, here is the exact command to get the abidiff: $ abidiff --d1 ginac-old/usr/lib/debug/ --d2 ginac-new/usr/lib/debug ginac-old/usr/lib64/libginac.so.6.0.0 ginac-new/usr/lib64/libginac.so.6.0.1 Kind regards, Till

2 3