为此,我尝试编写以下代码来获取矩阵中每个元素的锁:
std::map<unsigned int, omp_lock_t> ds_lock;
void Sparse_Matrix_FL::setValue(unsigned int rid, unsigned int cid, double value)
{
omp_set_lock(&ds_lock(rid,cid));
std::map<unsigned int, double>::iterator str = data_Matrix[rid].begin();
std::map<unsigned int, double>::iterator end = data_Matrix[rid].end();
while (str->first != cid && str != end) str++;
if (str != end)
if (value != 0)
str->second = value;
else
data_Matrix[rid].erase(str);
else
if (value != 0)
data_Matrix[rid][cid] = value;
omp_unset_lock(&ds_lock(rid,cid));
}
但是我得到了错误:
在没有适当的operator()的情况下调用类类型的对象 将函数转换为指针到函数类型
这是另一个错误。
term的结果不适用于带有2个参数的函数
我不明白我在哪里做错了。
答案 0 :(得分:2)
您可以通过下标运算符std::map
或使用[]
方法来查找find
元素。 std::map
是不可调用的,因此您不能使用ds_lock(rid,cid)
语法。
您的地图需要使用std::pair< unsigned int, unsigned int >
作为其键类型。
您将可以使用omp_set_lock(&ds_lock[std::make_pair(rid,cid)]);
但是,您需要确保已预先初始化映射的每个使用的元素,否则由于std::map
线程不安全,当下标运算符自动创建新元素时,您会遇到问题。需要在新创建的值上调用omp_init_lock
。
建议确保始终使用RAII模式调用omp_unset_lock
。解决上述问题的一个可以正常工作的示例如下所示:
class omp_locker
{
public:
omp_locker( omp_lock_t& lock )
: lock( lock )
{
omp_set_lock( &lock );
}
omp_locker( const omp_locker& ) = delete;
omp_locker& operator=( const omp_locker& ) = delete;
~omp_locker()
{
omp_unset_lock( &lock );
}
private:
omp_lock_t& lock;
};
std::map< std::pair< unsigned int, unsigned int >, omp_lock_t> ds_lock;
std::mutex ds_lock_mutex;
omp_lock_t& getLock( unsigned int rid, unsigned int cid )
{
std::lock_guard lock( ds_lock_mutex );
auto key = std::make_pair( rid, cid );
auto it = ds_lock.find( key );
if ( it == ds_lock.end() )
{
it = ds_lock.insert( it, std::make_pair( key, omp_lock_t() ) );
omp_init_lock( it->second );
}
return it->second;
}
void Sparse_Matrix_FL::setValue(unsigned int rid, unsigned int cid, double value)
{
omp_locker lock( getLock( rid, cid ) );
std::map<unsigned int, double>::iterator str = data_Matrix[rid].begin();
std::map<unsigned int, double>::iterator end = data_Matrix[rid].end();
while (str->first != cid && str != end) str++;
if (str != end)
if (value != 0)
str->second = value;
else
data_Matrix[rid].erase(str);
else
if (value != 0)
data_Matrix[rid][cid] = value;
}
您需要在程序结束时在omp_destroy_lock
的每个元素上调用std::map
。