November 2015 - Boost-users

Using boost on AIX with GCC
by Chris Franz 10 Nov '15

10 Nov '15

I'm trying to get started using boost on AIX with GCC. What I seem to be running into is that GCC has a "feature" on AIX such that system headers are wrappered with an implicit 'extern "C"'. I'm not finding a good way to include the boost headers because of this feature. As an example, using about the simplest boost header I could find: $ echo "#include <boost/function_equal.hpp>" > fe.cpp $ g++ fe.cpp -c -isystem/usr/local/include/boost-1_59 In file included from fe.cpp:1:0: /usr/local/include/boost-1_59/boost/function_equal.hpp:15:1: error: template with C linkage template<typename F, typename G> ^ /usr/local/include/boost-1_59/boost/function_equal.hpp:22:1: error: template with C linkage template<typename F, typename G> ^ $ One obvious possibility is to use -I instead of -isystem, but not treating boost headers as system headers results in a multitude of warnings using the verbosity we like to use (I've even run into some boost headers that emit errors if -pedantic-errors is specified). Do I need a different version of gcc and/or boost? How do other people configure this?

2 1

Newbie not able to compile mini_c example
by Dezai G 09 Nov '15

09 Nov '15

Hi boost usersIm newbye in using boost spirit.I have recenly downloaded boost 1.59 library. Im tryng than to compile the mini_c example both on windows with visual studio 2008 and on linux with gcc 4.4.6Im walking throught plenty of errors and didnt found in the documentation what is going wrong.I ve tryed several things throuh search in google but nothing seems to walk An example of error is for linux/gcc : boost_1_59_0/boost/utility/result_of.hpp:190: error: wrong number of template arguments (1, should be 3) Are my compilers incompatible with boost 1.59. In this case what should i use And what are the correct version of library should i use to be fully compliant with visual studio 2à08 or gcc 4.4.6 Best regardsDZ

1 0

[BGL] Connected components of a distributed graph with few edges
by Quaglino Alessio 06 Nov '15

06 Nov '15

Hello, I am trying to compute the number of connected components of a graph with only few edges (i.e. many vertices do not have edges). This works OK for a non-distributed graph. The simplest example of the issue I could come up with is the following: typedef adjacency_list <vecS, vecS, undirectedS> SerialGraph; typedef adjacency_list<vecS, distributedS<mpi_process_group, vecS>, undirectedS> Graph; typedef iterator_property_map<std::vector<int>::iterator, property_map<Graph, vertex_index_t>::type> LocalMap; #ifdef PARALLEL_GRAPH Graph G(nV+1); synchronize(G); std::vector<int> localComponent(nV+1); LocalMap components(localComponent.begin(),get(vertex_index, G)); num = connected_components_ps(G, components); #else SerialGraph G(nV+1); std::vector<int> globalComponent(nV+1); num = connected_components(G, &globalComponent[0]); #endif The distributed version causes a seg fault when calling connected_components_ps, while the non-distributed works fine and correctly finds nV+1 components. If I populate the graph connecting all vertices with edges, both do work fine. What am I doing wrong? Thank you in advance. Regards, Alessio

2 8

BGL, Johnson all pair shortest path and C
by Paolo Bolzoni 05 Nov '15

05 Nov '15

Dear list, I have to use johnson all pair shortest path[1] to have the all pair distance matrix and pass it to a C function. Here is the problem, the C function expect a row-order array, while the Johnson expect a BasicMatrix[2] that understands the m[row][col] syntax. I solved it as shown in the bottom of the email, but it seems so convoluted, is there a simpler way I missed? Yours faithfully, Paolo [1] http://www.boost.org/doc/libs/release/libs/graph/doc/johnson_all_pairs_shor… [2] http://www.boost.org/doc/libs/release/libs/graph/doc/BasicMatrix.html /* -------------- */ //p is the graph //distance_matrix is a std vector auto size = boost::num_vertices(p); Matrix distance_matrix1(distance_matrix, size, size); MatrixR distance_matrix_wrapper(distance_matrix1); //pass distance_matrix_wrapper to johnson algorithm. /* -------------- */ struct Matrix { Matrix(std::vector<int64_t>& storage, size_t rows, size_t columns) : m_{storage}, columns_{columns} { m_.clear(); m_.resize(rows * columns_, int64_t{0}); } int64_t& get(size_t row, size_t column) { return m_[row * columns_ + column]; } std::vector<int64_t>& m_; size_t columns_; }; struct MatrixR; struct MatrixRR { MatrixRR(MatrixR const& m, size_t column) : m_{m}, column_{column} {} int64_t& operator[] (size_t row) const; MatrixR const& m_; size_t column_; }; struct MatrixR { MatrixR(Matrix& m) : m_{m} {} MatrixRR operator[] (size_t column) const { return MatrixRR(*this, column); } Matrix& m_; }; int64_t & MatrixRR :: operator[] (size_t row) const { return m_.m_.get(row, column_); }

1 0

Cannot get boost::asio::read_until to properly read a line. [new to boost/asio]
by Nick Chambers 04 Nov '15

04 Nov '15

Hello everyone! I'm writing an IRC framework (still in very rough, early stages), and am having trouble reading a single line from the server (in this case irc.freenode.net:6667). Ideally, I would get output that looks like the following: > :wolfe.freenode.net NOTICE * :*** Looking up your hostname... > :wolfe.freenode.net NOTICE * :*** Checking Ident > :wolfe.freenode.net NOTICE * :*** Your forward and reverse DNS do not > match, ignoring hostname > :wolfe.freenode.net NOTICE * :*** No Ident response > :wolfe.freenode.net 001 foobarbazball :Welcome to the freenode Internet > Relay Chat Network foobarbazball > :wolfe.freenode.net 002 foobarbazball :Your host is wolfe.freenode.net[ > 91.217.189.44/6667], running version ircd-seven-1.1.3 > :wolfe.freenode.net 003 foobarbazball :This server was created Sun Aug 2 > 2015 at 18:58:22 UTC > :wolfe.freenode.net 004 foobarbazball wolfe.freenode.net ircd-seven-1.1.3 > DOQRSZaghilopswz CFILMPQSbcefgijklmnopqrstvz bkloveqjfI > :wolfe.freenode.net 005 foobarbazball CHANTYPES=# EXCEPTS INVEX > CHANMODES=eIbq,k,flj,CFLMPQScgimnprstz CHANLIMIT=#:120 PREFIX=(ov)@+ > MAXLIST=bqeI:100 MODES=4 NETWORK=freenode KNOCK STATUSMSG=@+ CALLERID=g > :are supported by this server > :wolfe.freenode.net 005 foobarbazball CASEMAPPING=rfc1459 CHARSET=ascii > NICKLEN=16 CHANNELLEN=50 TOPICLEN=390 ETRACE CPRIVMSG CNOTICE DEAF=D > MONITOR=100 FNC > TARGMAX=NAMES:1,LIST:1,KICK:1,WHOIS:1,PRIVMSG:4,NOTICE:4,ACCEPT:,MONITOR: > :are supported by this server > :wolfe.freenode.net 005 foobarbazball EXTBAN=$,ajrxz WHOX CLIENTVER=3.0 > SAFELIST ELIST=CTU :are supported by this server > :wolfe.freenode.net 251 foobarbazball :There are 153 users and 93407 > invisible on 27 servers > :wolfe.freenode.net 252 foobarbazball 24 :IRC Operators online > :wolfe.freenode.net 253 foobarbazball 8 :unknown connection(s) > :wolfe.freenode.net 254 foobarbazball 54832 :channels formed > :wolfe.freenode.net 255 foobarbazball :I have 6701 clients and 1 servers > :wolfe.freenode.net 265 foobarbazball 6701 7241 :Current local users > 6701, max 7241 > :wolfe.freenode.net 266 foobarbazball 93560 97293 :Current global users > 93560, max 97293 > :wolfe.freenode.net 250 foobarbazball :Highest connection count: 7242 > (7241 clients) (668628 connections received) > :wolfe.freenode.net 375 foobarbazball :- wolfe.freenode.net Message of > the Day - > :wolfe.freenode.net 372 foobarbazball :- Welcome to wolfe.freenode.net in > Stockholm, SE. > :wolfe.freenode.net 372 foobarbazball :- Thanks to > http://www.portlane.com/ for sponsoring > :wolfe.freenode.net 372 foobarbazball :- this server! > :wolfe.freenode.net 372 foobarbazball :- > :wolfe.freenode.net 372 foobarbazball :- WOLFE, GENE [1931-]. Prolific > writer of short stories and > :wolfe.freenode.net 372 foobarbazball :- novels. His best-known work is > the multi-volume novel The > :wolfe.freenode.net 372 foobarbazball :- Book of the New Sun. He has won > multiple awards including > :wolfe.freenode.net 372 foobarbazball :- the Nebula Award, the World > Fantasy Award, The Campell > :wolfe.freenode.net 372 foobarbazball :- Memorial Award and the Locus > Award. He was awarded the World > :wolfe.freenode.net 372 foobarbazball :- Fantasy Award for lifetime > achievement in 1996. > :wolfe.freenode.net 372 foobarbazball :- > :wolfe.freenode.net 372 foobarbazball :- Welcome to freenode - supporting > the free and open source > :wolfe.freenode.net 372 foobarbazball :- software communities since 1998. > :wolfe.freenode.net 372 foobarbazball :- > :wolfe.freenode.net 372 foobarbazball :- By connecting to freenode you > indicate that you have read and > :wolfe.freenode.net 372 foobarbazball :- accept our policies as set out > on http://www.freenode.net > :wolfe.freenode.net 372 foobarbazball :- freenode runs an open proxy > scanner. Please join #freenode for > :wolfe.freenode.net 372 foobarbazball :- any network-related questions or > queries, where a number of > :wolfe.freenode.net 372 foobarbazball :- volunteer staff and helpful > users will be happy to assist you. > :wolfe.freenode.net 372 foobarbazball :- > :wolfe.freenode.net 372 foobarbazball :- You can meet us at FOSSCON ( > http://www.fosscon.org) where we get > :wolfe.freenode.net 372 foobarbazball :- together with like-minded FOSS > enthusiasts for talks and > :wolfe.freenode.net 372 foobarbazball :- real-life collaboration. > :wolfe.freenode.net 372 foobarbazball :- > :wolfe.freenode.net 372 foobarbazball :- We would like to thank Private > Internet Access > :wolfe.freenode.net 372 foobarbazball :- ( > https://www.privateinternetaccess.com/) and the other > :wolfe.freenode.net 372 foobarbazball :- organisations that help keep > freenode and our other projects > :wolfe.freenode.net 372 foobarbazball :- running for their sustained > support. > :wolfe.freenode.net 372 foobarbazball :- > :wolfe.freenode.net 372 foobarbazball :- In particular we would like to > thank the sponsor > :wolfe.freenode.net 372 foobarbazball :- of this server, details of which > can be found above. > :wolfe.freenode.net 372 foobarbazball :- > :wolfe.freenode.net 376 foobarbazball :End of /MOTD command. > :foobarbazball MODE foobarbazball :+i > My output on the other hand looks more like: > :sendak.freenode.net NOTICE * :*** Looking up your hostname... > :sendak.freenode.net NOTICE * :*** Checking Ident > :sendak.freenode.net NOTICE * :*** No Ident response > :sendak.freenode.net 001 squawker :Welcome to the freenode Internet Relay > Chat Network squawker > Ibq,k,flj,CFLMPQScgimnprstz CHANLIMIT=#:120 PREFIX=(ov)@+ MAXLIST=bqeI:100 > MODES=4 NETWORK=freenode KNOCK STATUSMSG=@+ CALLERID=g :are supported by > this server > upported by this server > 93546, max 95827 > freenode.net 372 squawker :- SENDAK, MAURICE [1928-2012]. The author of > Seven Little > 372 squawker :- Welcome to freenode - supporting the free and open source > wker :- any network-related questions or queries, where a number of > sendak.freenode.net 372 squawker :- We would like to thank Private > Internet Access > :- of this server, details of which can be found above. > main.cpp: > #include <IRCClient/client.hpp> > /* > #include <IRCClient/logger.hpp> > #include <IRCClient/plugin.hpp> > */ > > int main(int argc, char *argv[]) { > /* > if(argc > 2) { > logger::log(logger::output_type::STDERR, "usage: {1} > [config_file]", argv[0]); > return 1; > } > */ > client::irc bot( > "irc.freenode.net", // server to connect to > "6667", // port to connect on > "", // password > false, // don't connect with ssl > "squawker", // nick > "squawker", // realname > "A c++11 Utility bot written by nchambers", // ident > "#cplusplus.com", // channels to join > "." // command prefix > ); > > bot.on_connect([&bot]() -> void { > bot.nickserv(client::ns::identify, <REDACTED>); > }); > /* > bot.add_handler(client::input_type::privmsg, [](std::string channel, > std::string user, std::string msg) -> void { > logger::log(logger::output_type::STDOUT, "[{1}/{2}] {3}", user, > channel, msg); > }); > > bot.add_handler(client::input_type::command, "load", > [&bot](std::string channel, std::string user, std::string msg) -> void { > bot.load_plugin(plugin::plugin(msg)); > > if(bot.plugin_operation_successful()) { > logger::log(logger::output_type::STDOUT, "successfully loaded > plugin {1}", msg); > bot.msg(channel, "{1}: successfully loaded plugin {2}", user, > msg); > } > > else { > logger::log(logger::output_type::STDERR, "failure to load > plugin {1}", msg); > bot.msg(channel, "{1}: failure to load plugin {2}", user, msg); > } > }); > > bot.add_handler(client::input_type::command, "unload", > [&bot](std::string channel, std::string user, std::string msg) -> void { > bot.unload_plugin(plugin::plugin(msg)); > > if(bot.plugin_operation_successful()) { > logger::log(logger::output_type::STDOUT, "successfully > unloaded plugin {1}", msg); > bot.msg(channel, "{1}: successfully unloaded plugin {2}", > user, msg); > } > > else { > logger::log(logger::output_type::STDERR, "failure to unload > plugin {1}", msg); > bot.msg(channel, "{1}: failure to unload plugin {2}", user, > msg); > } > }); > > bot.add_handler(client::input_type::command, "reload", > [&bot](std::string channel, std::string user, std::string msg) -> void { > bot.reload_plugin(plugin::plugin(msg)); > > if(bot.plugin_operation_successful()) { > logger::log(logger::output_type::STDOUT, "successfully > reloaded plugin {1}", msg); > bot.msg(channel, "{1}: successfully reloaded plugin {2}", > user, msg); > } > > else { > logger::log(logger::output_type::STDERR, "failure to reload > plugin {1}", msg); > bot.msg(channel, "{1}: failure to reload plugin {2}", user, > msg); > } > });*/ > > bot.run(); > return 0; > } > client.hpp: > #ifndef _ROOT_SQUAWKER_IRCCLIENT_CLIENT_HPP__ > #define _ROOT_SQUAWKER_IRCCLIENT_CLIENT_HPP__ > > #include <string> > #include <functional> > > #include <boost/asio.hpp> > > namespace client { > enum class ns: int { > identify=0, > }; > > enum class input_type: int { > privmsg=0, > }; > > class irc { > public: > irc(std::string, std::string, std::string, bool, std::string, > std::string, std::string, std::string, std::string); > > void on_connect(std::function<void()>); > void nickserv(client::ns, std::string); > void run(); > > void send(std::string); > > private: > void connect(); > std::string readline(); > > std::string server; > std::string port; > std::string pass; > bool use_ssl; > std::string nick; > std::string realname; > std::string ident; > std::string channels; > std::string prefix; > std::function<void()> on_connect_callback; > boost::asio::io_service io_service; > boost::asio::ip::tcp::socket socket; > }; > } > > #endif > client.cpp: > #include <IRCClient/client.hpp> > #include <iostream> > #include <sstream> > #include <array> > > client::irc::irc( > std::string server, std::string port, std::string pass, > bool use_ssl, std::string nick, std::string realname, > std::string ident, std::string channels, std::string prefix > ) : > server(server), port(port), pass(pass), > use_ssl(use_ssl), nick(nick), realname(realname), > ident(ident), channels(channels), prefix(prefix), > on_connect_callback([]()->void{}), socket(this->io_service) { > } > > void client::irc::on_connect(std::function<void()> on_connect_callback) { > this->on_connect_callback = on_connect_callback; > } > > void client::irc::nickserv(client::ns command, std::string arguments) { > switch(command) { > case client::ns::identify: > this->send("nickserv identify " + arguments); > break; > } > } > > void client::irc::run() { > this->connect(); > > while(true) { > std::array<char, 128> buf; > boost::system::error_code error; > size_t len = this->socket.read_some(boost::asio::buffer(buf), > error); > > if(error == boost::asio::error::eof) { > break; > } > > else if (error) { > throw boost::system::system_error(error); > } > > std::cout.write(buf.data(), len); > } > } > > void client::irc::send(std::string msg) { > this->socket.write_some(boost::asio::buffer(msg + > std::string("\r\n"))); > } > > void client::irc::connect() { > boost::asio::ip::tcp::resolver resolver(this->io_service); > boost::asio::ip::tcp::resolver::query query(this->server, this->port); > boost::asio::ip::tcp::resolver::iterator endpoint_iterator = > resolver.resolve(query); > > boost::asio::connect(this->socket, endpoint_iterator); > this->send("NICK " + this->nick); > this->send("USER " + this->realname + " * * :" + this->ident); > > while(true) { > std::string line = this->readline(); > std::cout << line << std::endl; > > if(line == ":squawker MODE squawker :+i") { > break; > } > } > > this->on_connect_callback(); > } > > std::string client::irc::readline() { > boost::asio::streambuf buffer; > size_t received = boost::asio::read_until(this->socket, buffer, > "\r\n"); > std::istream is(&buffer); > std::string line; > std::getline(is, line); > return std::string(std::begin(line), std::end(line)-1); > } > any help would be very much appreciated! thank you in advance!

2 2

[BGL] Correct idiom for working with vertices.
by Matthew Markland 04 Nov '15

04 Nov '15

All: I'm building a graph on which I will be running betweenness centrality. This graph consists of a number of edges represented as pairs. For example 1 20 15 10 would represent an edge from 1<->20 and one from 15<->10. These are bidirectional. My initial code created a property_map for the centrality results and then walked through that map using num_vertices(g) as the limit for the for-loop. What I've realized now, though, is that when BGL generates the values for num_vertices (and the iterators for vertices also), it "fills in" the gaps. So, in the graph above there would be 20 vertices that would be interated through even though most of those have no relation to edges in the graph. So, I'm wondering what the usual idiom is in this situation. Should I instead iterator over edges and pull information from the source and target vertices of the edges? Thanks in advance for the information! Matt ---- Matthew Markland mwmarkland(a)outlook.com

2 1

Re: [Boost-users] Crashing Issue
by swati upadhyaya 04 Nov '15

04 Nov '15

hi I am using like this where RunningThreads_ is vector containing boost thread object pointer(handler).I am using interrupt() for the interruption of thread interruption_point() catching that interruption then joining.everytime its crashing in join() for me but for others who have the same configuration its working fine all of us has windos 8. For me also it was working last week but this week the same piece of code is crashing. :( RunningThreads_[i].handle->interrupt(); RunningThreads_[i].handle->join(); delete RunningThreads_[i].handle; On Thu, Feb 5, 2015 at 9:31 AM, swati upadhyaya <swatiupadhyaya(a)gmail.com> wrote: > RunningThreads_[i].handle->interrupt(); > LOG4CPLUS_INFO(root, "Thread Interrupted for " << > RunningThreads_[i].app_name << " successful"); > RunningThreads_[i].handle->join(); > LOG4CPLUS_INFO(root, "Waiting for thread to exit " << > RunningThreads_[i].app_name << " successful"); > delete RunningThreads_[i].handle; > LOG4CPLUS_INFO(root, "Thread handle for " << > RunningThreads_[i].app_name << " successfully deleted"); > //RunningThreads_.erase(i); > > On Wed, Feb 4, 2015 at 12:10 AM, Vicente J. Botet Escriba < > vicente.botet(a)wanadoo.fr> wrote: > >> Le 03/02/15 09:10, swati upadhyaya a écrit : >> >> In my code join() timed_joined() is crashing, I am using these API >> after interrupt(),Same piece of code is working for others but in my system >> its crashing, >> Is there some dependency problem I am using >> Boost version: 1.55.0 >> OS: Windows 8 >> IDE: Visual Studio 2010 >> >> Hi, >> >> it could help if you post the example and the backtrace. >> >> Which version/OS/compiler use others for which this works? >> >> Vicente >> >> _______________________________________________ >> Boost-users mailing list >> Boost-users(a)lists.boost.org >> http://lists.boost.org/mailman/listinfo.cgi/boost-users >> > >

3 2

[boost-test] print list of names of failing tests on failure
by Steve Lorimer 03 Nov '15

03 Nov '15

I don't believe there is currently a way to print only the list of tests which failed. --report_level=short only prints the numbers of passed and failing tests --report_level=detailed prints a very verbose list of all tests Is there a middle ground which will give me the names of tests which have failed, but not those which have passed? Something along the lines of: Test module "Master Test Suite" has failed with: 1 test case out of 9 failed 1 test case out of 9 aborted 1 assertion out of 30 failed Test case "foo_test" was aborted with: 1 assertion out of 1 failed TIA Steve

1 0

[boost::asio] [stream::ssl] stackful coroutine yield from SSL handshake in servername callback
by SamJam Liddicott 03 Nov '15

03 Nov '15

I'm testing a spoofing MITM SSL proxy based on boost::asio and boost::asio::stream::ssl, using boost 1.55 I register the tlsext_servername_callback, so that I can observe the called name during SSL handshake with the client. In the callback, I have the information needed to either get/create an appropriate SSL_CTX with certificates to use with the client, or start an SSL handshake with the server, to obtain the server certificate for re-signing and then put that back into the SSL_CTX to be passed to the client, before returning the tlsext_servername callback. It is clear that I must block the tlsext_servername callback from returning until I have the new SSL_CTX ready, which could take some time if I need to wait on a handshake with the real server. I can block return by performing a blocking handshake with the server, but then I consume one of the io_service threads merely blocking, and I don't want to waste them like that. I can successfully run it all as a new thread using blocking IO, in this form: ssl_connect(...) { boost::thread([...]() ({ sni_callback = [...](std::string server_name) { server.handshake(...); SSL_set_SSL_CTX(client.native_handle(), server.ctx)); }; client_socket.handshake(...); handler(); // really done by io_service.post }).detach(); } But I want to use stackful coroutines in the belief that thay are less costly and that I will be able to have more connections in progress (though I may have to link with BoringSSL). Using a coroutine spawn, I can convert the handshake with the server, like this: ssl_connect(...) { boost::asio::spawn(my_strand, [...](boost::asio::yield_context yield) { sni_callback = [...](std::string server_name) { server.async_handshake(..., yield[ec]); // <---- works nicely SSL_set_SSL_CTX(client.native_handle(), server.ctx)); }; client_socket.handshake(...); handler(); // really done by io_service.post }); } but the handshake with the client must remain blocking, otherwise, with: client_socket.async_handshake(..., yield[ec]); I get segfaults in the coroutine with unhelpful stack traces, mentioning only some SSL cleanup functions. However if the client connection hangs, then the client handshake can still hang and block a worker thread. My guess is that if the handshake with the client is not blocking, the SNI callback does not occur on the coroutine, and so I cannot properly block the return and possibly the use of yield in the SNI callback is invalid. If there is a solution that doesn't involve patching either libssl or stream::ssl, it will involve some asio_invoke_handler hooks. However, for that to work, the coroutine would have to be activated, but not to complete the async_handshake call. Instead it would be activated to handle the handshake read completion, leading to the parsing, SNI callback & return, further handshake write, and yield back to the io_service. I don't know if the coroutine/strand implementation can even do that -- be activated for work and not for return. So my first question is: is it likely possible to make coroutine do this? Looking in boost source detail/io.hpp for case engine::want_input_and_retry I notice: next_layer_.async_read_some( boost::asio::buffer(core_.input_buffer_), BOOST_ASIO_MOVE_CAST(io_op)(*this)); which seems to make no attempt to use asio_hanlder_invoke to allow it to be hooked, so that might need fixing. What sort of asio_handler_invoke hook would cause the handler to be invoked on the coroutine? Could strand.dispatch() do it? If it would, I'm not sure how to write such a hook either to work with argument dependent lookup, or to invoke in the coroutine. Otherwise I may as well stick to starting a scope-limited thread, and using blocking handshake operations. But I like to know the limits... Samjam

1 0

boost::math::cyl_bessel_j_zero: few NaN
by campagne 02 Nov '15

02 Nov '15

double order = (double)lcur + 0.5; unsigned int n_zeros; vector<double> zeros; boost::math::cyl_bessel_j_zero(order,1,n_zeros, back_inserter(zeros)); I'am running Boost 1.59 on Mac OS X Maverick for lcur = 0, I expect zeros filed with Pi, 2 Pi, 3Pi.... but the first two zeros I get NaN I get also other nan as for the 2nd zeros for lcur=1 Have you experienced such problem or I have a memory leak somewhere ? Thanks JE

2 1